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1 Introduction 


This textbook on fluid mechanics is aimed at engineering students at Bachelor’s degree level. 
In Chapter 2 the fundamentals of fluid mechanics, needed for the description and analysis 
of flows in technology, are presented. One-dimensional stream filament theory, along with 
the momentum integral and angular momentum integral equations, already present us with 
one method of constructing fluid technological appliances and experiments. For example, the 
dimensions of a machine can be determined quite accurately in one first step and predictions 
made about the flow losses which occur. 


However, these methods fail when machines are to be optimized, or when appliances are to 
be developed to operate under extreme conditions, such as silent operation, high efficiency, 
small dimensions, greatly damped oscillatory behavior, etc. In addition to this, even for 
most simple cases of application, simple fluid mechanical fundamentals cannot determine 
the operating behavior of a machine accurately enough. This requires extensive experiments 
to be carried out and these can be very costly and time consuming. 


At this point we refer to more detailed textbooks which progress systematically from the 
fundamental fluid mechanical equations and their methods of solution to the application of 
fluid mechanical software. 


The methods used in fluid mechanics are analytical, numerical and experimental. Although 
numerical methods are increasingly used to replace experimental methods, all three are still 
needed to solve fluid mechanical problems. This book is restricted to theoretical, that is 
analytical, methods. After working through this book, students are qualified to understand 
the basics and phenomena of fluid mechanics, and will have obtained an initial insight into 
fluid mechanical software, even if further studies into the fundamental equations of fluid 
mechanics and numerical methods are necessary for the application of this software. 


The composition of this book is partly very theoretical. In order not to lose sight of the 
connection to technical applications during the extensive derivations, we have selected flows 
past airplane wings, flows past vehicles and flows in pipes of process engineering systems as 
representative examples to illustrate the basics of one-dimensional fluid mechanics. 


We now introduce the student to the great variety of fluid mechanical applications, demon- 
strating that flows are all around us in our technological environment, by describing selected 
examples of flows in the following introductory sections. 


1.1 Flows in Technology 


Flows are responsible for most transport and mixing processes, such as those which occur 
in the transport of pollutants in our environment, in industrial processes, and in living 
organisms. The combustion of the limited fossil fuels currently produces the majority of 
electrical and heat energy. Optimizing the flows in these combustion processes serves to 
reduce oil and pollutant emissions. Flows are also of interest in the propulsion of airplanes, 
ships and vehicles, in pumping oil and gas through pipelines, and in the production and 
coating of materials. It is the flow of oxygen and carbon dioxide in organisms which permits 
life. Flows are of importance in the construction of low drag automobiles and airplanes; 
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in the development of booster rockets and space craft for transport to space stations; in 
energy and environmental technology; from process engineering to the simulation of entire 
production plants; in geophysics and astrophysics; and from meteorology and climatology to 
medicine, where innovations frequently go hand in hand with fluid mechanical optimization 
of artificial heart valves, hearts and artificial vessels. 


We now turn to examples of technical flows. Our environment is characterized by flow 
phenomena in many different ways. The optimization of flows leads to a decrease in the drag 
of airplanes, trains and vehicles and thus to savings in fuel. In drive assemblies it leads to an 
increase in efficiency and a reduction in polluting emissions. When materials are constructed 
from a melt, flow optimization determines the inner structure and thus the strength and 
load capacity of the material. In chemical production plants and pipelines, the optimization 
of flows reduces losses, and thus also the pump power needed to manufacture and transport 
liquids and gases. 


The development of airplanes and trains over the last decades is shown in Figure 1.1. The aim 
here has been to find low drag body shapes corresponding to the transport velocity, in order 
to keep the engine fuel or the electrical power of the driving motors as low as possible. The 
development of civil air travel began in the 1930’s with the legendary Ju 52. It transported 
17 passengers with a velocity of 250 km/h and was driven by three piston air motors. The 
endeavor to fly as fast as possible from one place to another led to the development of the 
jet engine, which currently makes it possible to fly at an altitude of 10 km, at a speed 
of 950 km/h. Jumbo jets can transport up to 450 passengers, with the next generation 
carrying up to 650 passengers. The first member of this aerodynamic new generation of 
airplane was the Boeing 707 (center Figure 1.1). The decisive aerodynamic discovery was 
the swept wing in the early 1940’s at the aerodynamic laboratory in Gottingen. This only 
allowed low drag flight at so-called transonic velocities. The Airbus A310 is a representative 
of the new generation of civil aircraft. Here the body of the plane is much larger, in order to 
transport as many passengers as possible. And still a considerable saving of fuel is attained 
compared to the Boeing 707. As well as the improved aerodynamics of transonic airfoils, it 
is the lighter materials and improved production engineering, together with new fan engines 
and the automated two-pilot cockpit, which have led to this saving in fuel and thus to the 
reduction in polluting emissions. Fan engines have a considerably larger diameter than jet 
engines. One part of the cold air condensed by the fan is guided past the hot propulsion jet 
as an air coating. This has the additional useful effect of reducing the noise emissions of jet 
engines drastically while simultaneously increasing their efficiency. 


A similar course of development has also been seen in trains over the past decades. Because 
the power required increases with the cube of the velocity and the drag with the square, 
at velocities over 100 km/h it is necessary to adapt the aerodynamic shaping suitably. 
Whereas in the case of conventional steam engines, the flow optimization merely shielded 
the driver’s cab from smoke using side-mounted wind deflectors, in the case of Inter-City 
(IC) trains, a better low drag design of the locomotive and the air deflectors close to the 
carriage wheels were realized. It is only in the case of Inter-City-Express (ICE) trains, 
where speeds of 330 km/h are reached, that a thorough aerodynamic shaping has been 
technically implemented, although even here the current collectors, for example, could use an 
aerodynamic covering. The fluid mechanical development of trains is also not yet completed. 
At present there are projects planned for travel in pipes at speeds of up to 500 km/h. 
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In the past, the aerodynamics of airplanes and trains were developed entirely in wind tun- 
nels. Figure 1.2 shows the wind tunnel model of the Airbus A340 in its take-off phase. 
Scales integrated into the mounting support of the model measure six components of the 
aerodynamic forces. Because the model, at rest in the wind tunnel, is in an air flow with 
a wind speed of about 300 km/h, the floor of the wind tunnel has to be moved with the 
same speed. These are very thorough experiments and are a key factor in determining the 
time needed to develop an airplane, a process which can take up to eight years, from the 
definition of requirements (flight speed, payload), through the design phase, to the introduc- 
tion of the product. These very long and expensive development processes have now been 
reduced considerably with the use of fluid mechanical simulation on large-scale computers. 
Flow simulation permits reasonably simple variations in the geometry and flow parameters, 
without the need to build new wind tunnel models each time. In future projects, both fluid 
mechanical software on large-scale computers and the wind tunnel will be tools of develop- 





Fig. 1.1 : Development of airplanes and trains 
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ment for design engineers. Wind tunnel experiments will increasingly be used to verify the 
computer results. 


The low drag aerodynamic shaping of an automobile was already resolved technically in 
1938. The Mercedes-Benz racing car built in 1937 to beat the world land-speed record is 
shown in Figure 1.3. The present-day usual drag coefficient cw (dimensionless drag force) 
is 0.365, but by lowering the driver in the racing car and covering the wheels, a so-called 
streamline body (Chapter 2.3.2) was achieved, with the drastic reduction in drag to a cw 
value of 0.17. Figure 1.4 demonstrates that the optimal obtainable aerodynamic value is 
0.15, indicating just how remarkable was the achievement of those Mercedes-Benz engineers. 
These findings on motor vehicles were really only taken into account in the 1980’s, after the 
oil crisis made fuel savings necessary. At present, the automobile industry has reached a 
compromise in the drag coefficient with a value of about 0.26. This permits a comfortable 
passenger area with the necessary all-round view. 


This introduction to technical examples can be supplemented with an anecdote. In 1952, 
the only wind tunnel for the aerodynamic development of motor vehicles in Germany which 
was operational and equipped with the appropriate scales was in Schlichting’s institute in 
Braunschweig. So it stood to reason that the neighboring Wolfsburg factory should use the 
Braunschweig wind tunnel to measure the aerodynamic drag of the Volkswagen VW 11 and 
VW Xə models. The wind tunnel results are shown in Figure 1.5. The prototype VW Xə 
showed a considerably better drag coefficient of 0.22, while the VW Beetle which was fi- 
nally produced had a very poor drag coefficient of 0.4. Schlichting was so annoyed by the 
lack of attention paid to his work, that he presented the results in Figure 1.5 at the next 
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Fig. 1.2 : Model of an Airbus A 340 in a wind tunnel 
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Fig. 1.3 : Mercedes-Benz W125 in a wind tunnel 


international conference, much to the infuriation of the company involved. 


A further technical example of building aerodynamics is shown in Figure 1.6. The Tacoma 
Narrows bridge, which has now become famous for its erroneous aerodynamic construction, 
extends 1.81 km over the narrows of Puget Sound in Washington State. On the 7th of 
November 1940, the wind was blowing with a speed of about 68 km/h at right-angles to 
the bridge. This caused a periodic separation of the air, known as a Kármán vortex street, 
to form at the opposite side of the bridge. Unfortunately, the eigenfrequency of the bridge 
was the same as the frequency of the periodic flow separation, so that mechanical eigen- 
oscillations were stimulated, leading to the bridge’s collapse. 


The optimization of flows is also important in the construction of internal combustion 
engines. The well-known cycle of an Otto engine is shown in Figure 1.7. The fuel-air 
mixture is sucked in through the open intake valve by the backward motion of the piston. 
In order to attain a mixture which is as homogeneous as possible, a rotating flow, called a 
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Fig. 1.5 : Measurement of automobile drag 
coefficients in the wind tunnel 





tumble, is superimposed on the flow. In the second step, the valve is closed and the fuel- 
air mixture is compressed, so that, after igniting the combustion, the expanding hot gas 
moves the piston for the mechanical propulsion downwards. When the combustion cycle is 
completed, the exhaust is expelled through the outlet valve (step 4). It might be expected 
after more than 100 years of development of combustion engines that the flow processes 


Fig. 1.6 : Aerodynamically excited oscil- 
lations of Tacoma Bridge 
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Suction Compression Combustion Exhaust 





Fig. 1.7 : Cycle of an Otto engine 


of intake, compression, combustion and exhaust expulsion would already be optimized; yet 
the need for an additional catalyst to reduce emission of pollutants shows that this is still 
not the case. For example, intensive efforts are being made to realize the direct injection of 
fuel, usual in diesel engines, in Otto engines too. This is expected to yield a fuel saving of 
about 10%, with a simultaneous increase in efficiency. Figure 1.8 shows a direct fuel-injection 
Otto engine. In the compression phase, the combustible gas mixture introduced through the 
injection valve is ignited directly at the spark plug via the bend in the cavity. However, the 
fluid mechanical problem of optimizing the combustion with respect to pollutant emissions 
remains. 


Flows with chemical reactions are very important for many applications in energy, chem- 
ical and combustion technology. Figure 1.9 shows an example of a flame, whose structure 
is due to a flow. Flames are characterized by a wide spectrum of time and length scales. 
Typical length scales of the flow vary from the size of the combustion chamber right down to 
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Fig. 1.8 : Otto engine with direct fuel in- 
jection 
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Fig. 1.9 : Flow structure of a flame 


the smallest vortex in which kinetic energy is dissipated. The chemical reactions occurring 
in combustion determine the wide spectrum of time scales. Knowledge of the interaction of 
these combustion time scales with the spatial structure of the flow is required to predict the 
pollutants in combustion processes. 


Flows are also involved in the production of materials such as in the microscopic struc- 
ture of solidified steel in Figure 1.10. We see convection cells which have formed at the 
solidification front as the hot molten metal becomes rigid. They are due to lift flows and 
determine the material properties of steel, such as its strength. 


In process engineering and chemical production plants (Figure 1.11), it is pipe flows 
in bends and junctions which introduce losses. If there are liquid separators in use, as in 
refineries, multi-phase flows with drops and bubbles have to be taken into account, as they 
can complicate matters considerably in optimizing the flow processes. 


There are many more examples of flow in technology. If you have followed the text so far, 
your interest will be awakened to move on to the basics and methods of fluid mechanics in 
the following chapters, so that you too will be able to solve fluid mechanical problems in 
engineering. 


We conclude this introduction by referring to additional literature. To supplement this 


Fig. 1.10 : Microscopic structure of solid- 
ified steel 
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Fig. 1.11 : Production plant in chemical 
process engineering 





text on fluid mechanics, we recommend PRANDTL—Guide of Fluid Mechanics 2001 for 
additional material on the fundamentals of fluid mechanics. This text describes further areas 
of fluid mechanics, such as aerodynamics, turbulent flows, fluid mechanical instabilities, flows 
with heat and material transfer, two phase flows and flows with chemical reactions, flows 
in the atmosphere and in the ocean, biological flows, as well as flow machinery. The fluid 
mechanical phenomena derived from technical problems are to be found in our textbook H. 
OERTEL jr., M. BOEHLE 1999. Further information on analytical and numerical methods 
is to be found in the textbooks H. OERTEL jr., E. LAURIEN 2002, H. OERTEL jr., J. 
DELFS 1996, H. OERTEL jr. 1994. Analytical descriptions of the fluid mechanics basics and 
methods are to be found in G. K. BATCHELOR, 1994, K. GERSTEN, H. HERWIG 1992, W. 
SCHNEIDER. 1978, J. H. SPURK 1997, F.M. WHITE 1974 and the technical application of 
boundary layer theory in H. SCHLICHTING, K. GERSTEN 1999. The mathematical basics 
may be found in more detail in the books by K. MEYBERG, P. VACHENAUER 1997,1998. 


1.2 Flow Regimes 


The first contact one makes with a flow might be at a water faucet. If you hold your finger 
in the stream of water, you will feel a force F which the flow exerts on the finger. This 
force is called the drag, and is experienced by a body in a flow. The drag is dependent on 
the geometry of the body in the flow, the properties of its surface, the flowing medium and 
the flow variables. The drag has a different value for a stream of gas than for the water 
stream we have considered. In order not to constantly have to distinguish between gases 
and liquids, we will use the generic term fluid. 
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The moving fluid is considered to be a continuum. This means that we neglect the molecular 
structure of the fluid, since the mean free path of the molecules is small compared to the 
characteristic macroscopic dimensions of the flow field. The characteristic physical quantities 
of the flow field of Figure 1.12, such as the velocity vector V with components in three 
directions u,v, w, the pressure p, the density p and the temperature T, are assumed to be 
continuous functions of position X = (x,y,z) and time t. 


Let the finger in the stream of water which we considered initially be replaced by a horizontal 
plate, as in Figure 1.12. The unperturbed flow woo points in the vertical direction and 
has the index oo. In order to describe the flow, the three scalar field quantities p, p and T, 
as well as the three components (u,v,w) of the vector velocity V have to be calculated as 
functions of the coordinates (x,y,z) and the time t: 


ulg, y, z, t) 
p(z, y, 2, t) 9 p(z, y, 2, t) 9 78; y, 2, t) 9 v(x, yY, 2, t) = v(z, yY, 2, t) (1.1) 
wz, yY, z, t) 


In order to calculate these six flow quantities, we have the continuum mechanical funda- 
mental equations of mass, momentum and energy conservation as well as the ther- 
modynamic equations of state at our disposal. 


In Figure 1.13, the example of a horizontal plate in a flow is considered further, in order 
to introduce some fundamental terms needed for the description of flows. The left-hand 
picture shows flow paths made visible with small aluminum particles. One particular point 
is noticeable in the middle of the plate. This is called the stagnation point, where the 
streamlines branch to the left and right. In the stagnation point of a flow field, the velocity 
vector V is equal to zero, and the pressure p is at a maximum. 


The right-hand illustration in Figure 1.13 is a sketch of the flow. The no-slip condition 
holds at the surface of the plate. Here again the velocity is equal to zero, but the pressure 
generally varies along the coordinate x. The velocity perpendicular to the plate at the 
position under consideration varies from the value zero to the constant velocity of the outer 
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Fig. 1.12 : Forces acting in a flow 
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Fig. 1.13 : Liquid jet at a horizontal plate 


flow. This presents us with an initial division into different flow regimes: the boundary 
layer flow and the outer flow. If we take into consideration the flow properties, such as 
the viscosity u (see Chap. 2.1), responsible for the friction in the flow, the boundary layer 
flow is the viscous part of the flow field, and the outer flow is the inviscid part. 


The origin of the internal friction is the inter-molecular interaction forces of the fluid. 
Whereas two elastic spheres (Figure 1.14) exchange momentum and energy completely and 
instantaneously, and thus have the infinitely large interaction force sketched in Figure 1.14, 
the interaction between molecules of a moving fluid is characterized by repulsive or attrac- 
tive interaction forces (see Chap. 2.1), depending on their distance apart r. These interaction 
forces between the molecules determine the transport properties of the fluid, such as the 
viscosity (friction), heat conduction (energy transport) and diffusion (mass transport). 


Before collision 


m, m, 


After collision 


Yi Be Fig. 1.14 : Collision between two spheres 


m; m’, (point mechanics) 
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Fig. 1.15 : Interaction force in the collision 
= ; of hard spheres 


For the different regimes of the moving fluid we have the corresponding fundamental equa- 
tions of continuum mechanics: mass, momentum and energy conservation. These are valid 
for both viscous boundary layers flows, as well as for inviscid outer flows. 


If we consider the flow quantities velocity V and density p, we are led to a completely different 
division of the flow regimes. We speak of an incompressible flow if the density p in the 
flow field is constant for a given temperature, e.g. in water flows. The flow is compressible 
if the density varies in the flow field, e.g. in air flows. If the velocity vector v is equal to zero, 
the medium at rest is called hydrostatic (p = const.) or aerostatic (p variable). Similarly 
regimes where the fluid is in motion are called hydrodynamic or aerodynamic. 


Material at rest l Flow 
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Fig. 1.16 : Division of flow regimes 
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Examples are shown in Figure 1.16. For example, hydrostatics treats the linear pressure 
distribution in a standing column of water, aerostatics the pressure and temperature (or 
density) distributions in the atmosphere at rest, hydrodynamics the water flow past a plate, 
and aerodynamics the flow past a wing. 


Flow Past a Wing 


Let us consider the first technical flow problem, the flow past the wing of an airplane. This 
was already introduced in Chapter 1.1. Figure 1.17 shows the wing of an Airbus A310. 
The flow is from the left of the wing, with Mach number M (ratio of the free stream 
velocity Us to the speed of sound a, velocity of propagation of small perturbations). 
The free stream here has a subsonic Mach number of Mæ œ 0.8. A further dimensionless 
characteristic number which describes the viscous boundary layer regime of the flow past 
a wing is the Reynolds number, Rey, which is made up of the free stream velocity ua, 
the depth of the wing L and the kinematic viscosity v (v = p/p): Rep = Ux + L/v. For civil 
airplanes it has a value of about Re, & 7: 10". 


In flights, the flow losses must be kept low, so that the ratio between lift and drag attains as 
high a value as possible. In order to achieve this, the engineer must be aware of the different 
flow phenomena, so that the computational methods can be applied appropriately. 


However, it is not only the flow past a wing at altitude which is of interest. The design 
of the wing must be such that it generates enough lift at take-off and landing, with the 
use of additional high lift flaps. In addition, the development of an airplane must take into 





Fig. 1.17 : Wing of a civil airplane 
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account how the fuselage and engines of the plane affect the flow past the wing, and what, 
for example, the best position is for the engines. 


Analytical and, to an even greater extent, numerical methods are used to answer all these 
questions. During the design phase, the engineer endeavors, with the help of a few wind 
tunnel experiments, to keep the development costs and time as low as possible. In addition, 
the optimization of an Airbus wing and the investigation of the lift and drag behavior 
at different angles of attack, for example, would be inconceivable without modern fluid 
mechanical methods. 


The flow regimes at a section of a wing are shown in Figure 1.18, both with the dimensionless 
pressure distribution indicated, as well as the flow made visible with small particles. For 
the discussion that follows, we need the dimensionless pressure coefficient cp, which is 
defined as follows: 
P— Po 
Ge eh (1.2) 
2 Po ' Us 


p is the pressure at an arbitrary position in the flow field, with Pæ, pœ and us denoting the 
pressure, density and velocity respectively in the free stream. Figure 1.18 shows the values 
of —c, on the wing, indicating the negative pressure on the upper side (suction side) and 
the over-pressure on the lower side (pressure side) compared to the free stream. The free 
stream with the velocity us is decelerated along the stagnation line. At the leading edge 
of the wing, the flow comes to rest, reaching its maximum pressure coefficient of c, (—Cy 
minimum). This point on the wing is called the stagnation point. 





Flow made visible with small particles 


Fig. 1.18 : Flow regimes and pressure distribution on a wing 
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Moving away from the stagnation point, the stagnation line branches towards the suction 
and pressure sides. We will first look at —c, along the suction side. The flow is greatly 
accelerated from the stagnation point along the upper side (the —c, value becomes larger) 
and on the front part of the wing it reaches supersonic velocities. Further downstream, the 
flow is decelerated through a jump in pressure, called a shock wave, to subsonic velocities 
again (there is a sharp drop in the —c, value). The flow is then further decelerated towards 
the trailing edge. 


On the pressure side, the flow is also accelerated away from the stagnation point. However 
the acceleration close to the nose is not as great as that on the suction side, so that no 
supersonic velocities occur on the entire pressure side. At about the middle of the wing, the 
flow is again decelerated, and the —c, value approaches the —c, value on the suction side 
downstream. At the trailing edge, the pressure coefficients of the pressure and suction sides 
are approximately the same. 


A thin boundary layer forms on the suction and pressure sides. The suction and pressure 
side boundary layers meet at the trailing edge and form the wake flow further downstream. 
Both the flow in the boundary layer and the flow in the wake are viscous. Outside these 
regions, the flow is essentially inviscid. 


The different properties of the different flow regions result in different equations to compute 
each flow. In boundary layer flows it is the boundary layer equations which hold, to good 
approximation. On the other hand, calculating the wake flow and the flow close to the trailing 
edge is more difficult. In these regions the Navier-Stokes equations have to be solved. The 
inviscid flow in the region in front of the shock may be tackled with comparatively less 
difficulty using the inviscid equation. 


Fig. 1.19 : Flow traces on the surface of a 
wing in a wind tunnel 
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Figure 1.19 shows traces of color indicating the flows on a wing in a wind tunnel experiment. 
It is seen that the color follows straight lines in a large area of the width of the wing. The 
statements made above are valid in these regions. Close to the body of the plane, the flow 
lines are no longer straight, and they form a vortex on the rear surface of the wing. In the 
following chapters, this concept will be used together with the idea of flow separation, 
which has a considerable influence on the behavior of an airplane. 


As we conclude the example of a wing, we consider the question of why the wing of a civil 
airplane, in contrast to that of a glider, is swept. This is because of the Mach number depen- 
dence of the dimensionless drag coefficient ¢,, at high Mach numbers of 0.8. We introduce 
the drag coefficient cy, with 


W 


2 


(1.3) 


Cy = 


where W is the drag and A the cross-sectional area of the wing. The drag increases greatly 
for transonic flows. Since the aim in civil aircraft is to fly as fast as possible (high Mach 
number), but with as low a drag as possible to keep the fuel consumption down, the wing 
is swept to about @ = 30°. The geometric relation 


M = Ms: cos(&) 


lowers the local Mach number at that section in the profile of the wing by the value cos(@), 
reducing the drag coefficient by the corresponding amount. This means a civil aircraft can 
fly at Mach number M. = 0.8 at an altitude of 10 km with a velocity of 950 km/h. 


Flow Past an Automobile 


One of the first steps in developing an automobile is to decide on its contour. This is 
determined more by the designer than by the aerodynamic engineer. The limits within 
which the engineer determines the possible variations in the contour (Figure 1.20) are small 
compared to the length or width of the automobile to be developed. 


Taking these guidelines into account, the aerodynamic engineer optimizes the contour so 
that the flow drag is as small as possible. In recent years vehicles have been developed 
whose drag coefficients cw are smaller than 0.3. 


However, minimizing the drag is certainly not the only job of the aerodynamic engineer at 
the design stage. At the same time, all the forces and moments which arise in an air flow 
have to be taken into account in optimizing the contour. In particular, the lift force, the 
side-wind force and the moment about the main axis of the vehicle are important, since 
these influence the driving stability to the greatest degree. In addition, the aerodynamic 
engineer has to ensure that the wind noises are minimal, that the windows do not become 
dirty during travel, that the side mirrors do not vibrate at high speeds, and so on. 


In contrast to the flow past the wing of an airplane, the flow past an automobile can 
be described as incompressible to good approximation. This is because the variations in 
density are small. As in the case of the wing, in Figure 1.21 we differentiate between different 
flow regimes: inviscid flow, boundary layer flow and viscous wake flow. 
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inviscid flow 


Fig. 1.20 : Flow past an automobile 





The pressure distribution has a stagnation point at the radiator. The flow is accelerated 
along the hood, leading to a pressure drop. At the windshield, the flow is again blocked, 
leading to an increase in pressure. After the pressure minimum has been reached on the 
roof, the flow is decelerated, with a pressure rise. Downstream of the trunk the boundary 
layer flow turns into the wake flow. 


The flow in the wind tunnel experiment is made visible with smoke, and this shows that 
downstream from the back of the automobile, a backflow region forms. This is shown as the 
dark region into which no flow indicators (white smoke particles) can penetrate. 


If we consider the flow under the automobile we see that the lower ground flow can be 
considered to be a gap flow whose upper boundary is rough. The average height of the 
roughness is about + 10 cm for an automobile. This means that the flow directly at the 
upper wall has to be assumed to be swirled. 


inviscid flow 
Boundary layer 
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wake flow 


Fig. 1.21 : Flow regions and pressure force on an automobile 
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Fig. 1.22 : Structure of the wake flow of a vehicle 


In order to avoid the flow being swirled like this, front spoilers are placed on many automo- 
biles in the lower region of the vehicle, in front of the entrance to the gap. The result of this 
is that a large part of the flow does not pass under the vehicle, where the swirling would 
lead to an increase in the drag. The savings made in this way are greater than the losses 
due to the extra drag of the spoiler. 


The pressure force on the upper contour is considerably smaller in most regions than that 
under the vehicle, and this pressure difference leads to a lift. In the aerodynamic design, one 
attempts to keep both the lift and the drag small. As already mentioned in the discussion 
of the flow below the automobile, in many cases of application, changes are made to the 
contour and spoilers are introduced, which divert the flow so that the vehicle experiences a 
reduction in lift. 


Let us come back to the wake flow downstream of the vehicle. In Figure 1.22 the structure 
of the flow in a cross-section through the vehicle is considered. We recognise the backflow 
region with a stagnation point on the surface of the vehicle, together with a saddle point 
downstream in the the flow field, where the flow lines branch into the wake flow and the 
backflow region. Completing the three-dimensional structure of the wake flow, as in the 
right-hand picture in Figure 1.22, we see that in the upper region of the roof of the trunk a 
so-called horseshoe vortex forms, seen when snow is driven together at the rear of a vehicle. 


Liquid-Steam Separator 


In chemical production plants quite different flow regimes have to be taken into consid- 
eration. If we take the example of a liquid-steam separator, we have to differentiate between 
liquid flows in pipes and bubble flows, drop flows and steam flows. Liquid-steam separators 
are to be found, for example, in refineries for the extraction of heavy hydrocarbons from 
crude oil gas. 
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Figure 1.23 shows a basic sketch of such a set-up. It consists of a throttle valve, the separator 
(demister), the steam compressor and the pump which produces a pressure drop and under- 
cooling of the liquid. 


An incompressible liquid is present in flow region 1. This is guided via a throttle valve 
into flow region 2, a two-phase flow (liquid and steam). By means of an adiabatic throttle 
process, an obstacle such as a shut-off valve or a test diaphragm is used to produce a 
pressure drop Ap = p, — pe with po < py, with the temperature T, or the enthalpy h of the 
liquid remaining constant (log(p)-h diagram). The two-phase flow in 2 is then directed to a 
separator or demister, where isobaric separation of liquid and steam takes place. 


After separation, a compressible steam flow is present in region 3, and an incompressible 
liquid flow in region 5. These processes can be discussed by means of a thermodynamic 
log(p)-h diagram. Here the pressure p is plotted on a logarithmic scale against the enthalpy 
h of the medium. 


A mixture of liquid and steam is present in the two-phase or wet steam regime. This regime 
is defined to the left by the lower limiting curve and to the right by the upper limiting 
curve. Both curves meet at the critical point K. The lower limiting curve is the connecting 
line of all points at the start of evaporation, and separates the liquid from the two-phase 
regime. The upper limiting curve is the connection of all points at the end of evaporation 
and separates the two-phase regime from the steam. As is usual in thermodynamics, the 
variable X denotes the amount of steam in the vapor, so that the lower limiting curve is 
denoted by X = 0 and the upper by X = 1. 


By choosing different pressure differences Ap in the throttle, different flow regions 2 in the 
two-phase regime can be attained. This is indicated on the log(p)-A diagram by the thin 
lines. The steam flow 2 is compressed, leading to a compressible gas flow in flow region 4. In 
general, steam is any gas-like substance close to the limiting curve. The compressor ensures 
that the thermodynamical state of the medium in region 4 in the log(p)-h diagram is at a 
greater distance from the limiting curve than the steam in region 3. That is why the flow 
downstream from the compressor is called a gas flow rather than a steam flow. 


Inside the pump, where the incompressible liquid flow passes over from flow region 5 to flow 
region 6, the accelerating moving liquid can produce large pressure drops. Therefore it has 
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Fig. 1.23 : Sketch and pressure-enthalpy diagram of a liquid-steam separator 
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Fig. 1.24 : Flow regions in a two-phase 
flow 


to be ensured that the minimum static pressure Pin does not fall below the steam pressure 
pp. This would lead to damaging cavitation effects which could destroy the pump. 


The flow regime 1 as well as the regions 5 and 6 are typical examples of the incompressible 
flows through straight or curved pipes which frequently appear in technical applications. As 
a fluid flows through such pipes, friction effects lead to pressure losses Apy which have to be 
determined. For example, knowledge of these pressure losses is necessary in order to choose 
suitable pumps with appropriate efficiency. The flow region 2 of the two-phase flow, where 
the fluid is present in two aggregate states, is a typical example from process engineering. 
Flow regions 3 and 4 are examples of compressible flows, where, as well as the variations in 
pressure and velocity, the variation in the state variables density p and temperature T also 
have to be taken into account. 


Figure 1.24 shows the flow regimes of the two-phase flow. In the different pipe systems of 
the liquid-steam separator we distinguish between the liquid flow with X = 0, the bubble 
flow in in liquid with 0 < X << 1, the drop flow with fluid film at the walls of the pipe 
with 0 << X < 1, and the steam flow with X = 1. 


1.3 Product Development 


Whereas the aim of a scientist is to determine the mathematical and physical descriptions 
of flow processes, and thus derive the continuum mechanical fundamental equations, the 
engineer must implement this scientific knowledge in new products. The fundamentals of 
fluid mechanics described in Chapter 2 are necessary to do this, but so too are the analytical 
and numerical methods of solution and the fluid mechanical software of Chapter 3 needed 
to solve the fundamental equations. The formal development process of a new product is 
therefore always the same. 
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If we follow the development of a wing of an airplane, we first have to set out the require- 
ments. The airplane has to convey, for example, 250 passengers at an altitude of 10 km over 
a distance of 7,500 km at Mach number M% = 0.8. The requirements which, for example, 
the Airbus A300 satisfies lead to the following lift coefficient required for transport of the 
payload: 


A 


5 ‘ Poo * US Ar oe 


Ca = 


with lift force A, free stream velocity us and wing surface area Ap. It also has to be taken 
into consideration that the airplane will take off and land with a speed of about 250 km/h. 
This necessitates the integration of high-lift flaps into the wings, to increase the surface area 
of the wing at take-off and landing in correspondence to the lower velocity, so as to reach 
the required lift. 


The development engineer must design a wing which will have as low a drag as possible, 
in order to minimize the engine’s fuel consumption. The flight Mach number of M,, = 0.8 
leads to a sweep angle of @ = 30°. In addition, the integration of the wing into the cylindrical 
body of the airplane must be carried out, and the number of engines necessary for propulsion 
determined. 


The first job for the development engineer is the preliminary design, as shown in Figure 
1.25. The curvature of the wing profile is provisionally determined, using simple methods 
which we will meet in Chapter 2. The pressure distribution c, sketched in Figure 1.26 
has to satisfy the required lift coefficient c, and simultaneously have a drag coefficient cw 
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which is as low as possible. In the second step, the calculation of the swept wing designed 
with this profile is carried out. At this stage, the warping of the wing and its integration 
into the body of the plane also have to be taken into account. To do this we require the 
fluid mechanical software as presented in Chapter 3, as well as the simplifications of the 
fundamental equations in the different flow regions. These equations form a system of partial 
differential equations which have to be solved approximately using numerical methods. The 
first calculation of the wing generally will not attain the required lift coefficient ca, or the 
computed drag coefficient ¢,, will still be too large. A further iteration step is then necessary 
to get an improved preliminary design with the computed data. This design iteration is 
carried out in 2 to 3 steps. 





Preliminary design 


Calculation 


Wind tunnel experiment 





Fig. 1.26 : Wing design: preliminary design, calculation, wind tunnel experiment, verifica- 
tion and validation 
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If the required aerodynamic coefficients are satisfied, the second step of the design process 
takes place, the construction and the building of the wind tunnel model. This generally 
consists of stainless steel, and has numerous pressure holes to allow the measurement of the 
pressure distribution at different cross-sections of the wing. 


In Chapter 2 we will learn that the integral of this pressure distribution permits the cal- 
culation of the required forces. As shown in Figure 1.26, the wing model is fitted out with 
model engines, so that the airplane model in the wind tunnel is geometrically similar to the 
original. 


Extensive measurements are now carried out in different wind tunnels. The wing is in- 
vestigated at the free stream Mach number M, = 0.8 (950 km/h, 10 km altitude) and 
at different angles of attack, as well as in the take-off and landing phases with extended 
landing flaps, that is, with a larger wing surface and a lower speed of 250 to 300 km/h. 


The results of the calculation may not initially agree with the wind tunnel results. This is 
because of the mathematical and physical difficulty of the numerical solutions, and because 
of experimental errors and disturbances in the wind tunnel. 


The next important step in the design process is verification and validation. These require 
all the engineering skills of the developer. Verification is the comparison of the experimental 
results with the numerical results, as well as the adaptation of the numerical methods of 
solution and the instrumentation in the wind tunnel. Validation requires the further devel- 
opment of the physical models in the fundamental equations of the different flow regimes. 
This is a very time-consuming process which has a great influence on the development time 
of an airplane. 


As shown in Figure 1.26, in the verification and validation phases, the calculation, or the 
preliminary design, are corrected in 3 to 4 iteration steps, until the initial requirements are 
satisfied. In each iteration step, a new wind tunnel model has to be built, and the time- 
consuming measurements in the wind tunnels repeated. The fewer the number of iteration 
steps which have to be carried out, the more successful the design process. 


We have described the development steps of an airplane using the example of the aerody- 
namics of the wing. Similar development processes are also carried out for the structure 
of the plane, the development of the engines, the flight mechanics and the system integra- 
tion in the cockpit. Therefore different disciplines are interconnected in the different design 
processes. If we consider how many different companies in different countries are involved 
in the European project Airbus (Figure 1.27), for example, we can see how complex the 
development of an airplane is. The same is true of the development of an automobile, a flow 
engine or a process engineering production plant. For each of these, the development pro- 
cess described, using the example of a wing, has to be carried out for each of the disciplines 
involved. 


At the end of the development process comes the verification of the finished product. In 
the case of an airplane, this is its first flight and consequent registration. The airplane has 
to prove itself in certain extreme flight conditions. Only minimal changes to the airplane 
can now be carried out. Any development errors which were discovered during the design 
cycle can now only be corrected up to a point. For example, when developing the Airbus 
A320, the interaction between the relatively large engines and the body of the plane was 
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Fig. 1.27 : Construction elements of an Airbus A320 


not given enough attention. This manifests itself in an unpleasant drone in the body of the 
plane during take-off. The experiences gained in the final verification are retained in the 
databases of the airplane manufacturers and can be taken into account right at the start of 
the preliminary design of the next project. 
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2 Fundamentals of Fluid Mechanics 


2.1 Flow Properties 


We distinguish between the kinematic properties of the flow, and the transport proper- 
ties and thermodynamic properties of the fluid. While the kinematic properties velocity 
V, angular velocity w, acceleration b, and vortex strength WR are properties of the flow field 
and not of the fluid itself, we will see in Chapter 2.3.1 that the transport properties viscosity, 
heat conduction and mass diffusion, as well as the thermodynamic properties pressure p, 
density p, temperature T, enthalpy h, entropy s, the specific heats cp, €v, and the expansion 
coefficient a are properties of the fluid. First of all we present the definitions of the funda- 
mental concepts. For further details we recommend PRANDTL-Guide of Fluid Mechanics 
2001, and BIRD, STEWART, LIGHTFOOT 1960. 


2.1.1 Transport Properties 


One transport property which we have already met is the friction. It determines the mo- 
mentum transport in the viscous flow regimes and is linked to the gradient of the velocity 
vector V. For example, heavy oil or tar need a long time to flow out of a container, whereas 
light oil flows faster. 


In order to introduce the shear stress 7, we look at the one-dimensional flow problem in 
Figure 2.1. Between a lower plate at rest and an upper plate moving with constant velocity 
U, there is a constant shear rate with linear velocity profile u(z). This is called Couette 
flow. The boundary condition on the surfaces of the plates is the no-slip condition, which 
on the lower plate leads to uw = O and on the upper plate u = U. In order to maintain 
the constant velocity U in the presence of the friction, a constant force F is required. The 
applied force is proportional to the drag velocity F ~ U, proportional to the plate surface 
area A, F ~ A and inversely proportional to the height of the gap H, F ~ 1/H. 


Thus the force is 


Fig. 2.1 : Couette flow, definition of shear 
stress T 
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or with the constant of proportionality ps 


U-A 
F=u:- — 
LOH 
u is a material constant of the fluid which is called the dynamic viscosity. It has the dimen- 
sions [F - T/L?] with the force F, the characteristic time T and the characteristic length L, 
in our example the gap height H, and the units {N s/m?}. 


The shear stress 7 (shear rate) is 


F U 
T = — = . — 


A H 
with the units {N /ma}. For Couette flow we have the linear velocity profile 


U du 


H dz 


du 


If this linear relationship between the shear stress 7 and the velocity gradient du/dz holds, 
the fluid is a Newtonian fluid. Examples of Newtonian media are water, smoothly flowing 
oil and gases. 


This yields 


Using the relationship already derived, we can already discuss an important technical ap- 
plication. Air bearings have a particularly small friction drag. If, for example, we move a 
glass plate on a 0.1 mm thick air cushion with a constant velocity of 0.1 m/s, we have 


du U 1 
— = = 10? 
dz H 15} 


for air at 20° C, w= 1.71-10-° Ns/m?, and so 


d N 
T= uT — 1.71- 107? {aa} 


With a plate surface area of A = 0.01 m? we obtain the very small force 


F=7r-A=1.71-10*{N} 


In general flows are not, as assumed until now, one-dimensional. There are three shear stress 
components for each spatial direction, that is 9 components of the shear stress tensor 7;;. 
These characterize the friction in the three-dimensional flow field. In this terminology, the 
shear stress component Ty, in Couette flow determines the velocity profile u(z), with x the 
direction of flow and z the vertical coordinate 


Ou 1 Ow 
Tez — a a 
H Oz or 
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For the one-dimensional theory, to be treated in Chapter 2.3, with w and Ow/0x equal to 
zero, it suffices to set 


du 
dz 


Tez ZT =H 


As well as Newtonian fluids, we also have non-Newtonian fluids. In this case the functional 
relation in equation (2.1) is non-linear. Some examples of non-Newtonian fluids are shown 
in Figure 2.2. The curves for fluids which cannot withstand a shear rate must go through 
the origin. So-called yielding fluids have a finite shear stress even for vanishing velocity 
gradients. These fluids behave partly as rigid bodies and partly as fluids. The slope of the 
curves for pseudo-elastic fluids such as molten metal or high polymers decrease as the shear 
stress increases. In contrast to this, dilatant fluids, such as suspensions, have an increase in 
the slope. Toothpaste and mortar are examples of ideal Bingham materials. The finite value 
of rat du/dz = 0 follows the linear progression of a Newtonian fluid. In addition to this, the 
shear stress of some non-Newtonian media is time dependent. Even if the shear rate is kept 
constant, the shear stress changes. An ansatz frequently used for non-Newtonian media is 


Taz MKT, (2.2) 


where K and n are material constants. For n < 1 we have a pseudo-elastic fluid, n = 1 with 
K = mis a Newtonian fluid, and for n > 1 a dilatant fluid. Note that the ansatz (2.2) yields 
unrealistic values for the root in Figure 2.2. 


Numerous other laws have been derived for non-Newtonian media, mainly from experimen- 
tal results. In what follows we will no longer consider these and will restrict ourselves to 
Newtonian fluids. 


The viscosity uw of a Newtonian fluid is directly related to the interaction forces between 
the molecules of the flowing medium. If we consider the interaction force between two 
molecules in the air (nitrogen, oxygen) in Figure 2.3, we see that at large distances r the 
molecules have a negative attraction, and at small distances they have a positive repulsion. 
The attractive force between the molecules results from the Van der Waals interaction, whose 
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Fig. 2.2 : Shear stress 7 for Newtonian and 
non-Newtonian fluids 
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origin lies in the dipole moment produced by the deformation of the electron shells. The 
almost exponential repulsion is due to the electrostatic repulsion of shells with like charge. 
In a collision between molecules due to their eigen-motion, this interaction will first act to 
attract the molecules and then to repel them strongly. This interaction of the 10?° molecules 
per mole in the air causes not only friction, but also heat conduction and diffusion. Since 
the eigen-motion of the molecules, and thus also their collision probability, depends on the 
temperature T and pressure p, the viscosity u is also temperature and pressure dependent. 


Figure 2.4 shows the qualitative temperature dependence for liquids and gases at constant 
pressure. In liquids, the kinematic viscosity u decreases with increasing temperature, whereas 
in gases it increases. The viscosity of liquids and gases increases with increasing pressure. 


Following this brief digression into molecular physics, we look at the reason for the boundary 
conditions at fixed walls. The no-slip conditions V = 0 is due to the fact that the interaction 
between the molecules of the fluid and the crystal lattice at the fixed surface is considerably 
larger than that between the fluid molecules themselves. In continuum mechanical condi- 
tions, to which we will restrict ourselves in this textbook, each fluid molecule sticks to a 
fixed wall on collision. 


The energy transport through heat conduction can be developed analogously to the fric- 
tion. In Figure 2.5 the linear temperature profile T(z) in a fluid layer at rest between two 
horizontal plates with the temperature T} and To corresponds to the linear velocity profile 
u(z) of Couette flow. The heat flux q, which is the amount of heat transferred per unit time 
Q) per area A, corresponds to the shear stress T. 


F du Q dT 
SS a =—=-k-— , 2.3 
zer dz >) 
The heat flux q can be written as the Fourier law 
q=-k:- JT . (2.4) 


For the one-dimensional case under consideration, dT’/dz corresponds to the velocity gra- 
dient du/dz. This analogy is only valid for the one-dimensional case. For three-dimensional 
flow, we have already stated that the shear stress 7; is a tensor with 9 components, while 
q is a vector. 
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Introducing the kinematic viscosity with 


e 


the thermal diffusivity a = k/(p-c,), which has the same dimensions as v, can be used to 
form a dimensionless characteristic number 


V 
PSS 
a 


The Prandtl number Pr describes the ratio of momentum transport (friction) to energy 
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Fig. 2.5 : Analogy between friction and heat conduction 
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transport (heat conduction) in the fluid under consideration. Gases have the Prandtl number 
0.71, water 6.7, oils several thousand. 


The mass diffusion (mass transport) in the fluid may be treated equivalently if two media 
with partial densities p; (i = 1,2) mix together due to a concentration gradient. The con- 
centrations of the two components are C; = p;/p, with the total density of the mixture p. 
In analogy to friction and heat conduction, we postulate that the mass flux per unit time 
m, for the substance i is written as 


Mui 


A 
with the diffusion coefficient D[L?/T|{m?/s}. For mass concentrations C;, Fick’s law is 
written 


Mm; 


FG =P Ve Gi) 


In analogy to the Prandtl number, the dimensionless Schmidt number Sc and the Lewis 
number Le for the mass diffusion can be defined as follows: 


The Schmidt number describes the ratio of momentum transport to mass diffusion, the 
Lewis number the ratio of mass diffusion to energy transport. 


2.1.2 Thermodynamic Properties 


Classical thermodynamics, as taught in freshman courses, cannot be directly applied to 
fluid mechanics, since a viscous flow is not in thermodynamic equilibrium. However, in most 
technical applications, the deviation from local thermodynamic equilibrium is so small 
that it may be neglected. There are two exceptions: flows with chemical reactions and sudden 
changes in the thermodynamic variables of state, such as can occur in strong shock waves, 
to be treated in Chapter 2.3.3. 


The most important thermodynamic quantities are pressure p, density p, temperature T, 
entropy s, enthalpy h and the internal energy e. Two of these six quantities suffice to 
determine a thermodynamic state uniquely, as long as these are thermodynamic variables 
of state. The most important relations which we will need in the following chapters will now 
be discussed briefly. 


The first law of thermodynamics is written 
dE = dQ +d4W , (2.5) 


with dE the total energy of the system under consideration, dQ the heat supplied and dW 
the work done by the system. For a fluid at rest 


dW =—p-dV , dQ=T-ds , 
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where V is the volume. Thus (2.5) per unit mass becomes 
_ P 
de =T -ds + — -dp . (2.6) 
p 


Writing the total differential describing the change in the internal energy as 


de de 


de = — :ds + — -d 
e Ds ny, p , 

we obtain 

öde _ 2 Oe 

— Os p ’ P=P Op sS 
The enthalpy is by definition 

h=e+? , 
p 


and so with (2.6) the first law of thermodynamics can be written in the form 


1 
ah=T-ds+— dp) (2.7) 


The temperature T and 1/p are then 


= ðh 1 öh 


r=2 Z2, 
ðs? mr 


The thermal equation of state for the ideal gas is 


(2.8) 


with the gas constant R. This yields the speed of sound a 


o 
= L| =k RT , (2.9) 
Op 
with the dimensionless ratio of the specific heats ø 
Cp Oh de 
_ “p — „= — 2.10 
K Cy ’ Cp ƏT I» ’ C ƏT lo ( ) 
For flows with heat transport, the thermal expansion coefficient a is required: 
1 Op 
=—--—| ., 2.11 
a= -5 gT. (2.11) 


For ideal gases this yields 


Spe 


32 2. Fundamentals of Fluid Mechanics 


Liquids generally have thermal expansion coefficients which are smaller than 1/T. Negative 
values can also occur, such as in the case of water close to the freezing point. The thermal 
expansion coefficient can be used to write the dependence of the enthalpy on the pressure 
as 


d 
mege tasan (2.12) 
p 


For an ideal gas, the second term vanishes and the enthalpy depends only on the temperature 
h= Ty: 


In multi-phase flows, as met first in Chapter 1.2, the thermodynamic states close to the 
critical point are of particular interest. Figure 2.6 shows the isothermal lines of a liquid- 
steam mixture. The different thermodynamic behavior in the liquid and in the steam phase 
can be explained by considering the different inter-molecular interaction forces. If a gas is 
compressed isothermally, the mean translational energy of the molecules remains constant 
and the mean spacing between neighboring molecules decreases. If the specific volume 1/p 
of the gas becomes so small that the mean spacing is only a few molecular diameters, the 
attractive forces between the molecules become significant. If the temperature exceeds a 
critical value T}, a further reduction in the specific volume leads to an unstable state where 
the molecules are beyond the region of the inter-molecular attractive interaction forces and 
so begin to form molecule clusters. 


This intermediate state between the liquid and gas phases in unstable to the smallest dis- 
turbance. A minimal increase in the pressure will lead to complete condensation to a ho- 
mogeneous liquid with consequently large density, or a small drop in pressure will cause the 
homogeneous steam phase to occur with correspondingly small density. The almost constant 
pressure in the transition phase is called the saturation steam pressure. 


At temperatures above the critical temperature T}, the translational energy of the molecules 
becomes so large that the formation of molecule clusters is prevented. A continuous tran- 
sition occurs along the isothermal lines from the gas phase to the liquid phase, while the 
specific volume decreases. In this temperature regime, the Van der Waals equation describes 
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Fig. 2.6 : Isothermal lines of a liquid-steam 
1/p mixture 
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the thermodynamic state of real gases: 


R-p-T > 
— ie. 2.13 


where b and c are constants along the isothermal lines which characterize the attractive 
force between the molecules. 


A further property of fluids is the surface tension o of liquids and the interface tension 
between different liquids, or between liquids and solids. The temperature dependence of the 
surface tension can cause flows. Surface and interface tension is also due to the interaction 
forces between the molecules. Figure 2.7 shows the forces of a molecule in a liquid and 
of a molecule at the interface between liquid and gas. Within the liquid, the forces on 
the molecule under consideration balance each other out on average since the molecule is 
surrounded by the same number of molecules on all sides. At the surface of the liquid, the 
interaction between the liquid and gas molecules is considerably smaller than that between 
the liquid molecules. This causes the resulting force R, which causes the surface tension a. 
By definition this is 
|F| 

o=- (2.14) 
with the surface force F and the length of the surface L. For example, for the interface 
between water and air o = 7.1 -107° N/m at a given temperature. 


At a doubly curved surface with radii of curvature Rı and Rə, the balance of the forces at 
the surfaces leads to a pressure jump 


1 1 
Ap = — + — . 2.15 
p=o( iti] (2.15) 


This results in a higher pressure on the concave side of the curved surface. For a bubble or 
drop, with Ra = Rə = r, the pressure difference across the surface is 
2-0 


Ap 
7 


A soap bubble with an inner and an outer surface has the following increased pressure inside: 


4. 
r 





Fig. 2.7 : Surface tension 


34 2. Fundamentals of Fluid Mechanics 


This pressure difference in a drop causes, for example, a hole in a solid surface to be filled 
with liquid. The hole is then only filled when the angle of contact a between the liquid 
and the surface is smaller than 90°. 


This angle of contact between the liquid and solid surface is determined by the energy of 
the interfaces. It causes the raising and lowering of a liquid in a capillary. 


Let us consider the interfaces between different liquids in Figure 2.8. For example, on a glass 
surface a mercury droplet will not cause any wetting. The angle of contact a is larger than 
90° (about 150°) and the surface tension o of the mercury is larger than the adhesion force 
between mercury and glass. For a water droplet on a glass surface, the contact angle a is 
smaller than 90° and so there is wetting. The surface tension o of water is smaller than the 
adhesion force between water and glass. On the other hand, for a water droplet on a wax 
surface, there is no wetting and the contact angle is thus greater than 90°. The wetting of 
oil on glass is almost complete, so a — 0. The surface tension of the oil is vanishingly small 
compared to the adhesion force between oil and glass. 


The angle of contact œ between solid surfaces, liquids and gases can be computed using 
Young’s equation 


O solid/gas — 9 solid/liquid T O gas/liquid ' cos(a) ; (2.16) 
as long as the individual surface tensions are known. 


Because of the surface tension, the liquid tends to form minimal areas. This is seen in the 
experiment shown in Figure 2.9. A thread with a loop is laid onto a surface of soap suds. If 
the soapy skin inside the loop is burst, a circle forms, so that the remaining liquid surface 
forms a minimal area. Gradients of the surface tension Yo cause shear forces 7 in the 
surrounding media A and B, as, for example, at the interface between a liquid and a gas, 


VO =TA+TB 


The surface moves in the direction of greater surface tension, and by means of the shear 
stresses T4 and Tz, causes flows in the different media. Gradients in the surface tension can 
be caused by concentration gradients along the surface. This causes, for example, pieces 
of camphor to move back and forth sporadically on a water surface, since the camphor 
molecules locally reduce the surface tension. Another example are the tears in wine or 
cocktail glasses. Because of the concentration gradients in the water-alcohol mixture, the 
liquid at the glass is raised up and then flows back down into the liquid as regular drops. The 
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Fig. 2.8 : Angle of contact between solid, mercury, water, oil and air 
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Fig. 2.9 : Minimal areas 


evaporation of the alcohol causes a reduction of the alcohol content and thus an increase 
in the surface tension. The liquid is continually transported from the middle of the glass to 
the edge. 


Temperature gradients also cause gradients in the surface tension, bringing us back again 
to the thermodynamic properties of fluids. If a thin metal plate wetted with silicon oil is 
heated from below with a hot rod, a hole in the oil film occurs at the heated point. The 
increase in the temperature leads to a decrease in the surface tension. The liquid surface 
moves in the direction of the colder areas with greater surface tension. A piece of ice on the 
oil surface has the opposite effect. The liquid film causes a swelling in the colder region. 


The same effect can be used to transport bubbles in a liquid which is heated from one side. 
The cold side of the bubble has a higher surface tension than the warm side, and so it pulls 
the surface from the warm side, setting the bubble into motion. 


2.2 Hydro-Aerostatics 


In this chapter we will treat the properties and fundamental equations of fluids at rest. 
As in Chapter 1.2, hydrostatics describes incompressible fluids at rest and aerostatics 
compressible fluids at rest. When a fluid is at rest no shear stresses occur, so that the forces 
acting at each point on any volume element are normal to each surface. Therefore these 
forces can only be pressure or tractive forces. In a liquid, only pressure forces occur. Here 
pressure force F, on a surface element A increases with the size of the surface element. We 
introduce the scalar quantity pressure p as the pressure force per unit area. 

p= | with dimensions >| and units {ao} (2.17) 
The objective of hydro-acrostatics is to determine the pressure p(x,y, 2) at the different 
points of the fluid at rest. 


2.2.1 Hydrostatics 


In order to compute the pressure p(z) in a column of water at rest, we consider the balance 
of forces on some cubic liquid element dV = dx - dy - dz (Figure 2.10). Let the pressure 
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on the lower side of the liquid element be p, i.e. the pressure force on the surface element 
dx - dy is fk, = p- dx - dy. The pressure changes over the height of the fluid element dz. 
The pressure change can be written as a Taylor series, up to first order. Thus the pressure 
on the upper side of the fluid element is (p + (dp/dz) - dz +...) and the pressure force 
(p+ (dp/dz)-dz)-dx-dy. The pressure forces on the side surfaces of the fluid element cancel 
each other out, as they are the same all around at any given level and all act normal to the 
surface element. In addition, we have the gravitation G = dm-g = p-dV -g = p-g-dx-dy-dz 
which acts on the center of mass of the fluid element. 


The balance of forces on the fluid element at rest is therefore 
dp 
p- da dy — (p+ ——-dz)dx dy — p-g- dx -dy-dz=0 
z 


If we divide this equation by the fluid element dV = dx-dy-dz we obtain the fundamental 
hydrostatic equation for the pressure change in a column of water due to gravitation 


dp 
—=-—/p- . 2.18 
Teg (2.18) 
This is a first order ordinary differential equation which, on integrating once, yields the 
linear pressure distribution 


PO) = ng: 2+e 


The constant of integration c can be determined with the boundary condition of the given 
problem. For the vessel shown in Figure 2.11 with boundary condition p(z = 0) = po, € = po 
and the linear pressure distribution is 


P(2)=Po-—prg-2 . (2.19) 


Two important consequences follow from this relation. Figure 2.12 shown three containers 
with the same base area and height. The pressure on the base of each container is the same: 
D0 = Pot p-g:h. With the same base area, the pressure force is therefore the same, although 


(p+ dz) dxdy 


Fig. 2.10 : Balance of forces on a fluid el- 
ement at rest 
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Fig. 2.11 : Linear pressure distribution in 
Po Po a gravitational field 
the weight of the fluid in each of the containers is different. This fact is obvious when one 
considers (2.19) and is called the hydrostatic paradox. 


In communicating pipes, the pressure in both sides of the pipe is equal to the outer pressure 
Doo. Therefore the surface of the liquid must be at the same height in both, since, according 
to (2.19), the pressure depends uniquely on the height. 


A U-pipe can also be used as a pressure gauge. If a pressure vessel filled with gas is attached 
to one end of the U-pipe in the right-hand picture in Figure 2.12, an overpressure pı occurs, 
so that there is a height difference Ah in the two sides of the U-pipe. Since the density of 
the gas pc is considerably smaller than the density of the liquid pr, (2.19) can be written 
as 


D1 =Do tpr:g: Ah . (2.20) 
By measuring Ah, the overpressure pı in the gas container can be computed using (2.20). 


A further conclusion can be drawn from the hydrostatic fundamental equation. This is 
known as Archimedes’ principle. For a body of volume Vx completely submerged 
in liquid, the upthrust Fa is equal to the weight of the fluid displaced G. In 
order to derive this law, we consider a cubic fluid element of base area dA and height Ah, 
completely submerged in a liquid of density pr, as shown in Figure 2.13. 


Fig. 2.12 : Hydrostatic paradox, commu- 
nicating pipes, U-pipe pressure gauge 
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The pressure p> on the lower side of the body is larger than the pressure pı on the upper 
side of the body, due to the hydrostatic pressure distribution. The difference in the pressure 
forces Fy and F, results in an upthrust F a, directed vertically upwards. The magnitude of 
this upthrust is 


d|F4| = [Fa] - |Fıl = pp dA pi -dA= (mp — pi) dA 


Using the solution of the hydrostatic fundamental equation (2.19), po = pi + pr 9: Ah, it 
follows that 


dPal = pr: g: Ah: dA = prg: dVk = |a] = | pr:g:dVk=pr:9:Vk , 
VK 


upthrust |Fa| = pr -9 Vg| . (2.21) 


In a rotating vessel, moving with constant angular velocity œ (Figure 2.14), as well as the 
gravitation G, there is also the centrifugal force Z. This causes a liquid surface which would 
be horizontal were there no rotation to become deformed to a parabolic surface. The water 
surface is always perpendicular to the resulting force acting. For an observer rotating with 
the vessel, the liquid is at rest, so that the hydrostatic fundamental equation (2.18) has to 
be extended by the radially acting centrifugal force Z. The magnitude of the centrifugal 
force for a fluid element dV is 

Z| =p- dV -wr 


with r? = z? +2°. Integrating the hydrostatic fundamental equation for a uniformly rotating 


liquid yields 


1 
papeka PW a pre , (2.22) 


In the vertical direction, the pressure decreases linearly with the height z, as in a non- 
rotating liquid. In the horizontal direction it increases quadratically with the distance from 


Fig. 2.13 : Determining the upthrust 
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Fig. 2.14 : Pressure in a rotating liquid 
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the axis of rotation. On the surface of the liquid, with p = po, we find 


1 2 2 
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Therefore the surface of the liquid forms a paraboloid of rotation: 





2.2.2 Aerostatics 


The most important example of aerostatics is the computation of the pressure, density and 
temperature distributions in the atmosphere. Figure 2.15 shows that the pressure p in the 
earth’s atmosphere decreases continually with increasing height. In the different layers of the 
atmosphere, the temperature first decreases to about —56° C’, and then, because of chemical 
processes in the air, increases again. At higher altitudes, the temperature again drops with 
height, only to increase again at even greater altitudes with radiation from the sun. 


The lower layer of the atmosphere is called the troposphere, and can be between 9 and 
11 km thick, depending on the time of year. According to the ideal gas equation for air 
(2.8), the temperature and pressure decrease is associated with a decrease in density, and 
so the cold, heavy air is situated above the warm, lighter air. This air layering is thermally 
unstable, and causes the weather in the troposphere. 


The next layer is the stratosphere, at a height of between 11 and 47 km. Here the temperature 
initially remains constant and then increases again with increasing height. The ozone layer 
forms in this atmospheric layer. This absorbs the UV-radiation from the sun and thus leads 
to a temperature increase. The stratosphere is thermally stable, since the warmer lighter air 
is now layered above the cold, heavy air. This is reason why civil aircraft fly in the stable 
lower stratosphere whatever the weather. 


At altitudes between 47 and 86 km, the air chemistry dominant in the mesosphere again 
causes a thermally unstable temperature decrease. A small decrease in the density of the 
air with increasing height is connected with this chemistry. In this somewhat greater air 
density, small meteors glow and are seen as shooting stars. The dust from volcano eruptions 


40 2. Fundamentals of Fluid Mechanics 


can also reach these heights, and can influence the chemistry of the air considerably for 
many decades. 


The transition to the ionosphere occurs at heights above 87 km. Here the highly energetic 
sun and particle radiation leads to the ionization of the air molecules. This again leads to a 
thermally stable temperature increase with increasing height. ‘The stable temperature and 
density layering in the ionosphere ensures that the atmosphere of the earth does not escape 
from the earth. The electrically charged particles lead to electrical flows in the ionosphere 
which can affect short-wave radio signals considerably. 


In this chapter we restrict ourselves to the lower atmosphere (troposphere and lower strato- 
sphere), where 99 % of the mass of the atmosphere is, and whose temperature, density and 
pressure distributions can be found in the internet as the so-called US standard atmo- 
spheres. 


http: //aero.stanford.edu/StdAtm.html 


The standardized temperature distribution is shown in Figure 2.16. In the troposphere the 
temperature decreases linearly 


T(z) =To +a: (z-z) , (2.23) 
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Fig. 2.15 : Temperature and pressure distribution in the atmosphere 
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with T) = 288.16 K at the surface of the earth z) = 0, the constant a = —6.5- 107? K/m 
and the air pressure at the surface of the earth po = 1.013 - 10° N/m?. 


In the lower stratosphere, the temperature is assumed to be constant 


T = T, = 216.66 K = const. , with plz= 23 =11 km)=p Žž. (2.24) 
For the given temperature distribution, the fundamental aerostatic equation yields the 
associated pressure and density distributions p(z), p(z). Again the starting point is the 
fundamental hydrostatic equation (2.18) 


dp 
og 


However now the density p is a function of the height coordinate z. With the thermal 
equation of state for an ideal gas (2.8) 


p(2) 
=R- p T> = — 
we find the fundamental aerostatic equation 
dp g dz 
ae 2.25 
7 RT (2.25) 


This is again a first order differential equation which is uniquely solvable with one boundary 
condition and a given temperature distribution. 


For the temperature distribution in the troposphere (2.23), we obtain 


1 
dz=-:dT 
a 


Substituting dz into (2.25) and integrating the aerostatic fundamental equation yields 


dp g dT f” g [4 g T 
P I Za] ZI] Z> hP =- MmT 
p RaT Jyp Ra > u) ‚nn > 
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Fig. 2.16 : Standardized temperature dis- 
tribution in the troposphere and lower 
160 220 240 280 
TK stratosphere (US standard atmosphere) 
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, (2.26) 
p(z) follows from p(z), using the equation of state for ideal gases 
pa 2O me) _ OTO _ Tr (a) 
p(z) T(z) po To Po po Lo po T(z) To ' 
e- E 
po To To To í 
(2.27) 





Pressure and density decrease with increasing height in the atmosphere for the given linear 
temperature distribution, according to the power laws (2.26) and (2.27). 


In the lower stratosphere the temperature distribution is isothermal T = T, (2.24) and there 
is an exponential decrease in the pressure. 
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Fig. 2.17 : Temperature and pressure distributions in the US standard atmosphere 
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(2.28) 





The results of the pressure distributions are shown in Figure 2.17. p(z) follows from p(z) 
with the ideal equation of state p(z) = p(z)/(R-T}), 





2.3 Hydro-Aerodynamics, Stream Filament Theory 
2.3.1 Basic Kinematic Concepts 


Before we turn to the computation of flows within the framework of the simplified one- 
dimensional stream filament theory, we first set out the basic kinematic concepts 
for the mathematical description of flows. 


The kinematics of a flow describe the motion of the fluid without taking the forces which 
cause this motion into account. The aim of kinematics is to compute the position vector 
x(t) of a fluid element and thus the dependence of its motion on the time t with respect to 
the chosen coordinate system X = (æ, y, z) for a given velocity field V(u, v, w). 


As in Figure 2.18, we follow the path of a fluid element or of a particle in the flow in time. 
The initial position of the particle motion at time t = 0 is determined by the position vector 
Xo = (20, Yo, 20). At the time tı > 0, the particle has moved along the path sketched to the 
position xX(t,) and at the time tz > tı to the position x(t) etc. The instantaneous position 


particle path 





Fig. 2.18 : Particle path 
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x of the particle under consideration is therefore a function of the initial position X, and 
the time t. The particle path is therefore 


> 


X = f(Xo, t) 
The ordinary differential equation to compute the particle path for a given velocity field 
V(u,v, w) is therefore 

dX 

dt 
This is none other than the well-known defining equation of velocity. The differential equa- 
tions for the individual velocity components read 


-gR t). (2.30) 


dx d dz 

a ulz, y, z, t), 7 =u(x,y,2,t) , u wz, yY, z, t) . (2.31) 
This is a system of first order ordinary differential equations. The particle path is computed 
by integrating these differential equations with the initial conditions x» = x(t = 0). 


For a steady flow, the system of differential equations is not dependent on the time t and 
is written as 
dx 
— = V(X) . 2.32 
7 =) (2.32) 


Here it is to be noted that although 9/öt = 0, the total differential d/dt 4 0. 


Another way of describing flows is by means of streamlines. These indicate the direction 
field of the velocity vector V at a certain time tn (Figure 2.19). Since the tangents are parallel 
to the velocity vector at every position and every time, the determining equation for the 
streamline reads 


vxdx=0 . (2.33) 
The velocity components are then 
u dz v- dz — w -dy 0 v- dz =w dy 
v | x | dy | = |w- dr—u-dzļ=ļ|0 => w-dr=u-dz 


w dz u-dy—v-dı 0 u- dy =v -dz 


streamline 


<4 


x Fig. 2.19 : Streamline 
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This yields the system of first order differential equations for the streamline: 


dz _ w(ia,y,2,t) dz w(x,y,2,t) dy — v(x, y,2,t) (2.34) 





dy  v(z,y.z,t) de ulx,y,z,0) ° de ulx,y,z,0 


The streamlines are again determined by integration, using the method of separation of 
variables, and so are integral curves of the direction field of the given velocity vector v. 


In experiment, or even in a computed flow field, the path lines can be made visible by 
coloring a particle or a fluid element. By photographing the flow region using a long exposure 
times the particle path is made visible. Similarly a picture of the streamlines is obtained by 
marking many particles and photographing the flow field with a short exposure time. The 
picture then shows many short dashes, whose directions indicate the tangent field of the 
velocity vector at the time of the shot. The lines connecting the individual dashes are the 
streamlines. 


The third important way of describing flows is using streaklines. These are shown in Figure 
2.20 and, for a given time t,„, are the lines connecting the positions which the path lines of 
all particles have reached which passed the fixed position Xo at any time tọ < tn. If color or 
smoke is added to the flow field at the position Xo, the snapshots of the colored filaments 
or threads of smoke are streaklines. 


The equation of the streakline at the time t» is 


X = X(Xo, to, t) 5 (2.35) 


where fy is the curve parameter and Xp the family parameter. A parameterless representation 
of the streakline is found by eliminating the curve parameter to. 
For example, the following equation is known from a computation of a particle path: 


. ct) _ 4 _ 
z= RElo) = (A e t G) 


(yo — to +1) e70) +t- 1 y 


streakline 





X Fig. 2.20 : Streaklines 
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Let us consider the equation of that streakline in the (x, y) plane which passes through the 
point (£o, Yo) = (—1,—1) at the time t = 0. If we insert the above ansatz into equation 
(2.35), we obtain 


, y=-be-1 > r4tl1l=to-e , y+tl=-to:e® , 


z+1):-y+D)=-% > b=vV-(cH+l)- y+1) 


to 


to inserted into x = to -e~ — 1 yields an implicit equation for the streakline in the (x, y) 


plane 
r= /—(x +1). (y + 1) -exp (-Va=1)-@+1)) _1 


In steady flows, the particle paths, streamlines and streaklines are the same. In 
unsteady flows the curves are different. 


Let us return to the examples of flows discussed in the introduction in Chapter 1.2. Both 
the flow past a vertical plate and the flow past a wing and an automobile were presented as 
steady flow problems. Now we can interpret the flow lines in Figures 1.13 and 1.18 as particle 
paths and streamlines respectively. Small particles of aluminum are added to each flow in 
the water channel and a snapshot of these with a long enough exposure time characterizes 
the structure of the steady flow. In Figure 1.21, the wake flow of an automobile in a wind 
tunnel was made visible with smoke allowed into the free stream at a fixed position Xp. All 
smoke particles have passed through the same position, so that it is the streaklines which 
we see in this shot. Figure 2.21 supplements the sketches of particle paths, streamlines and 
streaklines, which all lie on the same curve in these three examples of steady flows. 


For unsteady flows, the particle paths are different from the streamlines and streaklines, 
making the interpretation of unsteady flows difficult. We will illustrate this using a simple 
flow example. In Figure 2.22 we move a sphere with constant velocity us through a fluid 
at rest. As the sphere is moved along, the particle path forms a loop, while the snapshot of 
the streamlines indicates closed curves. This is the flow field which we see as an external 
observer at rest. The streamline picture looks completely different if we move with the 
sphere. We then see the constant free stream us coming towards us and the flow becomes 
time independent. Instead of the closed streamlines, steady streamlines form from left to 
right, along the same curves as the pathlines and streaklines. Depending on our frame of 
reference, the flow field looks completely different. Expressed physically, this means that 
the streamlines and particles paths are not invariant when the inertial frame is changed 
(position transformation with constant translation velocity). 


Two further examples of shear flows serve to extend these ideas. Let us consider a plane 
wave in a plate boundary layer flow. The u component of the velocity displacement is 
written as 


u(x, 2,t) = i(z)- eh (ar—wrt) 


with the amplitude function &(z), which is a function only of the vertical coordinate z, the 
wave number a and the angular frequency w. The phase velocity c of the wave is c = w/a. 
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The observer at rest sees circles for the particle paths and streamlines of the wave moving 
past with the phase velocity c, as sketched in the snapshot in Figure 2.23. The observer 
moving with the wave sees the plate moving with the phase velocity c and a streamline 
picture which looks like cat’s eyes. 


£ N flow past a plate 


particle path 
streamline 
streakline 


flow past a wing 


flow past a car 





Fig. 2.21 : Particle paths, streamlines, streaklines of the steady flow past a horizontal plate, 
a wing and an automobile 
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The third example of shear layer flow is the wake flow of a cylinder, already met in Chapter 
1.1 as the Karman vortex street in connection with the collapse of the Tacoma bridge. The 
humming of high voltage wires in the wind is also due to the periodic flow separation of the 
Karman vortex street at the cylindrical cross-section. Figure 2.24 shows first the streaklines, 
particle paths and streamlines of a cylinder moving through a fluid at rest with constant 
velocity Us. for an observer at rest. The observer moving with the vortices periodically 
swimming downstream with phase velocity c again sees the perturbed streamlines as cat’s 


particle path 
observer at rest 


streamlines 
observer at rest 


streamlines 
moving observer 





Fig. 2.22 : Flow past a sphere, with ob- 
server at rest and moving 
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eyes. The historical photographs of PRANDTL 1929 and TIMME 1957 in the water channel 
and the theoretically computed streamlines of VON KARMAN 1912 are shown in Figure 
2.25. 


As we have already learned from the examples of unsteady flows, the description of the 
kinematics of unsteady flows in particular is very difficult. Much practice and experience 
is required to be able to interpret the experimental results in the wind tunnel or the flow 
simulations on a computer correctly. However it is precisely the kinematic description of the 
flow which yields important insights into the structure of a flow. 


Now that we have seen that the flow picture is dependent on the frame of reference, there 
are basically two different ways of treating a flow mathematically. In the Euler picture 
we assume a fixed observer. This manner of description corresponds to using measuring 
apparatus which is fixed in position to measure local flow quantities. 


The Lagrange picture assumes a frame of reference moving with a particle or fluid element. 
The mathematical relationship between the two picture is, for example for the acceleration 
of the flow b = dv/dt = d’X/dt*, the total differential of the velocity vector V(u, v, w). 
For the u component u(x, y, z, t) of the velocity vector we have 

_ Ou Ou Ou Ou 


ne aie ee 
Mahn 


So the total time derivative of u is 
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Fig. 2.23 : Wave in a boundary layer, observer at rest and moving 
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with 


and we have 


(2.36) 





In this equation 


S substantial rate of change, Lagrange picture, 

L local rate of change at a fixed position, Euler picture, 

C convective spatial changes due to convection from place to place, effect of the velocity 
field V = (u, v, w). 
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Fig. 2.24 : Kármán vortex street, observer at rest and moving 
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For the acceleration b of the flow field which we will need in the following chapters, we 
obtain 

dv OW OV OV OV OV 

ae Se ae, — prne v:-V)v 2.37 
io aa ee a0) 
with the nabla operator V = (0/0x,0/0y,0/0z) and (v-V) the scalar product of the 
velocity vector Vv and the nabla operator V. 


b = 


For Cartesian coordinates this yields 
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Fig. 2.25 : Streaklines, particle paths and streamlines of the Kármán vortex street 
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and for (V-V)-v 
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In the case of a steady flow, all partial derivatives with respect to time vanish, so 0/0t = 0, 
while the substantial derivative with respect to time d/dt can indeed be non-zero when 
convective changes occur. In unsteady flows both 0/0t 4 0 and d/dt F 0. 


2.3.2 Incompressible Flows 


The fundamental equations and methods of one-dimensional stream filament theory 
are still used today in industry for the preliminary design of new products. It is for this 
reason that one-dimensional stream filament theory is a good introduction to the theoretical 
treatment of flows. The solution software of the algebraic system of equations to be treated 
will be presented in Chapter 3.1. 


We will denote the one-dimensional velocity component with c(s). This is a function of 
s only, called the stream filament coordinate. In order to introduce this one-dimensional 
coordinate s it is first of all useful to introduce the concept of the streamtube. If the 


streamlines form a closed surface, this surrounding surface is called a streamtube (Figure 
2.26). 


stream tube 


stream filament 





Fig. 2.26 : Streamtube and stream fila- 
ment 
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Since the streamlines are by definition the tangents to the velocity vectors, no fluid mass will 
pass through the walls of the streamtube. This means that channels with fixed walls form 
streamtubes. If the changes in the flow quantities over the cross-section of the streamtube are 
small compared to the change along the streamtube, the (approximately) one-dimensional 
changes in the flow quantities along the hypothetical stream filament can be computed. 
The coordinate along the stream filament is called the stream filament coordinate s. Along 
a stream filament, for incompressible and steady flows, we have 


c=c(s) , p=pls) , A=Als) 


All flow quantities as well as the cross-section A of the streamtube are functions of the stream 
filament coordinate s only. For the problem of the flow past an automobile, for example, 
stream surfaces can be determined, corresponding to the streamtubes of channel flows. 
Figure 2.27 shows such a stream surface past an automobile. If the changes perpendicular 
to the stream surface are small compared to the changes along the streamline, as is the case 
along the mid-section of the vehicle flow, a stream filament can again be determined, along 
which the variation in the flow quantities is approximately one-dimensional. 


The fundamental equations for one-dimensional stream filament theory are as 
follows: 


Conservation of Mass: 
The mass flux entering a streamtube my, is equal to the mass flux exiting the streamtube 
my. Using the volume flux V; and V2, we find 


ry = 01-Vi = pi-C1 A = po: oa A=-m him , 


| (2.38) 





Fig. 2.27 : Stream surface and stream fil- 
ament 


54 2. Fundamentals of Fluid Mechanics 


Conservation of Momentum or Equation of Motion: 

We first formulate the equation of motion for a stream filament placed in the inviscid outer 
flow or the inviscid core flow of a channel. From the balance of forces along a selected stream 
filament element dV (Figure 2.28), the variation in cross-section along the stream filament 
can be neglected to first approximation. The equation of motion reads mass - acceleration 
= sum of all external forces acting. Therefore for the volume element dV we have 


dm-b=S oF, . (2.39) 


We have already considered the acceleration b in Chapter 2.3.1. For a one-dimensional stream 
filament, equation (2.37) is written as 


Ee LC race 
=a U mpg” 


for this flow which is assumed to be steady c- (dc/ds). The mass of the volume element 
dV considered in Figure 2.28 is dm = p-dAds. The forces acting on the volume element 
are the pressure forces and the gravitation, and their components along the stream filament 
coordinate are in equilibrium. This yields 


dc Oc Oc 
Asse eher = = 
p-dA-ds ; p-dA-ds ( er 3 


ð 
p:dA— (p+ 2 ds) -dA — p- g- dA- ds- costo) 


cos(p) = dz/ds and division by p-dA- ds delivers the Euler equation for the stream 
filament 


(2.40) 





For steady flows, all quantities are functions only of s and it follows that 


de d/e 1 dp dz © 1 


dz 





—— = COS 
3 ds x Fig. 2.28 : Balance of forces at a stream 
p dA p g dA ds filament element dV 
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Integration along the stream filament s from position 1 with c1, pı and s1, 2, to position 2 
with c 9, po and se, z2 yields 


N| = 


9 9 P2 1 
e-+) p Ptg (=)= 


Pl 


For the incompressible flow under consideration p = const., so that the factor 1/p can be 
written in front of the integral and we obtain the Bernoulli equation for incompressible, 
steady, inviscid flows. The dimensions are energy per unit mass: 


(2.41) 





Alternatively, the Bernoulli equation with the dimensions of energy per unit volume is 
frequently used 


1 1 
m+z P FP g a= p+ P tpg: a = const. . (2.42) 





At some given position, the Bernoulli equation for steady flows reads 


2 

P+ 3:0: 40:9: 2= const. or zt g+= const (2.43) 
Here the constant summarizes the three known terms at an initial state. It has the same 
value for all points along s of a stream filament, but can vary from one stream filament to 
the next. The Bernoulli equation is an algebraic equation and yields the relation between 
velocity and pressure. For unsteady flows, the partial time derivative Oc/Ot of the Euler 
equation must also be integrated along the stream filament s. Here the integration is to 
be carried out at a fixed time t from sı to sa. The Bernoulli equation for unsteady one- 
dimensional flows is found: 


2 oc 1 9 
p: — -ds + p2 +z: p +p: g: 22 = const. , 
S1 


ot 2 
F Oc ds +? + 2 +g: 22 = const (2.44) 
a Ot p 2 Joo | 


Application of the Bernoulli Equation 


There are many examples of application of the Bernoulli equation in the book of solved 
problems accompanying this textbook. We will look at two examples which are applied in 
practice. Using the Venturi pipe shown in Figure 2.29, the mass flux can be determined 
by measuring the pressure at the narrowest cross-section and applying the Bernoulli equa- 
tion (2.41). The narrowing of the cross-section causes an acceleration in the nozzle, and, 
according to the Bernoulli equation, a pressure drop. A widening in the cross-section causes 
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a deceleration of the flow with corresponding pressure recovery. If the pressure p is mea- 
sured at the narrowest cross-section A, and for known values of cı and p4, the velocity c is 
determined using 


Cp Pt 


Since in this example 2; = 2, the gravitation term is eliminated. For a known cross-sectional 
area at the narrowest cross-section A, the mass flux is determined to be 


m=p-c:A 


The application of the Bernoulli equation (2.41) therefore makes it possible to determine the 
flow velocity c from a measured pressure p. This is used in the case of aircraft, for example, 
to determine the flight velocity using a Prandtl tube. Before we look at the way a Prandtl 
tube works, we first have to introduce some pressure definitions. We consider the Bernoulli 
equation (2.42) 


P+ 2:0: 49-9: 2= const. 
Here we denote p = Pstat as the static pressure and (1/2) -p-c? = payn as the dynamic 
pressure. The static pressure Pstqz 18 the pressure measured when one moves with the flow 
velocity c in the fluid. Therefore the static pressure is that which is responsible for the 
pressure force which acts on a body in a flow. The dynamic pressure payn can be considered 
to be a measure of the kinetic energy per unit volume of a volume element of the fluid 
moving with velocity c. 


P, P. stream filament 
Te 


C] ĉi 
nozzle diffusor 
acceleration deceleration 


pressure drop 


pressure increase 





Fig. 2.29 : Venturi pipe 
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In the case of shear flows, such as the boundary layer flow past a wing, 21 = 2% and the 
gravitation term p- g- z is eliminated from the equation. The constant on the right-hand 
side of the Bernoulli equation can vary from one streamline to the next. It is a property 
of the streamline under consideration and is determined by suitable reference values. Such 
reference values can be, for example, the known values of the undisturbed free stream, such 
as Pæ and cæ. In the case of the flow past a wing, the constant on the so-called stagnation 
streamline can be determined. This stagnation streamline leads from the free stream in 
infinity, past a variable point 1 to the stagnation point 0 on the wing (Figure 2.30). 


The Bernoulli equation on the stagnation streamline reads 


1 2 1 2 
Boot 5° P* Coo = Pit 5° 0° C1 = Po = const. 


In the stagnation point c = 0, and so there is no dynamic pressure at this point. The variable 
po denotes the pressure in the stagnation point. This is also called the reservoir pressure or 
the total pressure. Therefore we have 


Po = Ptot = Pres = Pstat + Pdyn . (2.45) 


We have already used the dynamic pressure of the free stream (1/2)-p-c%, in the introductory 
chapters for the dimensionless pressure coefficient cp 
P— Poo 


1 2 
a °C 


Cy = 


The different pressure definitions are summarized in Figure 2.31. The pressures can be 
measured using classical hydrostatic methods. 


Measurement of the static pressure Pstat: 


The most simple measurement principle to determine the static pressure Psa, Consists of 
bore hole in a wall and the U-pipe pressure gauge introduced in Chapter 2.2.1. The static 
pressure Pstat Of the outer flow is imprinted onto the boundary layer, i.e. it is constant within 
the boundary layer in the direction normal to the wall. The static pressure of the outer flow 
is measured as follows. The relations in Figure 2.32 between pressure difference Ap and 
height difference Ah in the pressure gauge are valid, with pa the density of the air, p, the 
density of the liquid, and prer the reference pressure. 


Figure 2.33 shows the wind tunnel model of a wing. The pressure bore holes, from which the 
pressure distribution in Figure 1.18 was taken, are so small that they cannot be seen in the 


Pytat 





Fig. 2.30 : Pressure definitions in the flow 
stagnation point past a wing 
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static dynamic total Fig. 2.31 : The different pressure defini- 
Pisces Dr alle Pee tions: static pressure Pstar, dynamic pres- 
tat Payn Piot SUTE Payn, total pressure Prot 


figure. It is only the small pressure tubes inside the wing model which can be seen. These 
measure the pressure, and these days are piezoquarz pressure sensors rather than classical 
U-pipe pressure gauges. 


The static pressure Psat can be measured with a probe held in the flow (Figure 2.34). It 
works according to the same principle as the wall bore holes, but in this case the bore holes 
are distributed around the probe to measure the static pressure. Here too a boundary layer 
forms around the top of the probe and this is imprinted with the static pressure of the 
outer flow. In order to minimize measurement errors, the bore holes must be at a large 
enough distance from the top of the probe and the shaft of the probe. This is so that the 
disturbances caused by the holes die away and are not measured too. 





Deiat Pag: h= De} + pre gs Ah 
Paa = Mre =A g Ah pa gh 
frequently: pg: g-h<« py-g- Ah 
=> AP = Paat — Prep = pt. g: AN 


Fig. 2.32 : Measurement of the static pres- 
SUTE Pstat 
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Measurement of the total pressure Prot or reservoir pressure po: 


The measurement of the total pressure Piot or the reservoir pressure po is carried out using 
a Pitot tube. By placing this into the parallel flow, the tube will be filled with air for a few 
moments, until the air everywhere in the tube has come to rest, since it cannot expand any 
further. This is also true for the entrance cross-section, where the stagnation point forms 
with c = 0. It then follows that the same total pressure p;o is at hand everywhere in the 
Pitot tube, and this again is measured by a U-pipe pressure gauge. 


Measurement of the dynamic pressure Payn: 


In order to measure the dynamic pressure Payn, a combination of a static probe and a Pitot 
tube is used, called the Prandtl stagnation tube. This determines the dynamic pressure 
as the difference between the total pressure and the static pressure, and from it the velocity 


can be determined. 
Dias 2. ge Ah 
eng Sal), (2.46) 
PL PL 


An example of a Prandtl tube, as can be seen on any aircraft, is shown in Figure 2.35. 


Fig. 2.33 : Static pressure bore holes in a 
wing model 
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static pressure Pstat 
Pageh Kpg Ah 
Ap = Pstat — Pref 5 pt: 9: Ah 


total pressure po: 
Ptot = Pref + Pl -g- Ah 


dynamic pressure Payn 
Ptot = Pstat + pr g- Ah 

Pdyn = Ptot — Pstat 

Pay = £- pa E = p- g- Ah 


c — , /ZPim — , [2prg:Ah 
Pa Pa 


Fig. 2.34 : Measurement of the static pressure Pstat, the total pressure Piot and the dynamic 


pressure Pdyn 
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Fig. 2.35 : Prandtl stagnation tube 


Balance of Forces Perpendicular to the Stream Filament 


Up until now we have discussed examples of flows where the changes along the stream 
filament have been large compared to the changes perpendicular to the stream filament by 
definition. If the changes are larger than those along the stream filament, we carry out the 
balance of forces on the volume element dV perpendicular to the stream filament along the 
normal direction n. s now denotes the arc length of the stream filament, r is the local radius 
of curvature. The equation of motion normal to the stream filament reads 


dm- bn = XR, 


The mass element dm is found to be dm = p-dV = p-dA-dn. If c denotes the velocity along 
the stream filament coordinate s, the magnitude of the acceleration b, is determined from 
the ratio of the magnitude of the centripetal force F, and the mass element dm. Therefore 
we have 

o dm- e z oe 


= => n= = 
r dm r 








This acceleration b,, keeps the mass element on the curved path, and so its direction indicates 
the local center of curvature, opposite to the direction of n. The external forces are the 
pressure forces and one component of gravity p-dA-dn-g (Figure 2.36). Therefore the 
equation of motion becomes 


2 


ð 
dm: b)„ = p-dA-dn- (-<) =p-dA— (p+ an) -dA+p-dA-dn-gsin(y) , 
r n 


and after division by (—p-dA-dn) and with sin(y) = —dz/dn it follows that 


@ 1% dz 
ea i 2.47 
r pon ae dn ( ) 
For a plane shear flow at z = const., we find using dz = 0 that 
a- rg 
eng (2.48) 
r pon 


The pressure increases in the direction of the external normal direction n, or in the case of 
plane circular flows, in the radial direction r. The pressure force and centripetal force are 
in equilibrium. 


Vortex motions on concentric circular paths can be computed with the ordinary differential 
equation (2.48). For example, the pressure and velocity distributions of a tornado (Figure 
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dn 


dz Fig. 2.36 : Balance of forces at a volume el- 


PN ement perpendicular to the stream filament 





2.37) can be determined approximately using one-dimensional stream filament theory. The 
streamlines are concentric circles. For a velocity magnitude of c, on circular paths we have 
c(r) = ¢,/r with the circumferential velocity c(Ro) = co at a given radius Rp and the 
constant ¢, = co: Ro . With n =r, equation (2.48) is written 


2 1 dp 
ŽŽ =. , 2.49 
r? pdr ( ) 


Integration of this first order ordinary differential equation with the boundary condition 


p( Ro) = po at a fixed radius Ro yields 
2 
Pc 1 1 
A) 2.50 


This can be written in the following form 


-Ê (r) = po + 5 -c = const.| . (2.51) 








Fig. 2.37 : Circular flow in a tornado 
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This is the Bernoulli equation for vortex flows on concentric circles, derived from the balance 
of forces perpendicular to the stream filament. It can be shown that the flow on concentric 
circles is irrotational, with V x ¢ = 0. Figure 2.38, for r > Ro shows the pressure distribution 
computed with equation (2.50) as well as the assumed velocity distribution ce = ¢,/r. 


According to the Bernoulli equation (2.51), pressure and velocity behave differently with 
increasing r. For r > 0, the velocity for a potential vortex would increase indefinitely. Since 
this is not what happens in reality, the differential equation for inviscid vortex flows (2.49) 
is replaced by the differential equation for viscous flows for r > Ro. This equation will 
be treated at the end of this chapter. Here too the division of flow regimes introduced in 
Chapter 1.2 is seen. In the core of the vortex the flow corresponds to a frictional rigid body 
rotation with constant angular velocity w, and linear velocity distribution c = wp- r. The 
pressure drops for r < Ro again, and in the case of a tornado, reaches values of between 20 
and 200 mbar. 


Conservation of Energy 


The third fundamental equation needed for the complete mathematical description of flows 
with heat transport, which is also needed if the work done by flow machinery is to be taken 


p 
p 
0 
/ 
7 
7 
R, ‚ 
C rigid body 
rotation 
potential vortex 
C 





r Fig. 2.38 : Pressure and velocity distribu- 
tion in a potential vortex 
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into account, is the conservation of energy. In order to derive the balance of energy, we add 
to the sketches of the streamtube and the stream filament in Figure 2.26 an additional heat 
flux q (Figure 2.39). 


In general the conservation of energy of a steady inviscid flow states that the 
change in energy flux in a given volume element dV is equal to the work done 
by the forces acting and by the heat flux. The energy flux FE written in Watts ({W} = 


{J/s}) is 
E= (cS) (cS) pc 


with the internal energy e referred to the mass element dm = p -dV and the mass specific 
kinetic energy c*/2. For both cross-sections A; and A» of the streamtube considered, and 
with the continuity equation m = const., we have 


oy a a 
F = GIT ‘m= A Girt Ad y 


(eV. & 
N fat > -m = fat > - Px + C+ Ag 


The work done by the forces acting (pressure forces and gravity), as well as the power of 
the heat flux q-m, lead (neglecting the friction) to a change in the energy flux from 1 to 2 
according to the following balance equations 


Eo — By =p, +A, +61 — po: Ago +9: (4 —m%)-mt+aq-m 


E l e 
(a+ 2) “Mm — (a+) -m =p Ai e pr Az-e +g (a za) tam 
After division by m= p,-¢,- Ay = fo- Co- Ag it follows that 


Da Pı 
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Fig. 2.39 : Streamtube and stream fila- 
ment with heat flux q 





2.3. Hydro-Aerodynamics, Stream Filament Theory 65 


With the definition of the mass specific enthalpy h = e + p/p we find 


1 1 
ata Gtgmamt+s-atg-arta| 
The three quantities h,,c, and g- 2, at the cross-section A, which are known quantities 
according to hy + (1/2)-c + 9-2 = const., can be brought together to one constant. 
Considering the quantities at the cross-section A, we obtain 


1 
h+ 5: Č tgz q= const. i (2.52) 


If no heat is supplied or removed and so the internal energy does not change, the balance 
of energy and the Bernoulli equation are identical. This is only valid for the incompressible 
flow considered in this chapter. 


For flows with mechanical supply of energy (pumps) or mechanical removal of energy 
(turbines), such as are found in the night storage power station in Figure 2.40, the balance 
of energy (2.52), or, if the heat losses in the pumps and turbines are neglected, the Bernoulli 
equation (2.42), is extended by a term which describes the specific work done Alp by the 
pump. The same holds for the turbine with specific work done Aly (units {J/m?}). 


During the day, when the load is at its maximum, the water in the storage power station 
flows from the reservoir at height z2 along the pressure pipe down to the collection basin 
at height z,, and the flow drives the turbine. At night, when the network load is lower, the 
water is pumped back up from 21 to zo by the turbine which now acts as a pump. While 
the water being pumped from positions 1 to 2, energy is being supplied to it. The energy 
content per unit volume {J/m?} of the fluid is therefore larger at 2 than at 1 


1 1 
P+ sp Gtp- g a>pmtz pate gz 


2 




















Fig. 2.40 : Storage power station 
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Using the specific work of the pump Alp > 0, the Bernoulli equation (flow direction from 
1 — 2) reads 


1 1 
pet Pate g =p t5 tpg a+ Alp] . (2.53) 





If the fluid flows from 2 to 1, thereby driving the turbine, energy is removed from the fluid 
on the way from 2 to 1, and so the energy content of the fluid is lower at 1 than at 2 


1 1 
Dt zea tp ga <m+tz pP &+p 9% 


If we define the specific work transformed by a turbine into electrical energy to be positive 
too, i.e. Alp > 0, the Bernoulli equation in this case (How direction from 2 — 1) reads 


1 1 
Pitz PA +P Gg A=M+ Zp Qtp-g:2—-Alr (2.54) 





Note that the flow direction is different when the pump is used compared to when the 
turbine is used. 


Using the specific work Al for the pump or the turbine, we can obtain their power in 
{W} = {J/s} by multiplication with the volume flux V = A-c as 


L=AI:V 


Summary of the Fundamental Inviscid Equations of Stream Filament Theory 


The fundamental equations of one-dimensional stream filament theory for incompressible 
and inviscid flows can be summarized as follows: 


Conservation of mass p-c:A=const. (2.55) 


Conservation of momentum s 
. Oc p 1, 

Integral of Euler equation = J —-ds+—-+-=-@+g-2= const. (2.56) 

Bernoulli equation 


Ot p 2 


1 
Conservation of energy h + 5° C+g-2+q+Al=const. (2.57) 





These are 3 algebraic equations to determine the flow variables c, p, h. They are supple- 
mented by the thermodynamic relations 
OT 
h= c- T = —\-— 
p Zn Os 
The work done by flow machinery Al in general has to be measured. The solution of the 
algebraic equations (2.55) to (2.57) is found using either known algebraic methods, or, 
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if possible, analytically. The numerical solution will be presented in Chapter 3.1 in the 
form of the software package KAPPA (Karlsruhe Parallel Program for Aerodynamics) 
Stromfaden (stream filament). Examples of analytic solutions are found in the book of 
solved problems in Chapter 2.3. 


Navier-Stokes Equation 


Now at the end of this chapter on incompressible flows, we will consider two-dimensional 
conservation of momentum, that is the equation of motion, of viscous flows at fixed walls. 


Again we set out the stream surface and the stream filament as in Figure 2.27 in the 
boundary layer flow regime or the viscous wake of the flow past an automobile. Once more 
we select a cylindrical volume element along the stream filament and consider the streamtube 
of Figure 2.41 for the viscous flow. Here a cylindrical ring element of length ds and base 
area dA = 2-7-r-dr is considered. The velocity c is no longer only a function of s and 
possibly of t, but is now also dependent on the radial coordinate r. Since 0c/Or 4 0 holds 
for r Æ 0 there are shear stress contributions to the balance of forces. For the equation of 


motion 
dm -b = 2. F, 


and mass dm = p-dA-ds = p-2-7-r-dr-ds, acceleration b = 0c/O0t+ c¢- (Oc/0s) and forces 
acting Fy: pressure forces, shear stresses and the component of gravity along s, we have 


im. (Bee) pdemeredr-ds: (Tre. =p-2-7-r-dr— 


at’ Os al Os 
(p+ Bas) 2m rear po ge 2n-r dr ds-coslg) -7:2:7 r: ds+ 
(r+ ar) -2-a(r+dr) - ds 
-į (t+ ot dr) 27 (r + dr) ds Jo 





“pg 2nrdr ds 


Fig. 2.41 : Balance of forces at a stream filament element dV for viscous flow 
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With cos(y) = dz/ds and dividing by (p-2-7-r-dr-ds), while neglecting terms of the 
order (dr)?, the ansatz 7 = u: (dc/ör) with v = u/p yields the Navier-Stokes equation 
in cylindrical coordinates 


(2.58) 





This is a second order partial differential equation. In contrast to the Euler equation (2.40), 
the Navier-Stokes equation also takes the effect of friction into account, by means of changes 
in the shear stresses which give rise to the second derivative of the velocities. The left- 
hand side of the Navier-Stokes equation again characterizes the fundamental equation of 
kinematics(2.36) for one-dimensional flows, but is now supplemented by the pressure, friction 
and gravitation forces. 


The Navier-Stokes equation in stream filament coordinates s and n reads 


(2.59) 





The individual terms mean: 


Oc Oc. 
—-+c-— inertial forces per unit mass, 
ot OS 
1 Op . 
=. — pressure force per unit mass, 
p Os 
0c >. 
UV — friction force per unit mass, 
On? 
dz 
g: ds gravity per unit mass. 


The Navier-Stokes equation can be made dimensionless with suitable characteristic quan- 
tities. The dimensionless quantities are denoted by a star. All position coordinates appearing 
s, n and z are referred to a characteristic length L and the velocity c to a characteristic 
velocity Coo. The ratio L/c.. represents a characteristic time, used to make the time t dimen- 
sionless. The pressure p is made dimensionless with twice the value of the dynamic pressure, 
i.e. with p- cĉ. 

x C b+ ex ’ p 


? L ? Coo ? L ? P pc 








If we set these quantities into the dimensional Navier-Stokes equation (2.59), we obtain 


ade. 0 1 p-c Op* Co O’ L dž 


L æ L Os p L ðs TT Dn? 














Multiplication by the factor 4 leads to 


Oc | č Op v Oc gL dz 


ae Os As Co: lb On ce ds* 
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The combinations of characteristic quantities in front of each of the last two terms are the 
Reynolds number Re, = (cs : L)/v and the Froude number Fr; = c2,/(g- 1), each formed 
with the characteristic length L, which we already used in the introductory chapters. The 
dimensionless Navier-Stokes equation reads 


Oc 


dc“ Op* 
eR es + 2. 
ar l Ds ? 2) 





with the dimensionless characteristic numbers 














ä i Oc 2 
inertial force c- S C 
Froude number: Fr, = —————— = as eu 
gravity Ge geL 
inertial force C- 5 Cor 
Reynolds number: Re, = — ——— = Er 
viscous force p- E V 


For Fr > 1, the inertial force of the flow dominates and gravity may be neglected. For 
Re > 1 the inertial force also dominates. The effect of viscosity is restricted to a thin 
viscous layer close to the wall, which we recognize as the boundary layer. For the 
boundary layer thickness ô referred to the length L and shown in Figure 2.42, we have the 
relation 


1 


2 Nm ——— (2.61) 
L Re; 

The static pressure inside this boundary layer corresponds to the static pressure of the 
inviscid outer flow, imprinted onto the boundary layer. For Re < 1, the viscous force 
dominates in the entire flow field. Flows in this regime are called creeping flows (Figure 
2.42), and a division into an inviscid outer flow and a viscous flow layer close to the wall is 


no longer possible. 





Fig. 2.42 : Plate boundary layer flow and 
creeping flow past a cylinder 
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The orders of magnitude of Reynolds numbers which occur for living organisms and in 
technology are summarized in the following table. 


Nature Rer, motion 


bacteria friction dominates 
motion > ciliary 


amoeba (flagella) 


tadpoles inertial force dominates 
jet propulsion 


eel wave-like motion 


large Rer, numbers 
motion by means of vortex separation 


blue whale => tail fin 





Technology 


automobile combustion engine 
airplane 


submarine propeller 





If we integrate the dimensionless Navier-Stokes equation (2.60) at a fixed time t along the 
stream filament coordinate, we obtain 


O , 9 (c* . 
[= -ds fo GE = 


Op* 1 dz* 1 0? c* 
_ . d — —. . d + —. . d3* t. 
J Os* ° Fr; J ds* S oT Rer J On*2 ST CONS. , 


1 Ore 
J -ds* =const.| . (2.62) 




















zZ _— — 
Re; 


In the KAPPA software, equation (2.62) supplements the Bernoulli equation (2.56) by the 
viscous term (1/Re;,) - [(0?c*/On**) - ds*. In the software package KAPPA-Stromfaden 
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(Chapter 3.1), this is taken into account to compute the viscous layers. We are then no longer 
dealing with one-dimensional stream filament theory, rather with the general formulation of 
the fluid mechanical fundamental equations for three-dimensional flows. 


Analytical Solutions of the Navier-Stokes Equation 


We consider three analytical solutions of the Navier-Stokes equation. In a pipe with 
circular cross-section of radius R, a parabolic velocity profile e(r) is formed (Figure 2.43). 
This is a steady (Oc/0t = 0), fully formed (0c/0s = 0) pipe flow. The velocity profile does 
not change along the coordinate s, so that (1/p) -Op/d0s = const. This is a horizontal shear 
flow with dz = 0, and so the gravitational force g - (dz/ds) = 0 is eliminated. Under these 
assumptions, the Navier-Stokes equation in cylindrical coordinates (2.58) reads 


1 de „de 
r do dr? 
where the constant viscosity v has been added on to the constant pressure gradient (1/p) - 
Op/Os. Since the velocity c(r) is a function only of the radial coordinate r, we obtain a 


second order differential equation. With the two boundary conditions 


= const. (2.63) 


? 


reiR. = «che 
and the supplementary condition 


de 


aia — 0 
dr 


r=0 





the differential equation (2.63) can be solved using a power law ansatz for e(r): 
R? dp 1 
Se ee 
ay) 4-v-p ds ( =) 


With the maximum velocity Cmax = —(R?/(4-v-p)) - (dp/ds) it is found for the pipe flow 
that 





2 


c(r) = Cas (1 = =) | (2.64) 


A parabolic velocity profile e(n) is also found for the plane steady channel flow, called 

the Poiseuille flow (Figure 2.44). The Navier-Stokes equation to be solved (2.59), with 

Op/Os = const. and dz = 0 for the fully-formed channel flow 0c/0s = 0 can be written as 
2 


ð 
V: s = const. . (2.65) 


Ak L 


h 
UO A S 
0 = Fig. 2.43 : Hagen-Poiseuille pipe flow 
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After integrating twice and using the boundary conditions 
me =- e= 


we find the parabolic velocity profile 


h? dp n? n? 
en Seller leer. er 2.66 
AD) 2-v-p ds ( =) © ( =) en 


The shear stress of this viscous channel flow can be computed with (2.1) 





de 2- U- Cmaz 


The linear distribution of the magnitudes of the shear stresses shown in Figure 2.44 is 
therefore found. 


In the case of Couette flow, shown in Figure 2.45, in a channel with a lower wall at rest and 
an upper wall moving with constant velocity U, and the additional assumption Op/Os = 0 
the Navier-Stokes equation (2.59) is 


dc 
dn? 
After integrating twice and using the boundary conditions 


n=th => c-k)=0, c-+h)=U 


=0 . (2.67) 


we obtain the linear velocity profile 
u n 
=<-(1+5) . 2.68 
oa)= 5. (145 (2.68) 
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Fig. 2.44 : Poiseuille channel flow 
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2H c (n) 





Fig. 2.45 : Couette flow 


Drag Coefficients 


Now that we have looked at the basics of inviscid and viscous flow regimes, we can go back 
to the introductory examples in Chapter 1.2 and treat the drag of bodies in a flow more 
precisely. 


The total drag coefficient cw (1.2) 


W 


Giesen ar 
1 


? 


with drag W on the body, free stream velocity c,, and a characteristic cross-sectional area 
A, consists of terms corresponding to the inviscid and viscous regions in the flow field: 


, (2.60 


the part due to the pressure distribution c, called pressure drag cg and the friction drag 
cr. The drag coefficients are written as 


P Toy 


Cd , Cf T 
2 
2° Poo * x 


J 


1 


with the shear stress 7, at the wall. The pressure force P is computed from the pressure 
coefficient cp (1.1) 

= P — Po 

p= T 
9° Po* Ce, 


by integration along the wall streamline s, and the magnitude of the total wall shear stress 
Tw by integration of the local value of 7 on the surface of the body. This yields the total 
drag W of a body of length L to be 


L L 
W = ¢ crasina: ds | Cpi SiN Q + ds 
0 0 


i i 
1 
+f cjucosa -ds + | epeosa:ds) +5 ps CB : (2.70) 
0 0 
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2 | Fig. 2.46 : Pressure coefficient c, and drag 
0 coefficient cr of symmetric flow past a foil 
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where u and / are the upper and lower sides of the body respectively, and B is a characteristic 
depth, with A = L - B. The integration is carried out along the two surfaces. In splitting 
the drag up into the pressure drag and the viscous drag, it was assumed that, although the 
pressure drag depends on the shape of the body, the friction drag essentially only depends 
of the size of the surface of the body and not on the shape of this surface. 


Figure 2.46 shows the pressure drag coefficient cq and the friction drag coefficient cp for a 
symmetric foil in a flow with Cə. Here we note that, in contrast to the example in Chapter 
1, we now assume an incompressible flow with low Mach number, such as is found in the 
case of gliders. Figure 2.47 summarizes the drag contributions of bodies in flows. 


Se los 
w4 0% 


Ser Fig. 2.47 : Contributions due to the pres- 
g 3 sure drag cq and the friction drag cy of a 
SE 2 3 8 100%} 0% body in a flow 
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Fig. 2.48 : Friction coefficient cy of the 
1 s/L plate boundary layer 


The pressure is imprinted onto the boundary layer of a plate placed longitudinally in a flow. 
Here the pressure drag cg is equal to zero and the total drag Cw consists only of the friction 
drag cr, whose local value along the plate is shown in Figure 2.48. 


Because of its small cross-sectional area A, a slender foil only has a small pressure drag 
(Figure 2.47), and it is the friction drag which dominates. In the case of a cylinder in a flow, 
the relationship of the different drag terms is the opposite, and it is the pressure drag which 
dominates. A plate placed transversely in a flow has essentially only pressure drag, and the 
friction drag is vanishingly small. 


Let us return to the question of the body with the smallest total drag cw. The racing car 
build in 1938 and shown in Figure 1.3 with a c„ value of 0.17 is an example of the ideal 
geometry. It is streamline bodies, as shown in Figure 2.49, which have the smallest drag. 
However, in Chapter 2.4.5 we will see that even these small drag coefficients can be decreased 
even further by suitable influencing of the wall shear stress 7,,. 


— CD QO 
Fig. 2.49 : Streamlined bodies in a free 
Co = C—C O stream (Zeppelin) and close to the ground 


| — (automobile) 
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2.3.3 Compressible Flows 


Compressible flow is characterized by a quantity called the compressibility K 


K= relative change in volume _ dV oi | (2.71) 
change in pressure required V dp 


Since the pressure change dp > 0 for a volume change dV < 0, the definition for K gets 
an additional minus sign, so that K itself is positive. Its value for water, for example, is 
Kino = 5-10-° 1/bar. For gases at constant temperature, the Boyle-Mariotte law holds: 


m 
= t. = 2.72 
p= const: ç (2.72) 
m dV m 
V = const.. — > — = -const.: — , 
p d p* 
with (1/V) = (p/(m- const.) it follows that K is 
dV 1 
=- Z = const. Z. p K=- 
dp V p? m- const p 


The value for air at p = 1 bar is Kair = (1/p) = 1 (1/bar). A comparison between air and 
water yields 


K air 


H20 


= 20000 


Therefore air is about 20000 times more compressible than water. We have already used this 
fact by treating water flows generally as incompressible flows, and gas flows at higher flow 
speeds as compressible flows. As well as the characteristic numbers of the previous chapter, 
we now also introduce the Mach number M 

C flow velocity 


M=- = 


“Y 2.73 
a speed of sound ( ) 


as an additional dimensionless characteristic number. The speed of sound a corresponds 
to the velocity of expansion of small disturbances in the variables of state (e.g. pressure 
disturbances dp) in a compressible medium at rest (Figure 2.50). The velocity of sound is a 
signal velocity with which disturbances in the flow field are transmitted. The gas through 
which the sound wave has passed has a pressure perturbation dp, a density perturbation dp 
and a perturbation in the velocity de. 


For an observer moving with —a, the sound wave is at rest and behind the sound wave he 
sees the velocity dc — a. If we restrict ourselves to inviscid outer flows, we can write down 
the continuity equation for the sound wave at rest 


d 
m=p:c-A=const. => (pt+dp)-(-at+dc)-A=-p-a-A = a. 
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and the Bernoulli equation becomes 


9 Pd _ de)? ptdp q —q)2 P d 
Z+ f P= const. => ae | dp _ | a) +f 2 
2 o P 2 0 p 2 0 


_ aP 
p 





=> a-de 
The speed of sound a is therefore coupled with pressure and density changes in the medium. 


Small disturbances expand without losses, i.e. isentropically, and so for the square of the 
speed of sound we can write 
o 
= (2 
2). 


This corresponds to the equation of definition (2.9). Using the equation of isentropic change 


of state 
p p\* 
—=([(— 2.74 
Pı 4) ) 
it follows that 
ðp | (2) 1. (2) bP p 
= = pi k.: — .— = k.: —  ——— = k: — - — = k: —=— ; 
Op pı pı pı (2) p M p 
Pl 
slightly perturbed gas at rest propagation 
quantity (one-dimensional, steady) 
a 
p+dp — p 
p +% p 
dc c=0 
sound wave 
slightly perturbed unperturbed reference frame moving with -a 
quantity quantity 
p + dp -a p 
p +p — p 
-a + dc -a 


Fig. 2.50 : Expansion of a sound wave in 
sound wave at rest a reference frame at rest and moving with 
the sound wave 
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and with the ideal gas equation (2.8) 





(2.75) 
with the general gas constant R = 8.314 J/(mol- K) and the molar mass M {=}. We can 
write down the following important relations for the speed of sound a 

anvVvVT , an a (2.76) 
M 
The numerical values for air are 
J m km 
k=ls , R=28i = ; T293 15K = a=34320— = 12355 — 
kg: K S h 


Sound waves are all around us in both natural and technical areas. One impressive example 
is that of the crack of a whip. Figure 2.51 shows four snapshots of the end of a whip. In 1 
the end of the whip is just about to snap around. In 2, the whip opens up and the crack 
wave S is formed. This can no longer be considered to be a small disturbance. The sound 





Fig. 2.51 : Crack of a whip 
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wave is built up to a shock wave, which we hear as the loud crack. In the third and fourth 
pictures the wave spreads out into the compressible surrounding air. 


Consider the sound waves in Figure 2.52 which spread out from a source at rest or a moving 
source (such as the crack of a whip). In the case of the source which is at rest, the sound 
waves expand as concentric spherical waves. If the sound source is moving with a velocity 
U. smaller than the speed of sound as( Mæ < 1), the spherical waves are denser upstream. 
An external observer first hears a higher frequency (higher tone) and then, after the sound 
source has moved past, a lower frequency (lower tone). If the sound source moves with a 
velocity U” which is greater than the speed of sound as( Mæ > 1), the sound waves remain 
inside a characteristic cone, known as the Mach cone, with a cone angle of sin(@) = Ago /Us. 
If the source of the sound is a supersonic airplane, this Mach cone is built up to a shock 
wave (head shock wave), whose pressure distribution on the ground is sketched in Figure 
2.53. The shock wave generates a pressure jump Ap on the ground, which we hear as a 
sonic boom. In order to be able to attain the unperturbed thermodynamic state of the air 
Yoo behind the airplane again, a further shock is needed, the tail wave. This reverses the 
pressure increase due to the head wave. This explains the double sonic boom heard on the 
ground when a supersonic airplane flies past. 


Mach Number Regimes 
As well as characterizing viscous flows with the Reynolds number Rez; flows with heat 
transport with the Prandtl number Pr; and the effect of gravity with the Froude number 


Fr;, the Mach number M,, allows us to differentiate between regions of incompressible and 
compressible flow. An incompressible subsonic flow with 0p/0s < Oc/Os is present when 


Mx <1 subsonic flow, incompressible (flow past an automobile) , 
compressible subsonic flows are found in the regime 


0.2 < Mæ < 1 compressible subsonic flow (ICE or TGV trains) , 





u 70 u<a, 
a At 
UA t 
sound source moving moving 
at rest sound source sound source 


Fig. 2.52 : Wave propagation from a source 
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transonic flows with p/s ~ 0c/Os are found for 
Mel transonic flow (civil airplane) 
and supersonic flows with Op/ös > c/ðs for 
Mæ > 1 supersonic flow (supersonic airplane Concorde) 
Hypersonic flows occur for 


Mæ > 1 hypersonic flow (re-entry vehicle, Space Shuttle) 


Mo> 1 





front shock wave tail shock wave 


reflected 
shock waves 


Ap 


Fig. 2.53 : Supersonic flight and pressure distribution at ground 
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The thermodynamic equations of state for ideal gases are no longer valid in the hypersonic 
flow regime. In this Mach number regime, the chemical reactions of hot air have to be taken 
into account. These are treated in the textbook Aerothermodynamik, OERTEL 1994. For 
example, for the Mach number M,, = 10, the effect of compressibility dominates and 


1 1 
op 10 Oc 


= n 100-2 


p Os c Os 


Stream Filament Theory of Compressible Flows 


The derivation of one-dimensional stream filament theory for compressible flows is based on 
the Euler equation (2.40). In what follows, we consider a steady shear flow, such as inviscid 
outer flow or inviscid core flow of a nozzle. For the shear flow dz = 0, and the Euler equation 
is written in the stream filament coordinate s as 


„de Ady dp dp ak dey 
ds pds  p dp ds p ds l ? 
1 de___} 1 &p 
cds MÈ p ds ’ 


, (2.77) 





with (1/p)- (dp/ds) the relative change in density and (1/¢) - (dc/ds) the relative change in 
velocity. 


In the subsonic regime M2, < 1, and so, as already mentioned, the relative change in density 
in subsonic flows is much smaller than the relative change in velocity, and for small Mach 
number can frequently be completely neglected. 


In the supersonic region the opposite is true. Because M2, >> 1, the relative change in density 
is very much larger than the relative change in velocity. If a supersonic flow is accelerated, 
de/ds > 0, then, because of the factor —M2,, this acceleration leads to a considerable 
reduction in the density of the medium, dp/ds < 0. Therefore supersonic flows need space. 
Because of the continuity equation, the strong relative reduction in density in an accelerated 
supersonic flow means that the cross-section A increases along s. 


In transonic flows M2, = 1, and all changes, both density changes and relative velocity 
changes, are of the same order of magnitude. Again the integral of the Euler equation leads 
to the Bernoulli equation for compressible flow. If we consider the integral along the 
stream filament s from the position 1 to the position 2, neglecting gravity with z, = z2, we 


obtain 
p2 d 
-Aa Lo 
Pl p 


The equations of isentropic change of state (2.74) hold for the change in the variables of 
state (but not for shock waves!), 


K 4 4 
z- (2) > (Bot N 
Pı pı p P pı pı 
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P2 dp pr [P ı pr K nat |? pr K s1 Kol 
J 2.) pk - dp = L] | H, m" -= pi" | => 
p P pi pı 


[2-5 (2-2) 
pı pP K—1 p2 pı 


So the Bernoulli equation for compressible flows reads 





(2.79) 





Using the equation of state for the ideal gas (p/p) = R- T = (cp — c) - T and the isentrop 
exponent x = (Cy/C,) it follows that 


K p c 1 








1 1 
25o NT = Cp T +3: Ë = const. , (2.80) 


(2.81) 





This corresponds to the energy equation (2.52), without taking the heat flux and gravity 
into account. 


Determining the constants of the Bernoulli equation is carried out using the reservoir 
values, or the so-called critical values. 


The reservoir values po, po, ao, To with c = 0 satisfy equation (2.79) 











1, @ ao ofl ‚a 1 ao 
2 tI k] "153 FE) kI 
in a? u 1 
a 1+5 M? 


with a =«K-R-T and a2 =K-R- Tp it follows that 
T 1 


h 1+ M 


(2.82) 
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T is always smaller than 7%, since M2, > 0 holds at all times. The isentropic relation 


e- (F) 
po To 


leads to an expression for the reservoir pressure po 


p 1 


— = n ,. (2.83) 
K—l KL 
Po (1+ ye 
Since M2, > 0, p is smaller than pọ. Using the isentropic relation 
(m) 
Po To 
it follows that the reservoir pressure po satisfies 
1 
P —. (2.84) 


Pp (145.92) 


Equation (2.82) can also be used to determine the reservoir temperature To at the stagnation 
point of a projectile. If we assume a flow temperature of T = 300 K, in the stagnation 
point (c = 0), for a Mach number of M@, = 2, we find the stagnation point temperature 
To = 540 K. The stagnation point of the supersonic airplane Concorde therefore heats up 
during the flights. For Mæ = 5, the stagnation point temperature is already Tọ = 1,800 K. 
At such high temperatures, however, the assumption of isentropic change of state and the 
ideal gas law are not longer guaranteed. 


In order to determine the constants in the Bernoulli equation, the critical values can also be 
used (index x). Critical values are those values found in a flow when the velocity of sound 
M = 1 has just been reached: 








p_M=1)=p* ,„ T(M=1)=T* , pM=1)=p* , 
aM=1)=a , (M=1)=C=a 
Therefore 
2 *2 1 1 k+1 
2 __ *2 __ *2 zZ — *2 
at TROI (; =) 2. (k— I) 
or 
1 1 
z eto T= a to T" , 
with 
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Therefore there is a relation between the reservoir values (index 0) and the critical values 
(index x). To find this relation, we set the Mach number in equations (2.82) and (2.84) to 
M = 1, give variable quantities the index x, and leave the reservoir values as they are. We 
obtain 


(2.85) 





In particular, for air with the value x = 1.4 we find 


* * * 


T 
—_-083 , 2 =0.634 L — 0.528 
To Po Po 
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Stream Filament Theory with Variable Cross-Section A(s) 
For a variable cross-section A(s) the continuity equation reads 
m = p(s) - c(s) : Als) = const. 
Taking the logarithm of the continuity equation we find 
In(p(s) - e(s) - A(s)) = In(p(s)) + In(e(s)) + In(A(s)) = In{const.) , 
and differentiation d/ds leads to 


1 dp 1 de 1 dA _y 
p ds c ds A ds 


Using the Euler equation (2.77), the density term can be eliminated from the continuity 
equation and we obtain 


(2.86) 





Equation (2.86) allows one to determine how the cross-section A(s) of a nozzle must be 
shaped to accelerate a gas continuously from subsonic Mach numbers M < 1 to supersonic 
Mach numbers M > 1 (Figure 2.54). Continuous acceleration requires that de/ds > 0. If the 
Mach number is less than one, M < 1, a narrowing of the cross-section is needed: dA/ds < 0. 
If the Mach number is greater than one M < 1, a widening of the cross-section is required: 
dA/ds > 0, in order to accelerate the gas. For the Mach number M = 1, the differential 
equation (2.86) has a singularity. In order to ensure that de/ds > 0, dA/ds = 0 has to hold. 
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If one wishes to accelerate a flow continuously from the subsonic region to the supersonic 
region, the nozzle required must first have a narrowing in the cross-section, and then, down- 
stream from the narrowest cross-section, a further widening of the cross section. Such a 
nozzle is shown in Figure 2.54. It is called a Laval nozzle. 


At the narrowest point in the cross-section, the critical values (index *) introduced in equa- 
tion (2.85) are found at Mach number M = 1. The divergent part of the nozzle in the 
supersonic region can be explained by considering that the relative decrease in density in 
the supersonic region is much greater than the relative increase in velocity. For this reason, 
conservation of the constant mass flux m = p- c- A = const. requires that the cross-section 
A(s) increases along s. 


In what follows we will derive the differential equation which relates the relative change in 
cross-section, (1/A)-(dA/ds), to the relative change in Mach number (1/M)- (dM - ds). The 
logarithm of the defining equation for the Mach number c = M - a yields 


In(c) = In(M) + In(a) 


Differentiation d/ds leads to 
1 de 1 dM 1 da 
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M<1 M=1 M> 1 Fig. 2.54 : Laval nozzle 


86 2. Fundamentals of Fluid Mechanics 


Taking the logarithm of a? = «- (p/p) yields 2-In(a) = In(«) + In(p) — In(p). Differentiation 
d/ds leads to 


2 da 1 dp 1 dp 


a ds p ds p ds 
In the next step, the expression dp/ds has to be related to dp/ds: 


dp dp dp p dp 
2 _—_ dp — a’? -d L =a- kp. E.L 
a dp => p=a -dp => qe Fe Kee J; = 
1 dp x dp 
p ds p ds 
We find that 
2 da dp 
2. -(£- 1:2. 
a ds («= 1) s o’ 


and with the Euler equation it follows that 
1 da «-1 —M? de 


ad 2 C ds 





Inserting this equation into equation (2.87), taking equation (2.86) into account delivers 


1 1 dA 1 dM _ (s= 1I)(-M?) 1 1 dA 
M2—1 A ds M ds 2 M2—-1 A ds ” 


(2.88) 


1 dA (a) 


A ds M?-1 





This is a first order ordinary differential equation to determine M (s) for a given cross-section 
A(s). Using the boundary condition M = M* = 1 for A = Amin = A* with M* = 1 the 
solution reads 


(2.89) 





Equation (2.89) gives the Mach number as an implicit function of a given cross-section A(s), 
if the velocity is the velocity of sound at the narrowest point A*. In this case, the mass flux 
m through the nozzle can be determined as a function of the critical values 


* * 


m=p-c:-A=p-c-A*=p"-a*- A* =const. 


For the solution curve discussion of equation (2.89), we consider the direction field of the 
ordinary differential equation (2.88). First of all we solve equation (2.88) for dM/ds 


dM _,,.1 dA 1+ 44+. M? dM 
ds A ds 


Wo Fe M8) = FM A, 5) 
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The derivative of the Mach number W’(s) is therefore a function f of the Mach number M (s), 
of the given cross-section A(s) and of the coordinate s. The relation M’(s) = f(M, A,s) 
assigns a direction in the (s, M) plane to every point (s, M). 


A direction element to be specially defined is found at the narrowest cross-section Amin of 
Figure 2.55. For M #1, with dA/ds = 0, we find 


dM 
a 
ds 


? 


i.e. horizontal tangents . For M = 1, as long as dA/ds # 0, vertical tangents are found, with 


dM | 
ds 


The singular point at the narrowest cross-section Amin with dA/ds = 0 and with Mach 
number 1 is a saddle point. The singular point has no uniquely defined direction assigned 
to it. There are two possible directions of propagation. These three limiting cases can be 
used to draw the mathematically possible solution curves to equation (2.89) in Figure 2.55. 


Not all the solution curves are physically relevant for the (assumed) continuous acceleration 
in the Laval nozzle. The region where backflow occurs, for example, is not physically relevant, 
leaving those relevant solution curves as shown in Figure 2.56. The solution curve that 
ultimately is found in the Laval nozzle depends on ratio of the pressure p at the end of 
the nozzle to the reservoir pressure po. Depending on the application of either high back 
pressure pı or low back pressure pr, we obtain different flow forms, to be treated in what 
follows. 





Fig. 2.55 : Direction field of the Laval nozzle differential equation 
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For a high back pressure p4, we obtain pure subsonic flow M < 1 through the Laval nozzle. 
In the region where the cross-section narrows, the flow is accelerated (nozzle). Where the 
cross-section widens, the flow is then decelerated for M < 1. In this case the Laval nozzle 
acts as a diffusor. 


If the back pressure is pg, the Mach number 1 is just reached at the narrowest cross-section, 
and the critical values are obtained (index *). In the region where the cross-section of the 
Laval nozzle widens, a subsonic flow is again found, and the flow is decelerated. 


If the back pressure exceeds this critical value pp, as for the value pp, acceleration occurs for 
supersonic values M > 1, but continual flow through the Laval nozzle is no longer possible. 
A shock occurs in the subsonic region, leading to a discontinuity in the flow quantities. 
At the position s, the solution curve jumps from supersonic M > 1 to subsonic M < 1 
behavior. 


If the back pressure at the end of the nozzle is lowered to the value pp, the shock moves to 
the end of the nozzle. 


It is only for the back pressure pr that we can speak of an ideal Laval nozzle. The continuous 
acceleration of the flow follows the upper solution curve in Figure 2.56 from the subsonic 
regime M < 1 to the supersonic regime M > 1. At the end of the nozzle, as sketched in 
Figure 2.57, a free jet without a shock is found, and the pressure is the surrounding pressure 


Doo: 


For back pressures between pp and pr, oblique shocks are found at the end of the nozzle, 
followed by so-called expansion fans. This flow form of oblique shocks followed by expansion 
fans is then to be found periodically in the free jet, so that a characteristic node-structure 
appears. Such a supersonic free jet is used in a blowtorch to cut through metal. 


mT 
Oo 
m 





Pp M<i "Amin p Fig. 2.56 : Dependence of Mach number 
Po in the Laval nozzle on the back pressure p 
To at the outlet of the nozzle 
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If the back pressure at the end of the nozzle is further reduced to pr, the oblique shocks 
vanish, and pure expansion flow occurs at the end of the nozzle, seen as a free jet bell. This 
can be seen at high altitudes after a rocket starts. 


Figure 2.58 shows the mass flux density in the Laval nozzle. The mass flux density is the 
ratio between the mass flux m and the cross-sectional area of the flow A 

m 

a 
In the supersonic case, for the Laval nozzle with critical values at the narrowest cross-section 
A* = Amin we have 
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Since the cross-section A in a Laval nozzle is larger than A* everywhere except at the 
narrowest point, where Amin = A*, we have 
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Fig. 2.57 : Dependence of the flow structures at the end of the Laval nozzle on the back 
pressure p 
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Fig. 2.58 : Mach number and mass flux 
density in a Laval nozzle 





Therefore at the narrowest point in the Laval nozzle where A,,j;, = A*, the dimensionless 
mass flux density (p- c)/(p* - c*) has its maximum value (p- c)/(p* + c*) = 1. 


Shock Wave 


A shock wave is quite generally a sharp change in the flow variables velocity v, pressure 
p, density p and temperature T. These changes occur in an extremely thin layer of the gas 
with a size of the order of magnitude of several mean free paths. The mean free path is the 
mean distance which a molecule or atom moves between two collisions with other molecules. 
For air, the mean free path A in normal conditions is \ = 1077 m. In regions of this order 
of magnitude very large gradients of the variables of state occur, so that it is permitted 
to model the shock wave by a discontinuous change within the framework of continuum 
mechanics. A compression shock wave, which the expression shock wave generally refers to, 
is a discontinuous increase in the density p in the shock region. As well as the density, the 
temperature T and the pressure p also increase, while the magnitude of the velocity | v | 
decreases. 


A shock wave can in principle only occur in a supersonic flow region. In the special case of a 
normal shock wave, where the free stream direction and the shock front form a right angle, 
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the flow progresses from a supersonic region to a subsonic region. In the case of an oblique 
shock wave, seen for example in the Mach cone of the flow past the supersonic airplane 
Concorde, the free stream and the shock wave front form the Mach angle a, already met at 
the start of this chapter. In this case the flow can also progress from supersonic to supersonic 
flow, whereby the supersonic velocity after the shock wave must be smaller than that before 
the shock wave. 


The left-hand side of Figure 2.59 shows the velocity ratios on a schematic cross-section 
through a wing. ‘The supersonic regime on the wing is denoted by the Mach number M > 1. 
This region ends downstream with the shock wave, and here the flow is subsonic, with 
M < 1. This shock wave is slightly curved, and is almost perpendicular at the point where 
it touches the boundary layer. We will now write down the shock wave equations for such a 
shock wave. The same equations hold for the shock wave in the Laval nozzle. 


We assume, quite generally, a steady, inviscid supersonic flow. This is characterized by the 
values given for cı, pı, p1 and T;. Using the speed of sound (2.75), a; = W/k-pi/p1, we 
determine the Mach number of the free stream Mı = cı/aı. Here & denotes the ratio of the 
specific heats c,/c,. On passing through the shock wave these values undergo discontinuous 
changes in the direction of the normal to the shock wave surface. We consider the flow 
variables c2, po, po and T> downstream from the shock wave surface. The velocity cz is then 
smaller than the free stream velocity cı, while the other variables of state increase. On the 
right-hand side of Figure 2.59 this is shown by a smaller vector for ca after the shock wave. 
The change of state through the normal shock wave can be described using the conservation 
of mass, momentum and energy equations for a one-dimensional, steady, inviscid flow be- 
fore and after the shock wave. We start with the equations from the one-dimensional theory: 


Mass: Pi Cr = fo+ C2 (2.90) 
Momentum: Pit pr: E = P2 + 2° & (2.91) 
l > l 2 
Energy: hı + 5 a= ho + 3 (2.92) 
shock 


sonic line ae 





Fig. 2.59 : Shock wave on a transonic foil and change of state through a normal shock wave 
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For the enthalpy h we have the caloric equation of state: 
h=c,-T=c+2=c:T+2 
p p 
Solving the basic equations (2.90) - (2.92) for the four unknown quantities after the shock 
wave Co, P2, p2 and To, we obtain the shock wave equations. 


Taking the thermal equation of state for ideal gases p/p = R- T into account, we can write 
the dependence of the enthalpy on the following quantities: 


1 y 2 
n=a g E C a)? K P a 
Rp p Cp — Cy p k—1 p 











Thus the balance of energy (2.92) reads: 


k Pi lo. k pR 1 %3 
k-i m 2 UTRI p 2” 








Together with the conservation equations for mass (2.90) and momentum (2.91), we obtain 
a system of three algebraic equations to determine the three desired quantities co, pa and Pa 
after the shock wave. The temperature To can then be computed with the thermal equation 
of state using pọ and p>. Assuming the initial values cı, pı and pı are given, the system of 
equations can then be solved for the desired variables. We obtain 


e m Jt 

a a -A G-3) 
P2 _ l 2 

p (1+2 (22-1) 


? 


For given initial values in front of the shock wave, the system of equations yields two 
solutions. The upper solution with the value 1 is the identity solution, for the case where no 
shock wave occurs. The lower solution is the desired shock wave solution. Using the velocity 
of sound aı = \/K : pı/pı and the Mach number Mı = c/a we can bring the shock wave 
equations to a form where the only parameter on the right-hand side is the Mach number 
M,>1 
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The shock wave equations (2.93) - (2.95) yield the dependence of the values after the normal 
shock wave on the free stream Mach number. Whereas the pressure and the temperature 
after the shock wave can increase arbitrarily with increasing free stream Mach number, the 
density ratio p2/p1 for Mı — co tends towards the ratio («+1)/(K—1). For air with « = 1.4, 
the density after the shock wave can reach at most six times the free stream density. However 
this estimate is reached with the assumption of an ideal gas. 


We now determine a relationship between pa and ps» after the shock wave and to do this 
eliminate co from the equations (2.90) - (2.92). After several steps, we obtain a relation 
which describes a symmetric hyperbola in the (p1/p2, p2/pı) plane. This can be used to 
follow the thermodynamically possible changes in the variables of state pı and pı through 
the shock wave. This hyperbola is called the Hugoniot curve and reads 


p2 č k=l RIT bp 
Do Rel Boel (2.97) 





A further relation is obtained when we derive a relation for p/p, as a function of p1/Ppa 
from only the conservation of mass and momentum equations (2.90) and (2.91), without 
using the balance of energy. This then leads to the kinematically possible changes of state, 
described by an equation for a straight line. This line is called the Rayleigh line and reads 


Bise Mp (2-1) (2.98) 
Pı p2 

The Rayleigh line has the slope —K- M?, and meets the Hugoniot curve at two points, at 
the identity with pa = pı andpa = pı and at the shock wave solution after the shock wave 
(Figure 2.60). 







Hugoniot curve 


state after shock 


Rayleigh line 


kad Fig. 2.60 : Hugoniot diagram 
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The areas in the Hugoniot diagram can be interpreted as energies. The area below the 
Rayleigh line A’B’C'D is the internal energy e of the shock wave 


= 1 
oo 5. (2-1). (1-2) 41. (1-4) 
T 2 \pı p2 p2 
EERS a 


ABCD A'B'CD 





The triangular area AC'D above the Rayleigh line represents the kinetic energy c3/2 


1 
z (2-1) 6-2) | 
2 \Ppı p2 
ee 
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so that the total area A’B’C'D represents the increase of the total energy in the shock wave. 


In front of a blunt body in a supersonic flow M, > 1, a head shock wave as shown in Figure 
2.61 occurs. Around the stagnation streamline, the head shock wave can be approximately 


Ti, M} Ts 


M}= 2.9 





Fig. 2.61 : Leading shock wave of a sphere 
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considered as a normal shock wave. The temperature in the stagnation point Tə can be 
computed with the energy equation (2.92) and the caloric equation of state h = c,-T 


c2 


Cp Ts = Cp Ti+ > 
With M, = c,/a,, a? =K«-R-T, Cp —Cy = Rand k = c,/cy the stagnation point temperature 
Ts is 


Is k—1 
Ti 2 





.M? . (2.99) 


We have already computed Tg = 540 K for a supersonic flight with Mı = 2. For a hypersonic 
flight with M, = 10, a stagnation point temperature of Ts = 6, 300 K occurs, meaning that 
heat-shielding materials such as ceramic tiles are necessary for heat protection. Since the 
heat transfer depends on the radius of curvature, and is relatively low for large radii, i.e. 
for blunt bodies, re-entry vehicles, such as those used for the space shuttle, we designed to 
be blunt. Figure 2.62 shows the space shuttle in a supersonic wind tunnel. The head shock 
wave is an almost normal shock wave close to the stagnation streamline, which passes over 
to the oblique shock wave of the head shock wave. We have already mentioned that the 
downstream Mach number after the oblique shock wave can be M > 1. Therefore the wing 
of the space shuttle may also be in a supersonic free stream, leading to a second head shock 
wave in front of the wing. 


Figure 2.63 summarizes the possible flow forms from subsonic to supersonic free stream past 
a wing profile. In a subsonic free stream, with a Mach number smaller than Mæ = 0.75, the 
acceleration does not reach supersonic Mach numbers M > 1 at any point on the profile, 
and the flow is purely subsonic. In the case of the transonic Mach number Ma = 0.81 we 
obtain the supersonic region on the profile already discussed in Chapter 1.1 (Figure 1.18). 
At the end of the profile here is an almost normal shock wave. For subsonic Mach numbers 
larger than 0.85, a shock wave forms on the lower side of the profile, and for Mach numbers 
close to 1, this, together with the upper shock wave, passes over to an oblique shock wave 
of the tail wave. For supersonic Mach numbers M,, > 1, a separated head shock wave first 
appears in front of the profile. For the supersonic flight Mach number M,, = 2, an attached 
oblique shock wave appears as the head shock wave, which, together with the tail wave, 
leads to the double sonic boom of the supersonic airplane discussed in Figure 2.53. 


Fig. 2.62 : Head shock waves in front of 
the re-entry vehicle Space Shuttle, Mı = 3 
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Ma, M; 


0.75 


0.81 


0.98 


1.4 
Fig. 2.63 : Flow forms past a wing profile 
from subsonic to supersonic free stream 





Oblique shock waves can be computed using the fundamental equations of the normal 
shock wave (2.90) - (2.92) and (2.93) - (2.96), as long as these are applied to the normal 
components of the velocities. Figure 2.64 shows the change of direction of the velocity 
vector € = (cn, Cı) through an oblique shock wave with normal component c„ and tangential 
component c;. With 


Cni =G- sin(a) , Gı=cı:cosle) , 


Cn2 = C& sin(a — 8) , C2 =C: cosla— p) , 


the fundamental equations of the oblique shock wave (2.90) - (2.92) can be written 


S 





Fig. 2.64 : Oblique shock wave 
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Mass: P1 Cn = P2'm2 ; (2.100) 


Momentum: pı +1: a =P tp Go. (2.101) 





Pi? Cn © &t,1 5 P2 T P2° Cn2° Ct2 ; 


1 1 
Energy: hı + 2° È = ho + z` GS (2.102) 


Equations (2.100) and (2.101) yield the tangential component 
C1 = C2. (2.103) 


Equation (2.102), with c? = c? + c? delivers 


hy ba C2, =hn ba, . (2.104) 
2" 2" 

Therefore, the shock wave equations for the normal shock wave hold for the normal compo- 
nents of the velocity in front of and behind the shock wave, with the additional condition 
that the tangential components c; , and c;» must be the same. Figure 2.65 shows the possible 
shock wave angles a for different free stream Mach numbers M,, and it is seen that beyond 
a certain limiting value Gq of the downstream flow angle 8, an oblique shock wave is no 
longer possible. For 6 > Gg the separated head shock wave already discussed occurs. 


Unsteady shock waves are generated by a shock tube. A shock tube consists of a high- 
pressure section and an low-pressure section, separated by a membrane. If the high-pressure 
section is filled with propellant gas at over-pressure until the membrane bursts, an unsteady 
shock wave in the low-pressure section of the shock tube, filled with the test gas, moves 


OL . 
attached front shock wave ı separated 


90° front shock wave 





o Ba =50 B Fig. 2.65 : Shock wave angle a of oblique 
shock waves 
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shock 


expansion compression rest 


I 
high pressure «—————————— low pressure ——————> 











expansion compression 


X 


Fig. 2.66 : Shock tube, position-time diagram of a shock wave and an expansion wave 


with constant velocity c, (see the diagram in Figure 2.66). In the high-pressure section the 
corresponding expansion wave moves. If we move with constant shock wave velocity c, with 
the shock wave, we can compute the change of state through the unsteady shock wave with 
the fundamental equations (2.90) - (2.92) and (2.93) - (2.96) of the normal shock wave. 


l > 
Cq=-Cy , C2 = C2 — Ciha + 5° in 


2.4 Technical Flows 
2.4.1 Turbulent Flows 


Most flows which occur in nature and in technology are turbulent at high enough Reynolds 
numbers. In contrast to the laminar flows which we have discussed until now, turbulent flows 
are characterized by fluctuations in the flow variables, leading to an additional transverse 
exchange of momentum and energy. This results in more full time-averaged velocity profiles 
compared to the laminar profiles found in boundary layers, channels and pipes. 


Figure 2.67 shows the laminar velocity profile already discussed compared with the profiles 
of turbulent boundary layer and pipe flows which occur when a so-called critical Reynolds 
number Re. is exceeded. If a colored filament is added to a steady laminar flow, as in Figure 
2.68, a straight streakline will be seen, as already seen in Chapter 2.3.1. In the turbulent 
flow the colored filament breaks up because of the fluctuations superimposed on the flow 
and the consequent transverse exchange of momentum. 


2.4 Technical Flows 99 


u(r) U(x,r,@,t) 


p 





L laminar turbulent laminar turbulent 
boundary layer flow pipe flow 


Fig. 2.67 : Laminar and turbulent velocity profiles in boundary layer and pipe flow 


Turbulent flows are fundamentally three-dimensional and time dependent. We now leave the 
area of one-dimensional stream filament theory and return to denoting the flow variables as 


veya PEUR s PE YZ 1): 


The mathematical description of turbulent flows is derived from the experimental results 
in Figure 2.68. Reynolds was able to conclude from his experiments that the flow variables, 
such as the u component of the velocity, can be represented as a superposition of the time- 
averaged velocities u(xz,y,z) and the additional fluctuations u’(x, y, z,t). The Reynolds 
ansatz for turbulent flows is written (see Figure 2.69) 


V(2,y,2,t) = V(a,y,2)+ 0 (ay, 2,8). (2.105) 


laminar 










——————————————a eee 


color turbulent 





u(x,r,p,t) 





Fig. 2.68 : Reynolds experiment: laminar and turbulent pipe flow , REYNOLDS 1883 


100 2. Fundamentals of Fluid Mechanics 


The definition of the time-average at a fixed position for the example of the velocity com- 
ponent u reads 


1 T 
ü= =. / ulz, y, z, t) dt . (2.106) 
T Jo 

T is a time interval large enough that any increase in T does not lead to a further change 
in the time-averaged value %. From the definition of the time-average, we can see that the 
time-average of the fluctuating quantities must vanish, i.e. the velocity fluctuations satisfy 


© 


w=0 , v=0 , w= 


The proof of this for the u component of the velocity reads: 


1 fT ı fT 1 f7 1 f 


SI 


_ _ 1 T 
> uw=uUutwv > 10-5: | u -dt 
0 
The dimensionless turbulence intensity Tu is used to characterize turbulent flows. In the 
numerator of this quantity is the square root of the time-averaged square of the fluctuating 
variables, and in the denominator the time-averaged flow velocity at a certain position. The 
turbulence intensity for the velocity component u in the main direction of flow x is 


ww)? 


U 





Tu = 


Since turbulent flows are three-dimensional, the three-dimensional generalization of the 
turbulence intensity at a certain position in the flow field follows as 





1 (WP +P + we) yh (WP + wP+ WP) 
Pus 2 eae (2.107) 


Fig. 2.69 : Reynolds ansatz for the u com- 
ponent of the velocity 
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Because of the fluctuating motions wu’, v’ and w’ in a turbulent flow, there is an additional 
contribution to the flow drag. However, this additional contribution has nothing to do 
with the molecular viscosity jz, but rather is due to the additional transverse and longitudinal 
exchange processes which occur in a turbulent flow. In what follows we will look at the 
mathematical treatment of these processes. 


Our starting point is the Navier-Stokes equation (2.59). In passing from the stream filament 
coordinate system, s and n, to a Cartesian coordinate system (x, y, 2), the velocity c along 
the stream filament is replaced by the variable u, s by z and n by z. We obtain 


ðu ðu 1 Op o?u dz 


a t or p ðr OPT Tr 
This equation is also valid in principle for turbulent flows but has to be supplemented by 
the second and third Navier-Stokes equations for the u and v components respectively. For 
an incompressible flow with p = const., the turbulent flow quantities appearing in equation 
(2.108) are the velocity component u and the pressure p. Using the Reynolds ansatz (2.105) 
for u and p we find 


(2.108) 


Taw), D@t I OD), PUFA E 
Ot tatu) ðr p Ox ty 02? I` Tx 
Taking the computation rules for the time-average into account, this leads to 
_ oü +u’) 1 op ou dz 
I)e n = — _—_ . — 
(u+ u’) Ox p Ox tv: Oz? dx 


Note here that 0(u@ + u’)/Ot = 0 is only valid for flows whose time-average is steady. Such 
flows are called quasi-steady turbulent flows. The time-average of the non-linear inertia 
terms on the left-hand side of the equation is considered separately. We have 
(mtu) outu)  _ a t, OW a , w y. Ou ðu , Ou 
u+u’)-———~ =%-— +t - — u: — tu. — 
Ox Ox Ox Ox Ox Ox Ox 


In particular for the term wu’ - (Ou! /Ox): 
Ow 1 ft du 1 [708 (wP ð (1 [f WP 
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“ax T J "dr T J Ox ( 2 ) Ox (7 J 2 ) 
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Therefore the time-averaged Navier-Stokes equation reads 


(un? > 2, 
el -i2 Oru dz 
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If we multiply this equation with the constant density p and write the pressure and gravity 
terms on the left-hand side, we find 
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On the right-hand side of the equation are the terms responsible for the drag in the flow. 
As well as the shear stress 7,, due to the viscosity, a turbulent flow also has an additional 
drag term due to the velocity fluctuations, denoted here by the index ’ as 7/.,. In general, 
the additionally appearing parts of the stress 7’ in turbulent flows are called Reynolds 
apparent normal and shear stresses, since they are caused by the turbulent longitudinal 
and transverse exchange and not by the molecular viscosity u. 


The lower double indices in the stress variable 7 follow the same conventions used in solid 
state physics. The first index denotes the normal to cross-section and the second index 
denotes the direction in which the force on this cross-section acts. 


In the general three-dimensional case, 7’ is a stress tensor with 9 components, consisting of 
6 apparent shear stresses and 3 apparent normal stresses (trace of the shear stress tensor). 
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Because of the balance of moments, equivalent components of the stress tensor are equal, 
i.e. we have Tay = Tyg Or —p : W -wW = —p : w -u ete. The time-averaged products of the 
fluctuation variables and thus the components of the stress tensor 7’ are unknown and have 


to be described using model equations. 








Starting out from the Newtonian ansatz of laminar flows, Boussinesq made the assumption 
that the unknown fluctuation terms can be reduced to the known time-averaged quantities 
of the basic flow by introducing an unknown proportionality factor p,, called the turbulent 
viscosity. Using the Boussinesq assumption the following relations are found 


, — du ðu du 
Tey TTP U U = h: ar Br = M2 ’ 
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Here wu, is a function which has to be determined, and is not, like the molecular viscosity u, 
a material constant. 


One possible ansatz to determine jp; is the Prandtl mixing length ansatz. Figure 2.70 
shows a turbulent two-dimensional boundary layer flow in the (x, z) plane. The Reynolds 
ansatz yields 


/ 


u=Uu(z)+u 
w= w 

If we move a fluid element with fluctuation velocity from the height zo to the height zo + J, 
the change in ù with (zo +1) > (z9) is found writing out the Taylor expansion: 
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Neglecting higher order terms it follows that 


du 
u(2zo) — U(2o +8) = —-1- T 





This lower velocity —I- (du/dz|,,) at the height 2, + ! was considered by Prandtl to be the 


velocity fluctuation 


du 
l pa ee 
u (zo +1) Ap 





20 
at 2 +l. For continuity reasons, it follows for w’ that: 
du 

dz 


The mixing length / is therefore the distance a fluid element will travel before it is completely 
mixed with its surroundings and looses its identity. ‘This allows us to consider the velocity 
fluctuations u’ and w’ in terms of the mixing length / and the time-averaged velocity profile 
u(z). The apparent shear stress 7/,, = —p-u’-w’ can then be computed as 


——= dü dü dü\“ 
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Since we assumed a two-dimensional turbulent boundary layer flow with w = 0, we also 
have (OW/Ox) = 0, and it follows for the Boussinesq assumption (2.111) that 


w=: 





du 
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This yields an equation to determine the desired quantity pu, since we have 
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flow element Fig. 2.70 : Prandtl mixing length ansatz 
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and therefore 


dü 
=p. P. | — , 2.113 
wo) (2.118) 


The mixing length l is still unknown. It has to be determined from experiments, leading to 
empirical approximation formulae for the computation of J. 


Following these fundamental considerations of turbulent flows, let us return to the turbulent 
plate boundary layer flow in Figure 2.67. The order of magnitude of the turbulent apparent 
viscosity fl, permits the turbulent plate boundary layer to be divided up into regimes (Figure 
2.71). Directly at the wall u, < ps. This is the viscous sublayer and is of particular technical 
importance in reducing drag with so-called riblets. These will be discussed at the end of this 
chapter. 


In the viscous sublayer the velocity fluctuations u’ and w’ are very small, and for the mixing 
length we have | — 0. The total shear stress 7;., in the turbulent flow under consideration 
reads 


= dü aye ant 
Tit SH: Z TPU -W 


dz 


Because of u’ - w ~ 0, it follows that the wall shear stress 7,, in the viscous sublayer satisfies 
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Fig. 2.71 : Division of the turbulent boundary layer flow 
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Separation of variables leads to an ordinary differential equation for the desired velocity 
profile 


Integration yields 


U 1 zZ 
| m=. f mdz , 
0 H Jo 


that is, a linear velocity distribution u(z) for a constant shear stress 7, 


u(z) = 7 sz]. (2.114) 


Extending this expression with the constant density p delivers 


Defining the so-called wall shear stress velocity u, as a new quantity, with ur = \/(Tw/p), 
we obtain 


(2.115) 





with the new dimensionless coordinate 2* = (u, : z)/v. 


In the region of wall turbulence, outside the viscous sublayer but still close to the wall, the 
wall shear stress is still constant: 7, = const.. Prandtl assumed that the wall shear stress 
can be written as follows, with the mixing length | = k- z as a linear function of z (k is a 
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This leads to a differential equation to determine ü(z): 
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Indefinite integration yields 
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Summarizing the last two terms to a new constant of integration C, we obtain a logarithmic 
velocity profile in the region of wall turbulence 


(2.116) 





The time-averaged velocity profiles (2.114) - (2.116) close to the wall are shown in Figure 
2.72. The viscous sublayer exists in the region 0 < zt < 5. After this comes the transition 
region 5 < zt < 30 and then the logarithmic region, with 30 < zt < 350. 


The laminar-turbulent transition for the plate boundary layer flow takes place at the critical 
Reynolds number 


Re. = 5-10° 


For a pipe flow, the critical Reynolds number is Re, = 2300. 


The laminar-turbulent transition leads to an increase of the friction drag cy, shown in Figure 
2.73 as dependent on the Reynolds number Rez, formed with the length of the plate L. The 
local friction coefficient cp(x) at the position x = L is 


laminar boundary layer flow 


(2.117) 


turbulent boundary layer flow 





The transition from the laminar to the turbulent boundary layer flow does not take place 
suddenly, but rather over a transition region. The dimensionless integral friction drag coef- 
ficient Ç can be calculated from the local drag coefficient cp(x). The integral friction drag 
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Fig. 2.72 : Turbulent boundary layer profile 
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coefficients are defined as the ratio of the wall friction force Fr to the product of the dy- 
namic pressure and the plate surface area A = L-b. Here b denotes the breadth of the plate 
perpendicular to the plane sketched and L the length along the plate. 


Ou TaT) 2 
D = . — ; SS ; Fr=b w -d = 
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For the integral friction drag coefficient Ç at a distance L from the leading edge of the plate 
we have 


laminar boundary layer flow |(2.118) 


turbulent boundary layer flow 





Therefore the friction drag of a plate in a laminar flow is smaller than the friction drag of 
a plate in a completely turbulent flow under otherwise identical conditions. Therefore 


Ge > C1 


The different thickening behavior of the boundary layer thickness 6 of laminar and turbulent 
boundary layer flows can be seen in Figure 2.74. The starting point is the relation (2.61) for 
a laminar boundary layer flow: 






turbulent 


transition Fig. 2.73 : Friction drag cy of the laminar 


and turbulent plate boundary layers 
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In the case of the laminar Blasius boundary layer, the proportionality factor is 5: 
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Multiplication with y Rez, yields 
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For a turbulent boundary layer flow we have 
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Multiplication with y Rer, delivers 
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By manipulating the turbulent wall shear stress 7/,, the friction coefficient cp of the turbulent 
boundary layer flow can be reduced. The idea for this comes from nature. Fast-swimming 
sharks (up to 45 km/h) have microscopically thin grooves on their scales in the direction 
of flow. Figure 2.75 shows a close-up of the longitudinal and transverse grooves on a single 
scale of a blue shark. 


It can be assumed that less friction occurs on surfaces with longitudinal grooves than on 
smooth surfaces. This knowledge is applied in technology in the production of foils with 
longitudinal grooves, so-called riblets, with a height zt = 500 and a separation of y? = 100. 
Such foils are applied to smooth surfaces to reduce their friction drag. 


As a result, the fluctuations in the transverse flow v’ and therefore the transverse exchange 
in the viscous sublayer of the boundary layer are prevented. The outcome is a reduction 


laminar 





Fig. 2.74 : Boundary-layer thickness 6 of 
5-10° Re, laminar and turbulent plate boundary lay- 
ers 
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Fig. 2.75 : Shark scales and riblet foil 


Riblet foil 


in the friction drag cy of 8 %. In the case of an airplane, the friction drag cy is more than 
50 %. Since riblet foil cannot be applied to all parts of the aircraft, the real potential of 
drag reduction is 3 %. It has been seen that a 1 % reduction in fuel used by an Airbus A340 
is made when 30 % of its surface is covered with riblet foil. The drag reducing foil can also 
be used in next-generation high-speed trains, as well as to reduce losses in pipe flows. 


2.4.2 Balance of Momentum 


The balance of momentum is carried out on a control volume V to determine directly the 
integral forces for given flow variables at the edge of the control volume V. The momentum 
dI of a mass element dm = p- dV is defined as the product of the mass element and the 
velocity vector V 


a 5 dl, u u 
dl=dm-V=p:-v-dV > d= | dl, | =dm-|[v]=p-{v |-dv 
dl, w w 


Since the density p can be time dependent in compressible flows, p = p(t), conservation of 
mass m = p(t) - V(t) = const. requires that the volume V must also be taken to be time 
dependent, i.e. V = V(t). The momentum of the total mass m, or of the total volume under 
consideration, is calculated by integrating the differential momentum dI of the mass element 
over the volume V(t): 


Iy u 
t= f p-VdV > I=f|4 =) p-\v|-dV 
ve L ve w 


The balance of momentum states that the total time derivative d/dt of the 
momentum is equal to the sum of all external forces. External forces may be mass 
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forces Fm and surface forces F4. 


d á d 4 4 
_ vedV = F F 
dt dt vo" v > m+ >) A > 


dl, 
dt 
Fu Fax 
dl — 
=> Pa =X Puy +), Pay 
Fu: Pay 
dl, 


dt 
In what follows we will consider the time derivative of the integral more closely. Since both 
the region of integration V and the integrand p-v depend on the time, the simplest manner 
of procedure is to form the derivative d/dt as the limiting value of the difference quotients: 


d 
2 p-V-dV = 
1 V(t+At) V(t) 
lim — J pt + A0) st HAN) av - | PHO- Va) -dV 
At—0 At o ( )- v( ) 5 (t) - v(t) 


The first term is dealt with using the additivity of the integral 


V(t+At) 
p(t + At) -V(t+At):-dV = 


0 


V(t) V(t+At) 
J pr An vi AN AV + | pt + At) -V+ At) -dV 
0 V(t) 


Taylor expanding the integrand up to the linear terms yields 


p(t + At) vet At) = p(t) w(t) +A) are. 


We now insert the final two equations into the difference quotients, and find 


d 
2 o-¥-dV = 
1 V(t) d iw V(t+At) 
wo (| ae, plt + At) -F+ At) -dV 
t 


In the next step we reduce the volume integral over the difference V(t + At) — V(t) to an 
integral over the surface A(t) of the volume V(t). Using one-dimensional stream filament 
theory for the mass flux m we obtain the relation 


m=p:c:-A => ZVzeA 


Generalizing this to three-dimensional flows, we compute the volume flux V as a surface 
integral over the scalar product (V -n) of the velocity vector Vv = (u,v,w) with the unit 
vector normal to the external surface n = (Nx, Ny, Nz) 


v= | (wni)-aA 
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The volume flux satisfies 


. _ V(t+At) 
V= gm AADO = in f av= f (V-B) -dA 
At—0 At Ato At V(t) A(t) 


The total time derivative of the momentum yields 


d 
“| p-V:dV = 


=a olp- y / 5 oo ) 
lim [| — -dV + p(t + At)-v(t+ At)-(v-n)-dA 
am (Gf era ME AN) (Fa 





After taking the limit, we find the balance of momentum 





The first term describes the local time change of the momentum inside the control volume 
under consideration. In order to evaluate this integral, knowledge of the flow quantities inside 
the control volume is necessary. For steady flows (0/Ot) = 0. The second term describes the 
convective momentum flux through the surface of the control volume. In order to compute 
this integral, only flow data on the edge of the control volume are needed. 


For steady flows, the balance of momentum reads 


v= | ps dA =Y Fut Y Fa . (2.121) 
A 








With 
Fi U 
By == f pv (ii) aA > IM, =- fp v | (Fi) dA 

A Fiz A w 

we find 

Fı+) Fut) Fa=0| . (2.122) 

Fi Fu. Fax 0 
= Fiy | +X | Fuy |+) | Fay | = | 0 
Iz Fu: Pay 0 


The momentum force vector Fy is parallel to the velocity vector v, and the direction of Fy 
is always towards the inside of the control volume. The pressure force Fp, which is one of 
the surface forces F,, is defined as 


Fp x The 
Bp =— | p-it-da => Foy =- fo Ny - dA. 
A FD. A nz 
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Since the pressure p is a positive scalar quantity, and n is the external unit vector normal 
to the surface, the direction of the pressure force Fp, is, because of the minus sign, also 
towards the inside of the control volume. 


Let us now apply the balance of momentum (2.122) to laminar boundary layer flow. The 
function f which connects the boundary layer thickness with the Reynolds number Re; 
(2.61) can now be determined 


Ô 


mig f (Rer) 


The control volume used in Figure 2.76 is a rectangular block of length /, height ö(L) and 
depth b in the y direction. The pressure p is imprinted onto the boundary layer from the 
outer flow (Op/Oz) = 0, and, in the case of the plate boundary layer, is constant in the outer 
flow. Therefore it follows that the pressure must also be constant in the boundary layer, 
and so all pressure forces appearing must cancel out each other. At position 3, for reasons of 
simplicity, we assume a linear velocity profile u(z), since, compared to the Blasius boundary 
layer profile, this can be integrated analytically. For u(z) we then have 


Uoo 


At the left edge, at position 1, the flow passes through the cross-section A = b - ô&(L) with 
the constant velocity us. The momentum force Fy s, therefore reads 


Fie == | p: uv: (V-B): dA 
Ay 


Üss =l 
nel | ee 
0 0 





Fig. 2.76 : Forces on the control volume V for the laminar plate boundary layer 
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For the momentum force F7,,,, with dA3 = b- dz, it follows that 


Pin | p-u(z)- (vn) -dAs 
Asa 
UZ 1 ö(L) 
-- J p-ulz)- | 0]-{0 ddp | u’ (z) -dz 
As 0 0 0 


Therefore the momentum force F7,,, points in the +x direction and 77 „, in the —:r direction. 
Computing the integral yields 


ö(L) 2 








Eu 2 Uso 2 
Pis =p: f u (z) -dz = —p-b.- P 2". dz 
u 1 ö(L) 
= — = —— b- ur, -d(L 
nea El pew 2U) 


Before we compute the momentum force vector Fi, we will first use the conservation of 
mass. The mass flux m} entering through the surface A; = b- 6(L) is 


The mass flux mg exiting through the surface Ag is 


A3 ö(L) u ö(L) 
= | ulz)-dA3 = p-b- u) | z- dz 
i i ry I, 
Uoo 1 4 1 
— p.b. 2. |L. — =~. pe usb O(L 
mba [pe] spate bo 


Since m3 < m, and the plate is impermeable to the flow, the difference in mass fluxes 
M = Mı — Ma must pass out through the surface A> = b- L. 


fha = fh — tits = p -Uo = b: (L) — $ + p- uoo b: ACE) = Í pu b: 8C) 
Therefore the boundary layer has a displacement action, and the flow through the surface 
A> produces a momentum force Fy,. For the velocity component with which the mass flux 
mg passes vertically out through the surface Ag, we first assume an unknown component 
W(x) > 0 in the +2 direction. Further, ws(x) must satisfy the supplementary condition 
Wz K Ugo. The velocity of the fluid in the x direction along the surface Ag is of magnitude 
Uso. By definition, the momentum force vector Fy, is first of all computed quite generally to 


n=- pv (veri) day =~ | p-{ 0 0 |- foll -d4 , 
Aa A2 wlz) w(x) 1 
Fi. U 
Fy, = Fy yo = -/ p: 0 W(x) - dA» 
Fi z Aa walk) 
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For the x component Ft s, of the momentum force vector Fy, we therefore obtain 
1 9 
Firs = —Uso * p> W2(X) + dAg = Us: Ma => Pub ö(L) 
Asa 


The z component F7,,. has, because of w2(x) a very small magnitude, points in the z direction 
and is of no importance for the following considerations. 


The wall friction force Fw, is the force which decelerates the free stream velocity us to the 
value zero at the surface of the plate. Therefore it points in the negative x direction, and 
because du/dz > 0, we have 

L l 


L L du 
Fb: | Tw dx =—b- f u: — de= bu | —— : dı 
" 0 |e | 0 dz | zo o 9x) 


The balance of momentum in the x direction yields 





| Fra |> | Fres |> | Fres | — | Fwa 0; 


1 1 
P Uso ` b: ACL) = 3 P Uoo t b ACL) = 3 p Ua b: ACL) 
L ] 
0 6(x) 


1 L l p-u L l 
— 0 Us ÖL) = | — de > = .5(L) = | —— -dz 
or men o ôlz) 6- p (2) o O(#) 


If we differentiate the last equation on both sides by s and take the relation v = u/p into 
account, we obtain 





Us dd(x) 1 oy 
= 


dz 


4 
œn 
— 
È 
Qu 
OH 





6-v dg O(a 


Ne 


Integration yields 


ö(L) G. L 1 KL). 


Deve Lk È 12v 2 2 


Uoo I us: L et Re, 





> (L) 





The originally desired function 6/L = f(Rez) is thus found to be 


ö 12 3.464 


L Rer, Rer, 





3 
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The factor 3.464 is a consequence of the simplifying assumption of a linear velocity profile 
u(z). The exact value, found by using the real Blasius profile for the boundary layer, is 5, 
so that 

5.0 


Ô 
L V Re; 





The lift A and the drag W of a wing can also be determined directly with the balance 
of momentum (2.122) if the time-averaged wake profile a(z), p(z) is known. The time- 
averaged velocity and pressure distributions are measured on a fixed control volume V in 
a wind tunnel and numerical integration is used to determine the momentum and pressure 
forces. The balance of momentum in the x direction is written, using the notation of Figure 
2.77, as follows 


Fro — Fin — Fre + Fo. — Fpwn + Fp — W = 0 . (2.123) 


Here the index N denotes the time-averaged profiles in the wake, oo the unperturbed free 
stream, and Fr, Fp the momentum and pressure forces due to the displacement action of the 
wing. The minus sign in front of the drag force is because the drag enters into the balance 
law as a reaction force. The balance of momentum in the z direction yields 


Fis + Fp- 4A=0 . (2.124) 


Application of equations (2.123) and (2.124) is one of the usual methods to determine the 
drag and lift forces of bodies in a flow from the measured velocity and pressure profiles, 
without solving the fundamental fluid mechanical equations directly. The weight G has to 
be taken into account in a separate balance of the mass forces. 


—»5| 


He —— DN 
n «- (WI =N 


(œ) YA | 7 (N) Fig. 2.77 : Forces on the control vol- 
— F 
n l 


ume V for the flow past a foil 
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2.4.3 Balance of Angular Momentum 


For many applications, particularly in the area of flow machinery, a statement analogous to 
the balance of momentum is important for the moments. The balance of angular momen- 
tum can be used to determine the points of action of the momentum forces, or the work 
done or energy gained in the flow through a rotor. 


The angular momentum L is a vector lying in the plane spanned by a displacement vector 
r and the momentum vector I= m.-v. The angular momentum is 


Ly Ty u Ty W— Trz: V 
Ly | = ry | X [vu :m=|r- uU-r,.:w|-m 
L, ly w Tg’ U— Ty u 


The differential angular momentum dL of a mass element dm = p- dV yields 
dL = (£ x V) -dm = p- (É x V) -dV 


Therefore the angular momentum of a volume V(t) is 


pve) 
c= | p-(Ex¥)-dV 
0 


The balance of angular momentum states that the total rate of change d/dt of the 
angular momentum L is equal to the sum of all the external moments acting 
Y M.. 

These external moments >.M, are due to the mass and surface forces Fu + SEF A; 


already discussed in connection with the balance of momentum. Here they act on a lever r 
and we have 


SMa = S (Fx Fu) + S(t Fa) 
The balance of angular momentum reads 


dh d sv® 


_ (¢x¥)-dVv-S°M, . 2.125 
a al, p: (Tx V) > (2.125) 


The formation of the total time derivative is carried out completely analogously to the 
process described in connection with the balance of momentum. We obtain 


dL d , 4 
a a Po exw) dV 


= [COD ay fpe oa EA, 
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Just as in the case of the balance of momentum, for steady flows (0/0t) = 0 the volume 
integral drops away and we only require the surface integral and the flow data on the edge 
of the control region 


fo EXD a) dA = OM, 


A control volume at rest in a steady flow will of course satisfy the steady flow condition as 
long as the coordinate system is at rest too. 


Now a flow machine with a rotor in a coordinate system at rest generates an unsteady flow. 
Here we first carry out a transformation of reference frame to a coordinate system rotating 
with the rotor, in order to generate a steady flow. If we define the moment of momentum 
My. in analogy to the definition of the momentum force, as an inertial moment: 


Mi=- | pE x 0) (P) -dA , 
A 


we obtain the balance of angular momentum 


(2.126) 





The moment of momentum vector My is locally parallel to the vector product (% x vV), 
because the scalar product (V - n) only makes a contribution to the sign and the magnitude 
of the moment of momentum, and not to its direction. 


In order to illustrate the balance of angular momentum we will now consider an example 
of its application. Figure 2.78 shows an elbow bend, flanged onto a pipe. The elbow bend 
steers the flow from the vertical flow direction to the horizontal flow direction. The flow 
exits into free surroundings at the right end of the elbow bend. 


We want to consider how large is the moment Mpx which the elbow bend causes to act on 
the flange junction. Here we assume that the dimension / (Figure 2.78), the flow velocity c, 
the density p of the fluid and the cross-sectional area A, are known. 


Evaluating the integral 
Mi=- | p Ex) (n) AA (2.127) 
A 


for the control surface shown in Figure 2.78, we obtain the scalar magnitude of the vector 
M; as 


|My |=p-1-2- Ay. (2.128) 


In evaluating equation (2.127), the following points should be noted. At position 1, the fluid 
flows over the edge of the control space. The expression beneath the integral in equation 
(2.127) is equal to the zero vector for this section of the control surface, since r x V = 0. 
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control volume V 





Fig. 2.78 : Elbow bend 


At position 2, on the other hand, the cross product delivers a vector in the negative axial 
direction. It points into the plane of the sketch, and therefore its magnitude acquires a 
minus sign. The scalar product V -n is positive at the position 2 and evaluates as c- Aj. 
Taking these details into account, the scalar value for My as formulated in equation (2.128) 
is obtained . 


Other than this, no other resulting forces which produce a moment act on the control surface. 
The elbow bend induces the moment —M, on the fluid. The direction of rotation of Mk is 
initially assumed to be positive. The actual rotation direction is then calculated using the 
equation 


Mı+) M,=0 
The following equation for -Mx =< M, is obtained 
p-l-2-A-|M|=0 > |M|=p-1-2-A 


Therefore the fluid imposes a moment on the elbow bend which acts in the positive direction. 


2.4.4 Pipe Hydraulics 


The objective in this section is to determine the velocity distribution u(r) and, as a supple- 
ment to Chapter 2.3.2, the pressure loss Ap and friction loss cp for laminar and turbulent 
flows through circular pipes. 


We start out from the steady laminar Hagen-Poiseuille pipe flow in Figure 2.43. The 
flow is fully formed, i.e. the velocity profile u(r) depends only on the radial coordinate r 
and does not vary with x,(Ou/0x) = 0. The flow is driven by a constant pressure difference 
in the flow direction x, i.e. (dp/dx) = const. < 0. 


We have already met the parabolic velocity profile u(r) (2.63) which is the analytic solution 
of the Navier-Stokes equation (2.62). Our first step in this chapter on pipe dynamics will be 
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to determine the same result again with the force balance on a cylindrical volume element 
dV = 7-r?-dx as sketched in Figure 2.79. In a fully formed pipe flow, no resulting momentum 
forces act, so that the only forces are pressure forces. 


The pressure force at position 1 (pı > pa) reads 
| Fo. perea = paper 
The pressure force at position 2 is 
> d 
| Fp |= pe n:r” = (v4 Par) nr“ 
x 


The friction is 
Fr |=|r | -2-#-r-de 


Since the velocity distribution u(r) decreases from a maximum value in the center of the 
pipe Umax to the value zero at the wall of the pipe, for r 4 0 we have (du/dr) < 0. Therefore 
the magnitude of the shear stress is 


|7 |= =e: 


The balance of forces follows as 


| Ép, | — | Fp, | — | Fr |= 0 


, 


d 
= per? (pt Bode) m |r| 2emerede=0 
£ 


dp 5 dp r du 1 dp r 
-E m = Dir: zone a 
de“ rl | ee a We) dz 2 7 dr pu dz 2 


This equation corresponds to the first order ordinary differential equation (2.62) to deter- 
mine the desired velocity distribution u(r). After separating the variables and integrating 
indefinitely, we first of all obtain 





Fig. 2.79 : Balance of forces for Hagen-Poiseuille pipe flow 
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The constant of integration C is determined using the boundary condition u(r = R) = 0 as 


l dp i» 
C= ——— -— -R 
4-u da 


The velocity profile u(r) is therefore 


1 dp 2 1 dp 
4-u dx 4-u dz 4-u dz 


(2.129) 





We therefore obtain a parabolic velocity distribution for u(r) with the maximum velocity 


t aD pe 
4-u dz 


Umax 7 


The volume flux in the pipe V is then 


R R r2 
v= fuaa = f u(r) -2emerede= | Una (1-72) 2 werde 
A 0 0 


R r3 
= 2: tinas | (r-i) o 


R 
. 1 1 1 
V = 2 . T . Umax ; . r? — — > 5 = 2 . T . Umar . — > R? 
2 F 4 





Lte g p ME A yA 


The volumetric average of the pipe velocity Um follows as 


1 1 dp œ 
Um = 7° Umax = TD aT 
2 8-y dx 


The volume flux can then be determined as follows: 


(2.130) 





Therefore, for laminar Hagen-Poiseuille pipe flow we obtain the following proportionality 
relation at the position g = L: 


d 
Vr Ap=L-2 | VR 
dx 


Equation (2.130) clarifies the characteristic form of the volume flux. It is proportional to 
the pressure loss Ap = p; — po and to the fourth power of the radius R. 
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We now consider the size of the pressure loss Ap for a given volume flux. This pressure loss 
is a consequence of the effect of friction. From (2.130) 











gn LD u 
V=—.-.— -R Ap = 9, — 
it follows that 
DW m- Ri az z-R mu i R? = ~~ Reo 


In what follows we expand the term on the right-hand side of Ap so that characteristic 
quantities of the flow can be summarized: 
1 s oy d = Iover. 1 


Nuss ie ee ge 
pas Piar p = Pw pou 





Sl 


2 
m D\2 m ` üm: D 
Um: (F) v 

Defining the Reynolds number formed with the diameter of the pipe D as Rep = (um: D)/v 
and writing the factor 64/Rep as a loss coefficient Aşım, we obtain the following equations 


to determine the pressure loss 


(2.131) 





These equations hold for laminar pipe flows, i.e. for Reynolds numbers smaller than the 
critical Reynolds number Ree 


dD 
Re — < Re, = 2300 





For fully formed turbulent pipe flow the time-averaged velocity satisfies (OU/Ox) = 0, so 
that again momentum forces occur in the time average. In Figure 2.80 we apply the balance 
of momentum to a control volume V = 7- R?- L with pipe radius R, and determine the 
pressure force at position 1 (pı > p2) as 


| Fp. [=p r R? 





Fig. 2.80 : Balance of forces on a control volume V for turbulent pipe flow 
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The pressure force at position 2 is 
| Foe |= p-r- R? 
The wall friction force is 


| Fru |=| 7a | 2 m- R-L 


? 


and the balance of forces follows as 


| Ep1 | — | Ep2 |- | Fru |=0 , 
Por R- por R- |I| 2r R L=0, 


2- L 
(p-p) T: R = Ap r: RP =| a| 2r RL = Ap=| Tu | = 


There is no theoretical ansatz for the wall shear stress | 7,, |, and so we use an empirical 
ansatz, which determines the pressure loss equation Ap in analogy to the laminar case: 


Ww =. . .— => A zZ — >. . —_.—— — — . . —_. 


À 
™ 4 R 2 to? 


A, = (Rep) from experiment, Rep = 





From experimental results, we obtain the pressure loss coefficient A; in the form of the 
Blasius law 


0.3164 
A = - , valid for 3-10°< Rep < 10° , (2.133) 


(Rep)? 





and the implicit Prandtl representation 


1 
— = 2. logio (Rep . VA) -08 , valid for Rep < 10° . (2.134) 
VAt 


If the pipe is rough, values of A; can be read off from the Nikuradse diagram in Figure 
2.81. The roughness K, here is the spatial average of the surface roughness of the pipe walls. 
Some values for different materials are listed in Figure 2.82. 


The extension to the implicit equation (2.134) for rough pipes, determined from experiment, 


yields 
1 2.51 K, 
= =2- Jog, ( 22 | 2.135 

A 210 N 37. 5) (2.135) 


For Reynolds numbers Rep > 10° the loss coefficient >, becomes independent of the 
Reynolds number, since then the viscous sublayer of the turbulent pipe boundary layer 
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100A, D 


12 
10 


2,9 





10° 104 10° Re, 10° 
Fig. 2.81 : Nikuradse diagram 


covers the roughness of the pipe. Figure 2.83 shows ., for different pipe cross-sections. 


In order to compute the time-averaged turbulent velocity profile U(r) we start out from the 
ansatz for the wall shear stress F, 
_2 At 


_ 1 

| Tw |= z` 
Using the Blasius equation (2.133) 

0.3164 0.3164 


(Rep)! (Zu)? 


V 





glass, copper, brass 
EE fiber cement 


wood 


cast iron 





f Fig. 2.82 : Roughness of different materials 
10° 10% 10° 1 10 K, 10° 10° 
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and taking the proportionality R ~ D and Um ~ Umax into account, we obtain the relation 
| Tw I~ p’ Tnax " (man) 7 " Ri ° vi =p: (Tmax)? R. pī 
Initially we restrict ourselves to determining u(r) close to the wall for r — R and introduce 


the substitution z = R—r. A power ansatz with unknown exponent m can then be written 
down for the velocity profile u(r) close to the wall as follows: 


U(r) = Umar (5) 


R” 


> üm =) (Whar)? = WA) RE 


Therefore the wall shear stress follows as 


7m 7m 


Tu Im pa (2) RT 


_ Im 


1 1 
4.7 4 . pa 
Prandtl and von Karman introduced the hypothesis that | 7, | should be independent of 
the pipe radius Æ for a turbulent pipe flow, i.e. the exponent of R should vanish 

s Om fig s 

— — — — mM = — 

4 4 T 
Substituting r back into the expression for the time-averaged velocity profile, we obtain the 
(1/7) power law for turbulent pipe flow 





(2.136) 


For m = (1/7) the average velocity üm is 
Um = 0.816 - Umaz 


The region of validity of this law is the same as that for the Blasius equation (2.133), 
Rep < 10°. 





Fig. 2.83 : Pressure loss A, at different pipe 
H/B(D-d)/D cross-sections 
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Let us mention two unphysical disadvantages of this profile. At the pipe wall we obtain an 


infinitely steep increase in the velocity 
du 
dr r=R 


However this is unimportant, since the law is not valid in the viscous sublayer. 


At the center of the pipe there is a bend in the velocity profile and (dü/dr)(r = 0) is not 
defined. 


The parabolic velocity profile of laminar pipe flows (2.129), as well as the time-averaged 
velocity profile of turbulent pipe flows (2.136) are compared in Figure 2.84 for equal volume 
fluxes V. 


All that remains to do is determine the thickness of the viscous sublayer A. Using the ansatz 


ioe a At du 
me, m g7” dz) y 


we obtain a linear increase in the velocity from the value zero at the wall to the value 0.5 -üm 
at z = A, therefore 


diu\ _ 5° Um u du) _, a 2 & 
dz A a eu o ea 


w 








The thickness A of the viscous sublayer therefore follows as 
4. v A 4 V 4 


= => me gegen. 
üm Àt D A; üm: D Rep Ar 
Taking the Blasius law (2.133) into account 
0.3164 


t= 1 


(Ren)* 








it follows that 





A 
= , 2.1 
D (2.137) 


turbulent 


Fig. 2.84 : Velocity profiles for laminar and 
turbulent pipe flows 





126 2. Fundamentals of Fluid Mechanics 


2.4.5 Flows Past Bodies 


In problems involving flows past bodies, one of the most important phenomena is that of 
flow separation, which we met several times in previous sections. Depending on the size 
of the Reynolds number, flow separation can take place in a steady or unsteady manner. 


Let us consider the flow past a sphere. The left-hand picture in Figure 2.85 shows laminar, 
steady flow separation at low Reynolds numbers. The separation of the boundary layer 
on a sphere leads to a region of backflow. Because of the acceleration downstream of the 
stagnation point, the pressure p on the surface of the sphere decreases greatly and takes 
on a constant value in the backflow region. Separation of a turbulent boundary layer takes 
place at larger Reynolds numbers much further downstream on the surface of the sphere. 
Because of the flow deceleration, the pressure p first increases again beyond the widest point 
of the sphere, to then take on a constant value in the turbulent backflow region. 


The flow separation on a sphere can be made plausible with the following considerations. The 
separation of a flow from a wall occurs when the no-slip condition close to the wall causes 
the boundary layer fluid to decelerate and to be transported into the main flow region. If the 
pressure in the outer flow increases downstream, the slower fluid within the boundary layer 


laminar turbulent 


0° 90° ® 180° 0° 90° = 180° 


Fig. 2.85 : Separation and pressure distribution in the flow past a sphere 
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is no longer able, because of its lower kinetic energy, to flow downstream into the higher 
pressure region. Because the turbulent boundary layer flow has an additional longitudinal 
and transverse momentum exchange and so a higher kinetic energy, the turbulent boundary 
layer can remain attached to the surface of the sphere further downstream. The backflow 
region tapers off, as does the wake flow, so that the total drag Cw is considerably reduced. 


In order to achieve a mathematical description of the separation criterion, we look at a 
two-dimensional laminar or turbulent boundary layer on, for example, a cylinder. Because 
of the no-slip condition at the wall vu = 0 and w = O for r = R, with the radius of the 
cylinder R or z = 0, it follows from the Navier-Stokes equation in Cartesian coordinates 
(2.65) that 


1 dp u 


= UV: — 2.138 
p da 2 ) 





z=0 


Equation (2.138) and Figure 2.86 can be used to discuss the development of the boundary 
layer flow with the pressure gradient. If the pressure in the x direction decreases, i.e. if 
Op/Ox is negative, the flow outside the boundary layer is decelerated downstream. In this 
case (07u/027) < 0 too, and the curvature of the velocity profile u(z) at the wall is negative. 
Because of the acceleration of the flow, the velocity at the edge of the boundary layer 
increases, with the result that Ou/Oz increases with increasing z downstream. Because of 
Tu = u (Ou/öz)z = 0, the wall shear stress 7, then also increases with increasing x, so that 
(OF / On) > 0. 


If (Op/Ox) = 0, equation (2.138) implies that 02u/02* will also be zero at the wall, i.e. 
the velocity profile u(z) has a turning point at the wall. The velocity at the edge of the 
boundary layer remains constant because there is no pressure gradient present. Inside the 
boundary layer however, the flow is decelerated by the friction forces present. This means 
that close to wall, the velocity gradient Ou/Oz decreases with increasing x. This leads to a 
reduction of the wall shear stress Tẹ in the x direction, with (O7,/0x) < 0. 


The flow separation from the contour of the body begins at the position where the upstream 
positive wall shear stress Tą sinks so low that it assumes the value zero for the first time. 





Fig. 2.86 : Boundary-layer separation 
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This is the criterion for the start of separation. 
separation criterion: Ty =O} . (2.139) 


Figure 2.86 shows the boundary layer separation for the case of a positive pressure gradient 
(Op/Ox) > 0. A positive pressure gradient leads first of all to the flow outside the boundary 
layer being decelerated in the x direction. In the figure this is indicated in the way the 
velocity arrows at the edge of the boundary layer become ever shorter with increasing x. 


Because (Op/Ox) > 0, equation (2.138) implies that the curvature of the velocity profile 
at the wall must satisfy (0?u/0z7) > 0. At greater distances from the wall, the curvature 
of the velocity profile u(z) is always negative. Therefore, if the curvature at the wall is 
positive, with (0?u/0z?) > 0, there must be at least one point inside the boundary layer 
where (6°u/0z7) = 0. This point is a turning point of the velocity profile u(z). 


In contrast to the start of separation where the turning point is at the wall, the turning 
point moves inside the boundary layer downstream of the start of separation. Figure 2.85 
demonstrates the consequences of a positive pressure gradient (Op/Ox) > 0. In this case the 
boundary layer flow is not only decelerated by friction forces, but also by pressure forces, 
and the curvature at the wall is always positive. ‘The wall shear stress 7,, decreases in the x 
direction, at separation begins at Ty = O. 


In the two-dimensional case this is equivalent to (Ou/Oz) = 0. Further downstream, the 
wall shear stress becomes negative. This implies a reversal of the flow direction close to the 
wall with (Ou/O0z) < 0 and hence backflow. The backflow leads to a recirculation region 
downstream of the separation point. 


Now that we have worked out the separation criterion for two-dimensional boundary layer 
flow, let us return to the case of the flow past a sphere and discuss the Reynolds number 
dependence of the drag coefficient Cw = cy(Rep), where D is the diameter of the sphere. 


We start our discussion of the dependence of ¢, on the Reynolds number for the case of 
Reynolds numbers Rep < 1. In this case the friction forces are much greater than the inertial 
forces, and we speak of creeping flows, which can be described analytically. For the case of a 
sphere in a steady flow with Rep < 1, the analytic solution of the Navier-Stokes equation 
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leads to a drag force W of 
D 
WE ON Us i (2.140) 


One third of this drag force W is due to the pressure gradient, and two thirds to the friction 
forces. We also note that for creeping flows, the drag force W is proportional to the first 
power of the velocity u... Taking the definition of the c, value into account, we can used 
equation (2.38) to write down a relation for Cw = «w(Rep): 


W 24- u 24 


I pu, 2D? pus: D Rep 


(2.141) 


Cy = 


The relation c,, = (24/Rep) is also known as Stokes drag law and is valid in the Reynolds 
number region Rep < 1. If the Reynolds number is increased to a value of Rep = 300, 
steady separation takes places downstream. The fluid particles directly at the wall lose so 
much kinetic energy through the strong friction forces that they are unable to compensate 
the pressure increase in the rear half of the sphere. This leads to flow separation close to 
the equator of the sphere. A steady backflow region is found in the wake directly behind 
the sphere. In computing the steady wake flow the inertial terms can no longer be ignored, 
and the Navier-Stokes equations have to be solved. 
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Fig. 2.88 : Drag coefficient ¢, and recip- 
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A further increase in the Reynolds number to a value of Rep = 2000 initially leads to the 
formation of unsteady vortex separation of the laminar boundary layer on the surface of the 
sphere with a laminar wake. 


Above Rep = 2000, to about Rep = 3- 10°, the transition to a turbulent wake flow takes 
place. In the Reynolds number region 3- 10° < Rep < 4-10°, the boundary layer flow on 
the sphere becomes turbulent. The separation region moves downstream on the surface of 
the sphere and the wake flow tapers off. This leads to a drastic reduction in the ¢, value, 
as shown in Figure 2.87. The friction drag is larger in a turbulent boundary layer, i.e. the 
drop in the ¢,, value is due to a reduction in the pressure drag. 


In the region 4- 10° < Rep < 10° the laminar-turbulent transition region moves forwards 


PRANDTL 1927 
Start-up of a cylinder from rest, intermedi- 


ate states in the formation of the Karman 
vortex street 


HOMANN 1936 


Constant flow past a cylinder 


Fig. 2.89 : Steady flow past a cylinder and 
laminar Karman vortex street 
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on the surface of the sphere, increasing the friction drag, while the pressure drag essentially 
remains constant. Again the Cw value increases. In the Reynolds number region Rep > 
10°, the boundary layer on the surface of the sphere is turbulent downstream of the front 
stagnation point, fixing the separation point. This does not change for any further increase 
in Rep, and so the c, value of the sphere also becomes independent of Rep. 
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Fig. 2.91 : Flow separation at a wing profile 


In the case of the flow past a cylinder the dependence of the drag coefficient Cw on the 
Reynolds number Rep is equivalent. Figure 2.88 shows all known experimental values cy 
with the measured reciprocal values of the dimensionless separation frequency 1/Str. These 
are the values for the Karman vortex street already described in Chapter 1.1. Figure 2.89 
shows the flow past a cylinder for the steady flow separation region where 3 < Rep < 40 
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and for the laminar Karman vortex street region where 40 < Rep < 200. 


The periodic vortex separation of the Karman vortex street starts at the Reynolds number 
Rep = 40. With increasing Reynolds number, 1/Str decreases greatly, taking on an almost 
constant value of Sir = 0.21 for Reynolds numbers between 10° and 10+. After the transition 
to turbulent boundary layer flows on the cylinder, 1/Str drops greatly, corresponding to the 
drop in the drag coefficient cw. For Reynolds numbers greater that 10’ a constant separation 
frequency is found in the turbulent wake flow with constant Cu value. 


The drag coefficient c,, in the Reynolds number region from 10* to 10° are shown in Figure 
2.90 for different three-dimensional and two-dimensional body shapes, with the Reynolds 
number dependence for different axisymmetric bodies also shown. 


Flow separation also takes place in the flow past a wing upstream of the trailing edge 
of the wing and in the shock wave region. Let us return to Figure 1.18 of the introductory 
Chapter 1.2 and discuss the effect of the angle of attack a of the wing. The left-hand side 
of Figure 2.91 shows a subsonic wing profile at four different angles of attack a for a free 
stream with Mach number M. The regions where the friction effect is large are shaded grey. 
Case (a) shows the wing profile with a = 0°. The diagram on the right is the lift coefficient 
Ca aS a function of the angle of attack a. It can be seen that for a = 0° the lift coefficient 
satisfies c, > 0. As we have already discussed, this is a consequence of the profiling of the 
wing, which causes the fluid on the upper side of the wing to flow faster than on the lower 
side and thus generates a pressure difference. The friction effect mainly dominates in the 
boundary layer and in the wake of the wing. Directly behind the trailing edge of the wing a 
shear layer of the boundary layer forms above and below the wing, passing over to the wake 
flow downstream. 


A slight increase in the angle of attack a, as in case (b), leads to a considerable increase 
in the lift coefficient c,, without there being much effect on the wake, compared to case 
(a). In case (c), a further increase in the angle of attack leads to a maximum lift coefficient 
Ca for a certain angle a. However, here separation already occurs on the upper side of the 
profile, shown here as an extended shaded friction-dominated region downstream. A further 
increase of the angle of attack a leads to an expansion of the separation region downstream 
on the upper side of the wing, and thus to an increase of the lift coefficient. In the limiting 
case, separation will take place at the leading edge of the wing, and there is no lift at all. 


In the case of the transonic wing, the supersonic region on the wing, as already seen in 
Figure 1.18 and Figure 2.58, is ended by a shock wave. This shock wave causes a pressure 
jump in the wing boundary layer, which again can lead to flow separation if the boundary 
layer fluid cannot overcome the pressure increase. 


The conditions under which shock-boundary-layer interaction will lead to flow separation are 
shown in Figure 2.92. Velocity profiles in the boundary layer u(x, z), the sonic line M = 1, 
the pressure distribution p and the wall shear stress 7,, are shown. In the supersonic region 
in front of the shock wave, a supersonic flow already occurs inside the turbulent boundary 
layer. Below the dashed sonic line is a region where perturbations may spread out upstream. 
It is also shown how the shock wave fans out into different weaker compression regions due 
to the effect of the boundary layer. Behind the shock, the pressure in the turbulent boundary 
layer becomes so large that the wall shear stress 7, tends to zero. At this point separation 
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on the wing profile begins. 


Further downstream both the wall shear stress and and velocity take on negative values. 
Close to the wall backflow occurs, as shown in the middle velocity profile of the three shown 
in Figure 2.92. In the separation and backflow regions, the pressure is considerably greater 
than that without separation. This has an effect on the lift and the drag. Whereas the drag 
increases, the lift decreases. Because of the further acceleration of the flow on the profile and 
the dying away of the pressure increase in the turbulent boundary layer due to the shock 
wave, the flow reattaches again after as certain distance. At this point the wall shear stress 
again becomes zero, and takes on positive values further downstream. Therefore a bounded 
separation region forms on the wing. Again the expansion of the separation region is of 
interest, since this affects the global flow field, i.e. the outer flow. Similarity considerations 
for the length l of the separation region in the flow direction lead to the result that larger 
Reynolds numbers lower the length /, while higher Mach numbers increase J. 


Having now looked at the flow separation on a wing, we now return to the same phenomenon 
in the incompressible flow past an automobile. Whereas flow separation on the wing must 
be avoided in order to retain the lift, it is an essential component in the reduction of drag 
in the flow past an automobile. We have already got a first impression of the separation 
regimes in the flow past a vehicle in the introductory Chapter 1.2 in Figure 1.20. 


For the case of an automobile with notchback, let us compute the flow separation on the rear 
window and on the edge of the trunk cover. We are now immediately able to understand the 
flow behavior with one look at Figure 1.21. This shows the qualitative pressure distribution 
on the contour of a notchback vehicle. Positive pressure differences are seen in the regions in 
danger of separation mentioned above. Those regions, denoted with an © sign, are simply 
the regions where a positive pressure gradient Op/Ox is present, which, as we now know, 
leads to flow separation and backflow. 


The separated boundary layer on the back of the automobile in Figure 1.22 generates a 
part of the wake flow of the vehicle as a free shear layer after passing the edge of the trunk 
cover. In these region there is a point of departure to influence the flow deliberately in 
order to reduce the drag. The aim is to remove the shear layer behind the spoiler resulting 
from separation, and so to restrict the formation of a certain region of the wake flow which 
contributes greatly to the drag on the vehicle. 


Flow separation can also occur in pipes and diffusors. Starting out with the previous Chapter 
2.4.4, we consider the flow separation in the supply pipes. The flow separation here also 
causes additional losses. We consider the elbow bend in Figure 2.93, which steers a vertical 
flow to a horizontal flow. We assume a steady fully-formed pipe flow in the straight vertical 
piece of the pipe, in which a driving pressure gradient is present in the flow direction, and 
constant pressure in the radial direction transverse to the flow. 


The Bernoulli equation for curved stream filaments leads us to the statement that the 
pressure increases in the radial direction in order to keep the centrifugal force in equilibrium. 
A pressure gradient is formed transverse to the flow direction, leading to a pressure increase 
at the outer wall and a pressure decrease at the inner wall of the elbow bend. This works 
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against the pressure drop along the streamline coordinate s of the outer wall and boosts it 
at the inner wall. The streamline coordinate s denotes the arc length of a stream filament 
and is read positive downstream. 


In the previous examples we have already seen many times that a pressure increase in the 
flow direction leads to flow separation. Therefore separation will first occur at the outer 
wall at position A. At the exit from the elbow bend, the pressure again levels out transverse 
to the flow direction. The pressure at the inner wall increases and that at the outer wall 
drops again. This leads to reattachment of the flow A, at the outer wall, and to the start 
of flow separation at the point B at the inner wall. At the inner wall, the flow will again 
reattach with increasing arc length s some distance after passing the elbow bend in the 
straight horizontal pipe section B,,. Again a negative pressure gradient Op/Os is present 
here, keeping the friction forces in equilibrium. The pressure transverse to the flow direction 
is again constant in this straight section of the pipe. 


It can be seen in Figure 2.93 that, downstream of the separation points A and B at both 
the outer and inner walls, recirculation regions form, which cause additional energy losses in 
the flow. The lower picture in Figure 2.93 shows the pressure in the pipe for two streamlines 
in the outer and inner wall regions along the streamline coordinate s. The inclined straight 
line shows the linear pressure drop in a straight pipe section. Even if there is no separation, 
the energy losses in the flow due to the friction yield a pressure loss in the flow direction. 


Above the straight line, the thick curve shows the pressure of a streamline at the outer wall 
which would be found without separation. Below the straight line is the corresponding curve 
for a streamline at the inner wall. The separation at both point A and point B takes place 
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B ee Fig. 2.93 : Flow separation in an elbow 
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in a region of increasing pressure. The additional flow loss caused by separation is seen in 
the way that the dashed pressure curves at the outer and inner walls of the elbow bend lie 
below the curves without separation. 


As well as the flow separation, a secondary flow also occurs in the elbow bend. This is 
superimposed onto the main flow in the direction of the streamline coordinate s and causes 
considerable velocity components perpendicular to the main flow. The origin of this sec- 
ondary flow is the curvature of the pipe, as well as the deceleration of the flow from friction 
forces at the wall. The velocity at the inner side of the elbow bend is greater than that at 
the outer side. The fluid close to the wall has, because of the friction, a lower velocity than 
that in the center of the bend. The centrifugal forces which are larger in the center of the 
bend than at the side walls cause an outward motion. However, for reasons of continuity, 
this is only possible if a motion in the opposite direction takes place at the walls of the bend. 
Therefore a double vortex forms and is superimposed onto the main flow. The secondary 
vortices also lead to flow losses, so that we can divide the losses in an elbow bend into the 
following three components: friction losses and separation losses due to the curvature, 
and losses due to secondary flows. 
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3 Software and Applications 


It is not only in fluid mechanics that the field of activity of engineers has seen considerable 
changes due to the increased use of computers and networks. As well as having the analytic 
abilities needed to solve fluid mechanical problems, the engineer is increasingly required to 
be able to handle fluid mechanical software in practical applications in industry. In order 
to promote this development, we have supplied tutorial software in conjunction with this 
textbook, to ease the introduction to the use of commercial fluid mechanical software. It is 
essential for future engineers to be able to work with fluid mechanical software comfortably 
and independently. 


Fluid mechanical software accompanying the lectures is available on the internet, via the 
homepage of the Institute for Fluid Mechanics at the University of Karlsruhe. Interaction 
between students and teaching assistants is encouraged via email and by means of internet 
consultation. 


http: //www-isl.mach.uni-karlsruhe.de 


3.1 Fluid Mechanical Software 


The software accompanying the textbook concerning the fundamentals of fluid mechan- 
ics discussed in Chapter 2 is available in the software module 


KAPPA-Stromfaden 


http://www-isl.mach.uni-karlsruhe.de/stromfaden.html 





Here the algebraic equations of the one-dimensional stream filament theory presented in 
Chapters 2.3.2 and 2.3.3, as well as the two-dimensional Navier-Stokes equation for viscous 
flow, are solved iteratively for given examples. Examples of application are the flow past an 
automobile and a wing, as well as the unsteady flow out of a container (incompressible) and 
the flow through a nozzle (compressible). 


After calling up the page http://www-isl.mach.uni-karlsruhe.de/stromfaden.html or 
else clicking on KAPPA-Stromfaden, the menu shown on the next page appears: 
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KAPPA-Stream-Filament 





© incompressible © compressible 
© inviscid © viscous 


© steady © unsteady 


automobile: © upper side + lower side 


wing: © upper side © lower side 
nozzle: © nozzle flow 
exhaust: © 


free stream velocity in km/h (automobile), 
back pressure / reservoir pressure ratio (nozzle: | 


| start computation | delete input | 

















By clicking on the options offered, the steady, compressible, inviscid flow through a Laval 
nozzle, for example, can be computed. The results are obtained in the form of pressure 
and Mach number distributions, p(x) and M(x), along the axis of the nozzle for a selected 
pressure ratio of the pressure p, at the exit of the nozzle, to the reservoir pressure pp in 
the container to which the nozzle is connected. Therefore KAPPA-Stromfaden permits us 
to study the effect of the pressure ratio on the characteristic flow form in the nozzle. 


For example, it is seen that for a pressure ratio of (p4/po) = 0.98 a pure supersonic flow 
exists everywhere in the nozzle, with a maximum Mach number at the narrowest cross- 
section of Mmar © 0.37. If the pressure ratio is reduced to (p,/po) = 0.9, for example, 
a normal shock wave occurs downstream of the narrowest cross-section. This is seen in a 
discontinuous drop in the Mach number from M > 1to M < 1. For asmall pressure ratio of, 
for example, (p4/po) = 0.1, a continuously accelerated flow in the Laval nozzle is ultimately 
obtained. In this case the Mach number along the axis of the nozzle increases by a factor of 
10 from an initial value of M & 0.22 to about M #& 2.2. 


The KAPPA-Stromfaden software can be used to carry out the preliminary design in 
product development, as introduced in Chapter 1.3. Supplementary fluid mechanical soft- 
ware is available for the calculation of the preliminary design, whereby the mathematical 
basics necessary for this are not treated in this textbook. To this end we recommend the 
more extensive texts: H.OERTEL, M.BÖHLE, Strömungsmechanik, 1999 and H.OERTEL, 
E.LAURIEN, Numerische Strömungsmechanik, 2002. The university-developed software 
package KAPPA (Karlsruhe Parallel Program for Aerodynamics) and the commercial 
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software STAR-CD can be used for calculation, and both are available through the Insti- 
tute of Fluid Mechanics. Within the framework of product development, these packages can 
be used to compute laminar and turbulent, incompressible and compressible, steady and 
unsteady, three-dimensional flows for given product geometries. 


After calling up http://www-isl.mach.uni-karlsruhe.de/skappa.html, and download- 
ing, the following menu appears: 
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In order to use the software package, more than just the textbook is needed. The internet 
instructions given by the teaching assistants can also be studied independently with the 
training course in the following chapter. 


Downloading the page http://www-isl.mach.uni-karlsruhe.de/star-cd.html causes 
the following menu from the commercial software package STAR-CD to appear. A license is 
needed to use this software, and this license also includes training in the use of the software 
package. 
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The fluid mechanical software packages mentioned above can be used to compute the flow 
examples described in Chapters 3.2 and 3.3. 
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3.2 Training Course 


Now that we have already got to know some examples of stream filament theory software, we 
continue with the introduction to the software for the solution of fluid mechanical problems 
in industrial projects. Accompanying the introductory course, taught at the University of 
Karlsruhe (TH) as a practical training course in numerical fluid mechanics, we recommend 
our more extensive textbook H.OERTEL, E.LAURIEN 2002. The introductory course is 
structured so that it can also be studied via the internet. At the end of the introductory 
course, the student will be able to use the KAPPA software package independently to 
compute the flow past a transonic wing. 


The introductory course is to be found at the address 


http: //www-isl.mach.uni-karlsruhe.de/seinfuehrung. html 





Now that the first numerical computation with the software package KAPPA has been 
successfully carried out under the guidance of the teaching assistants, the training course 
for higher-semester students follows. This is an introduction to the application to industrial 
projects. In the practical application of fluid mechanical software, the first thing learned 
is that the computational grid and the model used (e.g. turbulence model), as well as the 
boundary conditions for each class of geometries, have to be set up from scratch, or else 
newly adapted. This first phase of the adaptation of the software is called the verification. 


steady flows past bodies steady flows through channels 


laminar plate boundary layer laminar pipe flow 


turbulent plate boundary layer turbulent pipe flow 
(turbulence models) (turbulence models) 


transonic foil (RAE 2822) backward-facing step 
(turbulence models) 


transonic wing (ONERA M6) 


unsteady flows through channels 
unsteady flows past bodies 


shock tube 
Rayleigh-Stokes problem 


laminar Karman vortex street 


turbulent Karman vortex street 
(turbulence models) 





flows past bodies flows through channels 


Fig. 3.1 : Test verification cases for steady and unsteady flows 
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It requires highly developed engineering skills and is generally only possible under expert 
guidance. In what follows we present some helpful points in the training course, which is to 
be found at the address 


http: //www-isl.mach.uni-karlsruhe.de/sverifikation. html 





and can be used independently, as long as software licenses are available on the computer 
used. 


Before we start computing an actual industrial project, let us run the software package 
through the test cases in Figure 3.1. Here it first has to be determined whether the problem 
concerns a flow past a body or a flow through channels. There are then further subdivisions 
into steady and unsteady test cases. 


Steady Flows Past Bodies 


Laminar Plate Boundary Layer 

The laminar plate boundary layer is a frequently used test case for spatial discretization in 
fluid mechanical software. In this case, a plate of length L at room temperature T,, = 293K 
is placed in a flow with Reynolds number Re, = 1105. 


We compare the numerical results with the analytical Blasius solution: 
0.664 


/ Re(x) 


Re(x) = p+ Uc + &/p is the local Reynolds number formed with the length x. Figure 3.2 
shows the dimensionless wall shear stress 
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For all discretization schemes tested, the numerical results agree very well with the analyt- 
ical Blasius boundary layer solution. 
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Fig. 3.2 : Local friction coefficient of the 
0 0.2 0.4 0.6 yi laminar plate boundary layer, Re, = 1105 
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Fig. 3.3 : Turbulent plate boundary layer, Rez, = 2 - 10° 


Turbulent Plate Boundary Layer 

The two-dimensional, turbulent plate boundary layer without pressure gradient is an impor- 
tant example of adapting turbulence models. A fine enough computational grid can be used 
to keep the effect of the discretization error on the numerical solution very small. To com- 
pute the numerical results an experiment is used. Here the Reynolds number is Rez = 2-106. 
The plate with the length L is in a flow with Mach number M,, = 0.1 and the temperature 
of the free flow is T,, = 293k. The intensity of turbulence in the free stream is given as 
Tus = 0.5%. The transition is fixed with a Reynolds number formed with the momentum 
thickness of Res = 770. The mean velocity profile at the position #/L and the total friction 
coefficient along the plate are compared with experimental data for these parameters. 


Figure 3.3 shows the computed and measured time-averaged velocity profile friction coeff- 
cients. It can be seen that the given turbulence model computes the boundary layer to be 
too full. The effect of the friction coefficients however is small, and there is good agreement 
of the computed total friction coefficients with the experimentally determined values. 


Transonic Foil 

The transonic foil RAE2822 is a test case for a compressible flow past a body. The foil is in 
a flow with Mach number M = 0.73, Reynolds number Re; = 6.5- 10° and at an angle of 
attack a@ = 3.19°, and a temperature of T,, = 300K. The transition is fixed at x/L = 0.05. 


Fig. 3.4 : Pressure coefficient of flow past 
a foil (RAE 2822), Re, = 6.5- 10°, Mx. = 
0.73 
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The intensity of turbulence in the free stream is Tu. = 0.3%. The turbulence length scale 
is estimated via the molecular viscosity and the Prandtl mixing length ansatz as 


j . 
laure © —— with ky = Tus: Us 
PV Ko 
The distribution of the pressure coefficient c, over the length x/L, as found in the numerical 
computation, is compared to the experimental data in Figure 3.4. The pressure distribu- 
tion and the computed shock wave position agree. The lift coefficient is also found to be in 
agreement with experiment with c, = 0.795, and the drag coefficient is c, = 1.7 - 1077. 


Transonic Wing 

The steady, three-dimensional transonic flow past a wing can be used to check the quality 
of the software package to compute compressible, three-dimensional flows past bodies with 
shock waves. The test wing used is the ONERA M6 wing. 
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Fig. 3.5 : Pressure coefficient, flow past a wing (ONERA M6), Rer = 1.17 -10°, Mœ = 0.84 
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A double shock wave occurs on the suction side when the wing is placed in a flow with 
Mo = 0.84. This moves towards the point of the wing and forms one shock wave. This shock 
wave branching makes this wing a particularly critical test case. The Reynolds number of the 
unperturbed free stream is Re; = 1.17: 107 and the temperature of the flow is T,, = 300K. 
The wing is placed at an angle a = 3.06° to the flow. The flow is assumed to be adiabatic and 
turbulent for the computation, i.e. the transition is fixed in the free stream. The turbulent 
kinetic energy is given by Tu, = 0.3%. The numerical solution is evaluated by comparison 
with experimental data in various cross-sections of the wing along its span. 


Figure 3.5 shows the double shock on the wave as the closely-spaced isobars on the surface. 
These two shock waves join together towards the point of the wing. The computed pres- 
sure distributions along the span s of the wing agree well with the experimental data. The 
positions of the shock waves on the wing and the branching of the shock, as well as the 
suction peak along the leading edge, are all in agreement with measurements, in as much as 
a shock-adapted fine grid can be chosen in the shock wave regime. The computed lift and 
drag coefficients are c, = 0.138 and ¢,, = 7.98 - 107? respectively. 


Rayleigh-Stokes Problem 

In order to verify the time precision of unsteady laminar flows, we will use the example 
of the first Stokes problem for a flat plate. Consider a flat plate at rest in the (x, y) 
plane, above which, in the z direction normal to the wall, is a viscous fluid at rest. At the 
time ¢ = 0 the plate is set instantaneously into motion at the constant velocity Us. The 
effect of friction causes the fluid above the plate to be set into motion, and this fluid forms 
an unsteady boundary layer. The flow is now computed for the Mach number M, = 0.1. 
The velocity profiles are similar at different times. Figure 3.6 shows the computed velocity 
profiles drawn against the dimensionless coordinate 7. This is the same as the results of 
the analytical solution for the Rayleigh-Stokes problem, and so is a proof of both the time 
precision of the numerical solution and the spatial discretization. 


Laminar Karman Vortex Street 


A further indication of the time precision of the software is the example of the Karman 
vortex street found at the plane flow past a cylinder at super-critical Reynolds numbers. 
The laminar Karman vortex street is computed at the Reynolds number Rep = 500. The 
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test quantity chosen is the dimensionless frequency of separation, the Strouhal number 
Str = fD/ugo. The computed Strouhal number is Str = 0.209, compared to that found 
by experiment Str = 0.21. There is also proof of the time precision of the software in this 
unsteady test case. Figure 3.7 shows the (computed) snapshot of the pressure distribution 
in the wake of the plane cylinder. 


Turbulent Karman Vortex Street 

Another test case is that of the turbulent Karman vortex street, with which the turbulence 
model for unsteady flow can be chosen. The numerical calculation is carried out using a 
Reynolds number of Rep = 1.4-10* (Fig. 3.8). The computed Strouhal number of Str = 
0.206 is slightly above the experimental value of Str = 0.2. Here again it is seen that the 
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Fig. 3.7 : Snapshot of the velocity distribution of the laminar Karman vortex street, Rep = 
500 
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Fig. 3.8 : Snapshot of the velocity distribution of the turbulent Karman vortex street, 
Rep = 1.4- 10+ 
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turbulence models available in the software are adequate. 


Flows Through Channels 


Laminar Pipe Flow 

Here we consider the laminar flow in a pipe with length L/D = 50 and diameter D as a test 
case. The pipe is hydraulically smooth. For the Reynolds number, formed with the mean 
velocity and the diameter of the pipe, we choose a value of Rep = 660, while the temperature 
is 7x = 293 K. The numerical solution is compared with the analytical solution for Hagen- 
Poiseuille flow. The computed velocity profile shown in Fig. 3.9 fits the parabolic profile of 
the analytical solution. 


Turbulent Pipe Flow 

Turbulent pipe flow can again be used to select the turbulence model for flows through 
channels. As in the case of laminar flow, the hydraulically smooth pipe with a length L/D = 
50 is computed. The Reynolds number is Rep = 1- 10° and the temperature T,, = 293 K, 
with an intensity of turbulence of Tus. = 0.01 used. Figure 3.10 shows a comparison of the 
numerical solutions with Prandtl’s logarithmic wall law. The variation from the analytical 
solution is minimal, indicating that turbulence models which use a logarithmic wall law are 
adequate for pipe flow. 
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Fig. 3.10 : Velocity profile of turbulent Fig. 3.9 : Velocity profile of laminar pipe 
pipe flow Rep = 1- 10° flow Rep = 660 
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Backward-Facing Step 

A critical test case is that of the reattachment of separated turbulent flow in a negative 
pressure gradient. The backward-facing step is a simple geometry which permits investiga- 
tion of reattachment of a flow. Again the numerical results are compared with experimental 
data. The Reynolds number formed with the height of the step H is Rey = 3.7-10*, at a 
temperature of T,, = 293 K. The computed time-averaged velocity profile in the backflow 
region matches the experimentally determined profile well. However it is seen that the reat- 
tachment point varies between «/H = 5.4 and 6.1, depending on the turbulence model used. 
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Fig. 3.11 : Turbulent flow separation behind a backward-facing step 
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In particular, the variations in the turbulent kinetic energy indicate that the isentropy of 
turbulent viscosity, as assumed in the turbulence model, is not satisfied for separated flows. 


Shock Tube 

In order to verify the computation of time-dependent, compressible flows through channels, 
we consider the unsteady shock tube flow. In a shock tube, the high pressure region (driving 
tube) and the low pressure region (running tube) are separated by a membrane. At the time 
t = 0, the membrane is burst and a flow process takes place. A shock wave moves into the 
running tube, and a steady expansion fan moves through the driving tube. The aim in this 
test case is to determine the shock velocity and the analytical Hugoniot values across the 
shock wave. In addition, it is to be shown that the shock computation is independent. of 
the grid size used. The simulation results are compared to the one-dimensional analytical 
solution. 


Figure 3.12 shows a snapshot of the density, pressure, temperature and velocity behind the 
shock wave, at the boundary between the driving gas and the gas in the running tube, and 
in the expansion fan in the driving tube. At the time ¢ = 0 the membrane is at the position 
x/L = 0.5. It can be seen that the changes of state from the shock wave (right), to the 
boundary between the media (center), and on to the expansion fan (left) are represented 
very well. However the discontinuities across the shock wave and the boundary between the 
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Fig. 3.12 : Shock tube flow 
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media are no longer sharp. 


3.3 Industrial Applications 


Now that we have successfully tested the software package in the introductory training 
course, we look at how to approach an industrial project. This second part of the training 
course is to be found at the address 


http: //www-isl.mach.uni-karlsruhe.de/sindustrieprojekt.html 
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Fig. 3.13 : Transonic wing with contour variation (bump) 
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Transonic Wing 


As we have already seen in the introductory Chapter 1.2, the development of the wing of 
an airplane aims to deliver a drag coefficient c, which is a small as possible for a given 
lift coefficient c,. The first step would be to increase the stretch for which the boundary 
layer on the wing is laminar. ‘This leads to transonic laminar wings with a maximum sweep 
angle of about ® = 20° (see Figure 1.26). However for realistic sweep angles of civil aircraft 
the laminizing effect is nullified, because of the appearance of so-called transverse flow 
instabilities in the three-dimensional boundary layer. For this reason other measures are 
needed to reduce the drag. One possibility is the so-called adaptive wing which adapts itself 
optimally to the current state of flight. Another possibility is the bump, a change in the 
contour on the surface of the wing in the shock regime which influences the shock-boundary- 
layer interaction on the wing to such a degree that a drag reduction of up to 9% becomes 
possible. 


Figure 3.13 shows the effect of such a variation in contour on the wing. First of all the 
computational grid on an Airbus A320 model wing is shown. The free stream Mach number 
is M,, = 0.78, the Reynolds number is Re; = 26.6 - 10°, the angle of attack a = 2° and 
the sweep angle ® = 20°. With the variation in contour, the section of the solution in the 
shock-boundary-layer interaction region shows that the shock wave fans out and a wake 
expansion region (shaded in gray) is formed. The change in contour in the shock region and 
the consequent expansion reduces the tendency towards separation. The boundary layer 
thickness is reduced and the fanning out of the shock simultaneously leads to a reduction 
in the drag. These effects both lead towards the desired drag reduction. 


Flow Past an Automobile with Mirrors 


The structure of the wake flow of an automobile has already been introduced in Chapter 
1.2. Of particular interest is the effect of the outer mirrors on the aero-acoustics of the 
automobile. The noise effects due to the flow past the mirrors are locally concentrated on 
the side windows and are carried on the one hand, into the far field and, on the other hand, 
through the side windows and the sealing around the doors into the passenger area. Figure 
3.14 shows the CAD geometry of the automobile as well as a sketch of the wake flow of the 
mirror. Close to the surface of the automobile, the no-slip condition at the wall causes a 
horseshoe vortex to form. Above the horseshoe vortex shear layers are formed at the edge 
of the mirror and these pass over into the backflow behind the mirror and into the wake 
flow downstream from the saddle point. A high level of noise is formed in the shear layers 
of the horseshoe vortex and the wake flow, and this disturbance is heard in the passenger 
area. The numerical calculation of the flow past the mirrors is carried out with a free stream 
velocity of 140 km/h and the Reynolds number Rep = 5- 10° of the outer mirror. If we 
recall the discussion of the flow past a cylinder in Chapter 2.4.5, no dominant separation 
frequency of the Karman vortex street, which could cause an additional sound level at the 
Kármán frequency, occurs at the Reynolds number 5- 10°. The numerical computation and 
the experiments confirm this fact for the semi-cylindrical shape of the mirror. Figure 3.14 
shows the local noise level, computed from the numerical calculation, in two planes vertical 
and horizontal compared to the mirror. These are compared to experimental results. ‘The 
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Fig. 3.14 : Aero-acoustics of an automobile mirror 


measurement of the local sound sources here is carried out using aero-acoustic holography. 
It can easily be seen that high noise levels occur in the shear layers already discussed and 
as sources of the change of shape of the mean flow. These can be seen particularly close to 
the saddle point of the wake flow. 
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Coriolis Revolution Sensor 


In order to be able to control the stability of a vehicle, the Coriolis force is measured using 
micro-motion sensors in EPS systems. The sensor consists of a closed container of height 
240 um, in which fingers with micrometer dimensions oscillate against each other with 
frequency f. This gives rise to an unsteady gap flow whose flow losses are kept as small as 
possible, and so the surrounding pressure of the periodically oscillating fingers is lowered to 
the gas kinetic region. The Reynolds number Rey formed with the height of the gap H is 
so small that the no-slip condition does not hold at the finger wall. 


As the motion sensor is used, gas leakage causes the pressure in the enclosed sensor container 
to increase continually, and the damping characteristics of the motion sensor deteriorates. 


Figure 3.15 shows a snapshot of the geometry and the moving finite-volume computational 
grid for a section of the motion sensor. The gap height H between the oscillating fingers and 
those at rest is 2.7 um. At an oscillation frequency of f = 1450 Hz, the Reynolds number of 
the gap flow is Rep = 0.01 at a pressure of p = 0.01 bar. Figure 3.15 shows isotachic lines in 
different sections of the area of computation. Horizontal close-ups around the leading and 
trailing edges of the central of the three moving fingers are shown at different heights. The 
different sections are at the lower edge of the finger (z/H = 0.1), in the middle (z/H = 0.5), 
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Fig. 3.15 : Gap flow of an oscillating Coriolis revolution sensor 
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and at the upper edge (z/H = 0.9). The three-dimensional flow past the leading edge can 
be seen. The gap flow met in Chapter 2.4.4 is found in the gaps between the fingers at rest 
and those moving. 


The damping D is made up of the parts due to the pressure and friction forces, shown for 
one period T. 
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Fig. 3.16 : Gap flow in a rotary valve pump 
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Rotary Valve Pump 


The mode of operation of a rotary valve pump is based on the fact that the position of a 
sliding piston is changed during rotation, so that the chambers filled with fluid continually 
increase their volume from 0° to 180° as they coat the inlet slit, and then, after passing 
through the inlet slit, continually decrease their volume until the compressed medium exits 
through the outlet slit. As the fluid passes through the compression gap, its pressure and 
temperature rise. 


Figure 3.16 shows the simplified geometry of a rotary valve pump as well as the computa- 
tional grid of the gap flow. It consists of an inlet ring over the entire angular range, and 
the matching outlet element with corresponding angles. The first segment has the starting 
position as shown, whereby the cells are initially switched off at the inlet collar, and then 
successively switched on during rotation. In this type of calculation, the interaction effects 
between neighboring cells in the suction and expulsion phases are not taken into account. 
The grid motion is initially fixed, with one row of cells at the inlet, at one degree per time 
step. In the first steps, the cell layers are successively activated until the entire segment is 
in contact with the inlet. The segment then moves clockwise until one row of cells reaches 
the outlet edge. At this point, the cell layers are successively deactivated. Figure 3.16 shows 
the temperature distributions at six different angular positions of the sliding piston. The 
inlet temperature of the fluid is 293 K and the wall temperature of the pump 333 K. In the 
inlet region, vortex formation takes place in the flow past the piston. In the outlet region a 
jet-like expansion flow is formed in the chamber. The temperature of the compressed flow 
in the gas increases to values of up to 600 K and then drops again greatly in the expansion 
area of the outlet chamber. 


Flow Through an Engine Cylinder 


Another example of application in automobile technology is the flow through an engine 
cylinder. We have already seen the cycle of the moving piston in an Otto engine in Figure 
1.7. Therefore in order to compute the flow through an engine, it is necessary to model the 
following processes physically correctly: the mixing of the air-fuel mixture in the suction 
phase; the production of the highly turbulent kinetic energy in the compression phase; the 
combustion after ignition; the expansion phase; and finally the expulsion of the combustion 
products. The problem in doing this is the turbulence models employed in the software 
packages. 


The suction phase of a 4-valve Otto engine will first be modeled. Here we consider to what 
degree the turbulence models implemented in the commercial software correctly represent 
the mixing of the air-fuel mixture. Ultimately this is a requirement for the computation of 
the compression and combustion in the Otto engine. Figure 3.17 shows a picture of the com- 
putational grid at a certain instant of the piston motion in the suction phase of the motor. 
Here the motor is taken to be mechanically driven with open valves. The computational grid 
is continuously extended as the piston motion causes the volume of the cylinder to increase. 
We assume that the engine is symmetric through the middle plane, so that two inlet valves 
are taken into account. The numerical results show the coherent vortex structures typical to 
suction pipe injection, such as rotating and tumble flows during the suction phase and well 
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into the compression phase. The turbulent mixture of the inlet flows of both valves shown in 
Figure 3.17 (shown dark and light) is, however, not in agreement with experimental results. 


The drive for the turbulent momentum exchange through the mean velocity gradients of 
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Turbulent mixing in the suction phases of a four-valve Otto engine 


Fig. 3.17 
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the shear layers is not computed precisely enough. The distribution of the turbulent kinetic 
energy through turbulent diffusion due to velocity and concentration gradients on the one 
hand, and of the vortex strength concentration in the main vortices on the other hand, was 
predicted to be too small. The numerical simulation does not yield much mixing between 
the fresh gas and the exhaust gas until well into the compression phase. These are effects of 
the isentropy assumed in the turbulence model, which is not present in the rotational inlet 
and compression phases. The development of new turbulence models is necessary in order 
to compute the turbulent mixing more precisely. 


Valve Flow 


The fluid mechanical elements of an ABS system include hydraulic pipes, and the pumps 
and control valves which cause the periodic stuttering brake force of an automobile. In 
what follows we consider the steady three-dimensional flow in a valve in a given state. The 
numerical computation of the force on the valve head is of interest to predict the size of the 
reaction force needed. 


The geometry and the computational grid of the hydraulic valve are shown in Figure 3.18. 
The hydraulic fluid flows from the left past a throttle position and through a ring gap which 
is formed by the moving part of the valve (diaphragm) and the rigid base of the valve. Up 
until this point, the configuration is rotationally symmetric with respect to the axis of the 
free stream bore hole. Because the bore hole for the downstream flow is directed upwards, 
the flow is turned in this direction, so that a three-dimensional separated flow is formed in 
the valve space. 


The computed streamlines show the three-dimensional structure of the valve flow. A vortex 
forms around the valve head leading into the downstream flow pipe. In the front end of the 
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Fig. 3.18 : Structure of the flow in a hydraulic valve 
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valve space a secondary vortex is formed, and in the rear end, an additional vortex which is 
connected with a diagonal separation line on the diaphragm. If we begin the computation 
with a rough grid, a reversal of the orientation of this vortex point takes place as we move 
to the next finest computational grid. Therefore a system involving a total of three vortex 
points, with the corresponding saddle points at the fixed edge, forms along the sides of the 
pipe. The point of separation is in the form of a focus. This can also be represented on the 
roughest grid, and moves downstream with better resolution of the grid. A system of saddle 
points and nodes forms above the diagonal separation line in the finest grid shown. This is 
clearly bounded towards the front and includes the rear part of the diaphragm. 


At first we can see that the computed structure of the flow depends on the computational 
grid chosen. Even for the same grid, different turbulence models lead to different structures of 
this three-dimensional separated flow. Turbulence models assume isentropy in the turbulent 
flow, and this is not true in this flow example. 


This application example was chosen for this textbook for precisely this reason. Even if the 
fluid mechanical software is used physically correctly, we see that the skills and experience 
of an engineer are also required, and these are difficult to teach in a textbook. 


Volume Flux Sensor 


The separation frequency of the Karman vortex street behind a blunt body can be used 
to measure the volume flux of liquids, gases and steam in a pipe flow. The volume flux 
measurement is based on the proportionality between the flow velocity and the separation 
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Fig. 3.19 : Karman vortex street behind a cylinder in a pipe flow 
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frequency of the Kärmän vortex street behind a blunt body, shaped as a cylinder and set 
into the pipe. The cylinder with diameter d in the pipe of diameter D is shown in Figure 
3.19. The diameter ratio D/d is 4. The Reynolds number of the numerical computation 
is related to the time-averaged free stream of the turbulent pipe flow and has the value 
Reg = 5300. 


Figure 3.19 shows vertical and horizontal cross-sections of the normal components of the 
vortex strength vector. In the upper picture we see a snapshot of the Karman vortex street 
in the pipe. Large vortex strengths are found in both the periodically separating vortices 
and the turbulent boundary layer at the edge of the wall. The large vortex strength of the 
approaching wall boundary layer can also be seen in the horizontal cross-section. Because 
of the no-slip condition at the wall, this boundary layer becomes a steady horseshoe vortex 
around the cylinder. The effect of the wall integration into the pipe in the wake of the 
cylinder is seen in a characteristic expansion region of large vortex strength, downstream 
from which the turbulent boundary layer forms again, and the Karman vortex street expands 
undisturbed. The computed Strouhal number of 0.3 agrees with the experimental values. 


Smoke Detector in a Pipe System 


For safety reasons, smoke detectors are installed into the pipe system of an air-conditioning 
system of an airplane. On the one hand they are intended to warn the crew of smoke in 
the airplane, and on the other hand, they activate an automatic valve which ensures that 
the smoke is sucked away, particularly in the cockpit and cabin. In order to be able to 
determine the correct position for such smoke detectors, a steady-state computation of the 
air-conditioning system was carried out. The aim was to compute the degree of saturation 
of the air with smoke particles at each detector position which would be necessary to cause 
the smoke detector to go off. 


t=0.2 





t=0.4 





t=0.6 
computational grid, 200000 grid points 





Fig. 3.20 : Concentration fronts in a smoke detector 
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In order to construct the system safely, two redundant detectors are also integrated into the 
system. Figure 3.20 shows the positions of installation and the surface network. The pipe 
system was widened at the detector positions, in order that no obstruction be caused by 
smoke sensor jutting into the cabin or cockpit. The surface network shows how the detectors 
are linked together. As well as the detector region, a suction pipe was also included into the 
computational region. 


The flow is turbulent at a Reynolds number Rep = 1.9- 10°. The given intensity of turbu- 
lence for the medium at the entrance to the simulation region was Tu,, = 5%, and the pipe 
flow was matched to the cabin air at an altitude of 11.2 km. The volume flux was given 
as V = 118 /s. Starting out from the steady numerical solution, the smoke distribution is 
computed with a passive scalar, i.e. the transport of the inert massless smoke is represented 
based on a steady solution. The particle paths for the massless, frictionless smoke particles 
are integrated along the steady streamlines. Figure 3.20 shows a snapshot of the smoke dis- 
persion at three different points in time in the smoke detector region. The rate of dispersion 


of the smoke shows that the smoke sensors will go off within 1 second. 


Clean Air Conditioning in an Operating Theater 


Pure air covers, such as those used in clean rooms, are used in operating theaters. A laminar 
vertical flow ensures that germ-free air is constantly blown in from above and dirt removed 
from the operating table. Cooled clean air is introduced to the operating area and is sucked 
away below and to the sides. The colder air sinks downwards with a small flow velocity. On 
the one hand this leads to the formation of a shear layer between the normal air in the room 
and the germ-free clean air, and on the other hand to the dirt introduced in the operating 
area being wiped away. The essential quantity considered to determine the effectiveness is 
the concentration of the clean air on the operating table. 


Because the flow is strongly asymmetric, due to the position of the hot operating lights and 
supports and the air suction, a three-dimensional computational grid is generated to cover 
the entire operating theater, including operating table, lights and supports. 
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Fig. 3.21 : Clean air conditioning in an operating theater 
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The clean air is supplied with a Reynolds number of Rep = 6100, where D is the diameter 
of the operating lights, and at a temperature of 291 K. The temperature of the room is 
293 K. 


The warm surfaces of the lights are modeled as isothermal walls whose temperature dis- 
tribution has been determined experimentally. The temperature expansion is computed as 
convective heat transport. 


After the steady numerical solution is computed, a transient expansion calculation is car- 
ried out for the concentration of the germ-free air. This corresponds to the particle path 
integration in the smoke detector example. Figure 3.21 shows a longitudinal cross-section 
of the clean air concentration in the operating theater, after germ-free air of concentration 
1 has been allowed to enter through the inlet. Regions where the clean air concentration is 
higher are light, and those there the air is dirtier are dark. The result of the computation 
shows that, although the surface of the operating table is supplied with clean air, the effect 
on the flow of the operating lights may not be neglected. 
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